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14.  ABSTRACT 

Purpose:  The  objective  of  this  USAMRMC  training  and  researching  project  is:  (1)  to  train  the  PI  in  the  optics,  biology  and 
biochemistry  fields  towards  a  multi-disciplinary  prostate  cancer  researcher;  and  (2)  to  understand  and  investigate  the  receptor  target 
contrast  agents  and  energy  transfer  from  contrast  agents  to  prostate  cancer  cells  using  ultrafast  time -resolved  spectroscopy. 

Specific  tasks  for  the  report  year  (04/15/201 1  -  04/14/2012):  (1)  Training  on  biological  knowledge  and  technique  for  cell  culture  and 
tissue  studies  in  Memorial  Sloan-Kettering  Cancer  Center  (Mons.1-6);  (2)  Training  on  time-resolved  polarization  and  steady-state 
spectroscopic  measurements  for  contrast  agents  using  streak  camera  and  the  femtosecond  Ti:Sapphir  laser  system  (Mons.7-12);  and  (3) 
Following  the  training,  applying  Pi’s  optical  knowledge  into  the  training  program  to  generate  researching  outcomes  (Mons.1-12). 

Major  accomplishments  and  findings:  (1)  PI  has  been  trained  on  biological  knowledge  and  technique,  and  applied  for  IRB  approval 
for  prostate  cells  and  tissues  measurements;  (2)  PI  has  been  trained  on  time-resolved  polarization  and  steady-state  spectroscopic 
measurements  for  contrast  agents  using  streak  camera  and  a  femtosecond  Ti:Sapphir  laser  system  to  explore  time-resolved  dynamics  of 
fluorescein  dye-labeled  polymer  chains  with  different  molecular  volume;  (3)  Followed  the  training  program,  PI  has  applied  his  optical 
knowledge  to  the  following  research:  (3a)  use  a  prostate  cancer  receptor-targeting  contrast  agent,  namely  Cytate,  to  improve  the 
detection  depth  of  the  cancerous  prostate  areas  using  optical  scanning  imaging  technology  study  to  improve  the  detection  depth  of  the 
cancerous  areas;  (3b)  study  the  fractal  dimension  parameters,  absorption  and  scattering  coefficients  for  cancerous  and  normal  prostate 
tissues  systematically;  and  (3c)  evaluate  the  chemotherapeutic  effects  on  malignant  cells  using  their  nature  fluorescence. 

15.  SUBJECT  TERMS 

Prostate  cancer,  cell  culture,  steady-state  and  time-resolved  spectroscopy.  NIR  scanning  polarization  imaging,  independent  component 
analysis,  absorption  and  emission  of  cancerous  and  normal  prostate  tissues,  and  contrast  agents. 
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Annual  Report 

For  Army  Grant  W81XWH-1 1-1-0335 
(CUNY  RF  47204-00-01) 

Period  covered  04/15/11  -  04/14/12 

Introduction 


The  objective  of  this  Army  Medical  Research  and  Materiel  Command  research  project  (# 
W81XWH-1 1-1-0335)  is  to  train  the  PI  in  the  optics,  biology  and  biochemistry  fields  towards  a 
multi-disciplinary  prostate  cancer  researcher.  The  proposed  training  plan  will  also  establish  a 
prostate  cancer  research  program  at  City  College  of  New  York  (CCNY)  through  close 
collaboration  with  Memorial  Sloan-Kettering  Cancer  Center  (MSKCC)  and  Washington 
University  School  of  Medicine  at  St  Louis  (WUSM).  The  focus  of  the  training  component  of  the 
project  is  to  familiarize  the  PI  to  the  biological  aspects  of  cancer  research  through  attending 
relevant  courses,  and  cancer  research  practicum  through  laboratory  rotations.  The  objectives  of 
the  research  component  of  the  project  are  to  develop  optical  approaches  to  (1)  validate  the 
potential  over-expressed  receptors  on  prostate  cancer  cells;  and  (2)  study  the  energy  transfer 
from  the  QD-based  contrast  agents  to  the  cancer  cells.  While  the  emphasis  of  the  program  is  on 
the  training  component  during  the  first  year,  significant  progress  has  been  made  in  both  training 
and  research  during  the  first  reporting  period  (April  15,  2011  -  April  14,  2012)  covered  by  this 
report. 


The  tasks  for  the  first  year  (04/15/2011  -  04/14/2012)  were:  “Task  1,  Training  of  biological 
knowledge  and  technique,  and  apply  for  IRB  approval  for  prostate  cells  and  tissues 
measurements  (Months  1-6)”;  and  “Task  2,  Training  of  biochemistry  knowledge  and  techniques 
(Months  6-12).” 


The  tasks  for  the  first  year  were  performed  successfully  and  significant  progress  was  achieved 
as: 


1.  Accomplishment  of  the  planned  training  program  on  biological  knowledge  and  technique 
through  laboratory  rotations  and  attending  group  meetings  and  seminars  at  MSKCC  [1]; 

2.  Accomplishment  of  the  planned  training  program  on  time-resolved  polarization  and  steady- 
state  spectroscopic  techniques  through  -  hands-on  laboratory  training  at  Institute  for  the 
Ultrafast  Spectroscopy  and  Lasers  (IUSL)  in  CCNY,  and  development  of  research  on  time- 
resolved  polarization  spectroscopic  model  to  explore  ultrafast  dynamics  of  fluorescein  dye- 
labeled  polymer  chains  with  different  molecular  volume  in  both  experiment  and  theory  [2]; 

3.  Starting  of  the  planned  training  program  on  synthesis  of  receptor-targeting  contrast  agents 
through  laboratory  rotations  and  hands-on  laboratory  training  at  WUSM; 

4.  Development  of  research  on  non-  or  less-  invasive  optical  near  infrared  scanning  optical 
imaging  approach  for  prostate  cancer  detection  using  receptor- targeting  contrast  agents  [3]; 

5.  Development  of  research  on  systematic  studies  of  fractal  dimension  parameters,  and 
absorption  and  scattering  coefficients  for  cancerous  and  normal  prostate  tissues  [4];  and 

6.  Development  of  research  on  spectroscopic  approaches  to  evaluate  the  chemotherapeutic 
effects  on  malignant  cells  [5], 

Report  Body 

Since  this  grant  is  training-concentrated  and  budget  of  this  training  grant  is  limited,  the  training 

program  is  also  jointly  supported  by  other  related  grants  numbered  W81XWH-08-1-0717  (for 
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which,  Mentor:  Professor  Alfano  is  Co-PI)  and  W81XWH-07- 1-0454  (for  which,  Co-mentor: 
Prof.  J.  Koutcher  is  PI),  in  training  facilities  such  as  equipments,  instruments,  biological  and 
biochemical  materials  and  supplies  (cells,  dyes,  cultural  medium  etc.).  During  the  training,  the 
PI,  as  a  trainee,  gains  the  experiences  and  knowledge  in  biology,  biochemistry  and  biomedical 
optics.  Furthermore,  the  PI,  as  a  postdoctoral  research  fellow,  develops  programs  to  achieve 
researching  purposes  under  the  supports  of  these  three  grants.  The  following  highlights  our  major 
accomplishments  during  the  first  year: 

(1)  Laboratory  cell  culture  training  and  animal  handling  training 

The  activities  reported  in  this  section  fulfilled  Task  1  described  in  the  proposed  S.O.W  for 
training  program  on  biological  knowledge  and  technique. 

During  the  first  year,  the  PI  (trainee)  has  attended  the  following  biological  research  training  and 
small  animal  training  sessions  ( Training  Task  1),  held  by  Research  Animal  Resource  Center 
(RARC)  at  MSKCC: 

•  RARC  Orientation  training  session  (Certificate  as  attachment  1); 

•  RARC  Basic  Mouse  training  session  that  included  handling  and  restraint;  injectable  anesthetics 
and  anesthetic  monitoring;  parenteral  administration;  and  euthanasia  (Certificate  as  attachment 
2); 

•  RARC  Rodent  Survival  Surgery  training  session  that  involved  surgical  equipment  and 
instruments,  disinfection  and  sterilization,  animal  preparation,  surgeon  preparation,  suture 
materials  and  wound  closure,  and  post-operative  care  (Certificate  as  attachment  3). 

The  trainee  rotated  through  the  following  laboratories  and  received  training  on  cell  culture,  flow 
cytometry,  fluorescence  microscopic  imaging  techniques,  bioluminescence  assay  application, 
cells  viability  monitoring,  and  other  key  techniques  used  in  cancer  research. 

Laboratory  PI:  Prof.  Jason  Koutcher  (role:  Co-mentor) 

Key  person  in  Laboratory  training:  Dr.  Ellen  Ackerstaff  (role:  the  trainer),  and  Ms.  Natalia 
Kruchevsky  (role:  the  trainer,  Laboratory  technician) 

Laboratory  background:  One  focused  goal  of  the  current  research  in  Prof.  Kouthcer’s  lab  is  on 
distinguishing  metastatic  from  non-metastatic  clinical  breast  cancer  using  non-invasive 
techniques. 

Training  Topic:  Cell  culture,  flow  cytometry,  fluorescence  microscopic  imaging  techniques, 
bioluminescence  assay  application,  and  cells  viability 

•  Cell  culture :  The  trainee  learned  the  basic  cell  culture  techniques.  The  trainee  had  the 
responsibility  to  maintain  one  metastatic  prostate  cancer  cell  line:  PC3  in  training  program  and 
some  other  kinds  of  cancer  cell  lines:  4T1  [1],  HTB22  and  HTB126  [5].  Routine  maintenance, 
such  as  subculture  and  cell  counting,  was  performed. 

•  Flow  cytometry  technique,  fluorescence  microscopic  imaging  technique,  and  bioluminescence 
assay  application:  The  trainee  learned  to  use  bioluminescence  assay,  flow  cytometry  technique, 
and  use  fluorescence  and  conventional  microscopic  imaging  techniques  to  monitor  and  analyze 
cell  life  cycle. 
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•  Cells  viability.  The  trainee  learned  how  cell  metabolically  adapts  to  the  dynamic 
microenvironmental  changes,  such  as  major  nutrient  supply  vs.  deprivation  (glucose  and 
glutamine)  and  oxygen  availablilty  vs.  hypoxia  [1]. 

In  this  research  and  training  program,  4T1  (metastatic  cell)  and  67NR  (non-metastatic  cell) 
cancer  cell  lines  were  studied  in  an  NMR-compatible  cell  perfusion  system  and  investigated  their 
difference  in  dynamic  metabolic  response  to:  (1)  oxygen  availability  vs.  hypoxia;  and  (2) 
glutamine  (Gin)  supply  vs.  Gin  deprivation. 

In  this  study,  4T1  and  67NR  cells 
were  cultured  and  studied  to 
determine  the  response  of  cell 
growth  to  modified  cultural  media 
composition  (DME,  supplemented 
with  1%  P/S  and  10%  FBS: 

DMEcompi).  (1)  DMEeompi  contains 
25  nM  glucose  (Glc)  and  6  mM 
Gin;  (2)  DME  with  2  mM  Gin;  (3) 

DME  without  Gin;  and  (4)  DME 
with  Glc.  Cells  were  incubated  for 
48  hours  in  each  medium,  harvest, 
and  counted  (ViaCount  Assay, 

Guava  Technologies).  Three 
independent  experiments,  with 
duplicate  samples  per  experiment, 
were  performed  for  each  cell  line. 

As  shown  in  Fig.  1,  the  in  vitro 
cell  growth  is  dependent  on  the  availability  of  both  Glc  and  Gin,  in  4T 1  and  67NR  lines:  cell 
growth  drop  by  85-89%  in  the  absence  of  either  Gin  and  Glc,  and  remains  at  35%  -  45%  when 
Gin  alone  is  decreased  by  2/3;  no  significant  difference  were  observed  in  growth  between  the 
two  cell  lines  in  response  to  each  nutrient  deprivation. 

Although  no  significant  difference  in  cell  growth  were  observed  in  vitro  between  4T1  and  67NR 
cells,  when  exposed  to  the  defined  metabolic  stress  conditions,  each  cell  line  adopts  different 
metabolic  strategies  to  cope  with  said  changes  depending  on  the  level  of  oxygen  availability. 
Mostly,  non-metastatic  67NR  cells  reveal  higher  glycolytic  activity  than  metastatic  4T1  cells. 

(2)  Laboratory  training  of  time-resolved  polarization  spectroscopy  of  contrast  agents 

The  activities  reported  in  this  section  fulfilled  Task  3  described  in  the  proposed  S.O.W  for 
training  program  on  time-resolved  polarization  and  steady-state  spectroscopic  study  of  contrast 
agents. 

During  the  first  year,  the  PI  (trainee)  has  accomplished  hands-on  laboratory  training  of  the 
ultrafast  spectroscopy  and  lasers  techniques  to  investigate  the  rotational  motion  of  fluorescent 
molecules  in  solution.  In  this  study,  time-resolved  fluorescence  spectroscopy  of  Indocyanine 
Green  (ICG)  and  fluorescein  dye-labeled  polymer  chains  with  molecular  volume  ranged  from 
-410  angstrom  to  -5000  angstrom  at  parallel  and  perpendicular  polarization  directions.  A 
theoretical  fluorescence  depolarization  model  developed  by  G.  Weber  [6]  was  extended  from 
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micro-second  to  picosecond  regimes  to  investigate  and  demonstrate  the  origin  of  the  different 
decay  rate  of  the  two  fluorescence  components.  Experimental  tests  of  this  theory  were  verified 
by  using  streak  camera  optical  multi-channel  detections  for  the  time -resolved  parallel  and 
perpendicular  components  of  the  fluorescence  emission  [2]. 

We  assume  dye  molecules  placed  at  the  center  and  along  the  x-axis  of  the  co-ordinates  excited  by 
light  pulse  of  an  arbitrary  linear  polarization  propagating  at  the  YO  direction  at  time  to.  G.  Weber 
has  pioneered  the  theory  of  the  fluorescence  polarization  from  an  isotropic  molecular  system  due 
to  molecular  rotation  in  nanosecond  range  [6].  By  extending  this  theory  from  nanosecond  to 
picoseconds  [7],  the  total,  parallel  and  perpendicular  intensity  of  the  fluorescence  excited  by  a 
linear  polarized  light  is  obtained  as: 

I(t)  =  I0ex p( — —)  l  (t)  =  —  exp(-— )(l  +  2r0exp( — — )|  G(?)  =  xrexP(  )  l-r0exp(  )  (1) 

Tr  3  Tf  {  Trol  J  3  ^  Tr°>  ) 

The  average  decay  time  of  the  polarized  components  (t>u  and  (t)j_  are: 


Wrjr  ro,l(Trot  +  TF)\- 
!-r0[rro(  l(jrol+TF)] 


+rF)]2 

l  +  2r0  Vrot  /(T„,+^f)] 


(2) 


where  Tf  is  fluorescence  lifetime  of  the  molecule,  Trot  is  rotation  time  of  the  molecule,  and  ro  is 
initial  fluorescence  anisotropy  at  time  to.  The  time  dependent  fluorescence  anisotropy  is  usually 
used  to  describe  fluorescence  depolarization  property  using  the  equation  of  [8,  9]: 


r(t) 


(t) 

I\\(t)  +  2I±(t) 


(3) 


where  the  initial  anisotropy  indicates  the  initial  orientation  distribution  of  the  fluorescent  dipoles 
and  is  reported  to  be  0.4  theoretically  [6-9]  and  the  rotation  time  t,-o1  can  be  expressed  in  terms  of 
the  diffusion  coefficient  (D),  or  the  solvent  viscosity  (q)  and  the  molecular  radius  (a)  by: 


xrot  =  (6D)_1 


47iqa3 

3kT 


(4) 


where  k  is  the  Boltzmann  constant  and  T  is  the  absolute  temperature  [6-9]. 

In  the  time-resolved  fluorescence  experiments  for  ICG,  optical  pulses  of  130  fs  at  800  nm  (at 
400  nm  for  fluorescein  derivatives  measurements)  with  a  repetition  rate  of  82  MHz  from  a 
Coherent  Mira  900  Mode-locked  Ti: Sapphire  laser  were  used  to  pump  the  samples  of  ICG  and 
fluorescein  dye-labeled  polymer  chains  [2,  7].  The  temporal  profiles  of  the  fluorescence  emitted 
from  ICG  solution  for  two  polarization  directions,  parallel  and  perpendicular  to  the  polarization 
of  the  excitation  at  800  nm  were  measured  and  shown  in  Fig.  2(a).  There  is  a  distinct  difference 
between  the  decay  slopes  of  the  two  components,  as  shown  by  the  dashed  and  thin-solid  line  in 
Fig.  2(a),  respectively.  In  order  to  see  clearly,  the  peaks  of  the  perpendicular  and  the  parallel 
components  were  normalized  to  the  total  intensity  (which  was  obtained  by  using  equation 
/(t)=/±(t)+2/||(t)  shown  as  the  thin-dot  line  in  Fig.  2(a)).  The  thin  line  in  Fig.  2(b)  displays  the 
time  evolution  of  polarization  anisotropy  r(t).  The  obvious  phenomenon  is  the  observation  of 
depolarization  as  mentioned  by  the  previous  researchers  [8,  9].  Here  we  try  to  explore  the 
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understanding  in  first  time  based  on  rotation,  initial  orientation  and  lifetime  of  the  excited 
molecules. 


The  fluorescence  lifetime  Tp  is  acquired  from  fitting  the  total  temporal  fluorescence  intensity  data 
shown  in  Fig.  2(b)  using  equation  (1).  The  polarization  anisotropy  peak  value  ro  and  the  rotation 
time  Trot  of  ICG  molecule  in  solution  can  be  obtained  by  fitting  the  experimental  data  of  r(t)  in 
Fig.  2  using  equation  (2). 

Time-resolved  fluorescence  of  ICG  solution  Time  Resolved  Polarization  Anisotropy  of  ICG  solution  The  rotation  effect  for  parallel  and  parpcndicular  components 
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Fig.  2(a)  Profiles  of  the  time-resolved  emission  jj  and  _L  compoi  its;  (b)  Time-dependent  polarization 
anisotropy;  (c)  The  schematic  temporal  contribution  of  rotation  to  ||  and  _L  component.  All  experimental 
data  are  shown  as  thin  line,  but  calculation  as  thick-line. 

The  dimensions  of  ICG  molecules  were  estimated  as  oblate  with  semi-axes  as  13A  and  2 A  from 
bond  length  and  van  der  Waals  radii.  The  theoretically  rotation  time  of  ICG  is  calculated  using 
equation  (4)  and  these  values  of  q=1.003xl0'3  Pa-s  (viscosity  of  water),  T=293  (20°C  -  the  room 
temperature),  and  the  Boltzmann  constant  k«  1 .380  x  10'  .  Using  equation  (2),  the  theoretical 
decay  times  of  Ij_(t),  I(t)  and  I||(t)  can  be  quantified.  The  experimental  data  and  calculated  results 
are  concluded  in  Table  I. 


Table  I:  Comparison  of  the  experimental  and  theoretical  data  of  ICG 


Type  of  Data 

ro 

^rot 

<t)|| 

tf 

<t>± 

Experimental 

0.391T0.0238 

353.37+10.55 

169.35+4.58 

192.46±5.28 

239.15+6.48 

Theoretical 

0.4 

351.02 

163.79 

— - 

244.91 

To  illustrate  the  temporal  contribution  of  the  molecular  rotation  to  ||  and  _L  components 
schematically,  Fig.  2(c)  is  generated  by  using  31  |(t)/I(t)  and  3I_i_(t)/I(t)  with  equation  (1)  showing 
the  influence  of  the  rotational  component  for  the  change  of  Ij_(t)  and  I||(t).  The  trend  of  the 
influence  of  the  rotation  to  Ij_(t)  evidently  rises  up  while  that  to  I||(t)  drops.  Even  I  (t)  and  I||(t) 
emitted  from  same  type  of  molecules  with  same  lifetime,  the  influence  of  the  molecules’ 
Brownian  rotation  causes  decay  of  parallel  component  faster  than  the  total  decay  and  the 
perpendicular  component  when  the  lifetime  and  the  rotation  time  of  the  molecule  are  on  same 
order  of  magnitude.  Our  study  extends  Web’s  theory  to  picoseconds  range. 


The  time-resolved  fluorescence  polarization  studies  were  also  performed  for  fluorescein  dye,  and 
the  fluorescein  dye-labeled  dextrans  with  different  Molecular  Volume  (MV)  from  400  to  5,000 
angstrom  .  The  dimensions  of  fluorescein  were  estimated  as  an  oblate  ellipsoid  with  semi-axes  of 
7  and  2  angstroms,  respectively  [10].  As  an  example,  Figs.  3  (a)  and  3  (b)  show  the  measured 
temporal  profiles  of  parallel  and 
perpendicular  components  of  the 
emission  from  the  conjugated  dye- 
polymer  with  different  MV  of 
-707  and  -3,200  angstrom’, 
respectively.  The  solid  line  profiles 
show  the  parallel  component,  and 
the  dot  lines  show  the 
perpendicular  one.  Fig.  3  (a)  and  3 
(b)  show  that  I||(t)  is  greater  than 
Ii(t)  at  all  of  the  decay  times 
indicating  that  this  dye-polymer 
keeps  property  of  the  polarization 
preservation.  The  time-resolved 
polarization  anisotropy  r(t)  can  be 
obtained  by  substituting  the  data  of 
I|](t)  and  Ij_(t)  shown  in  Figs.  3  (a) 
and  3  (b)  into  the  time  dependent 
equation  (3).  The  calculated  results 
are  shown  in  Figs.  3  (c)  and  3  (d) 
displaying  the  temporal  profiles  of 
polarization  anisotropy  for 
Fluorescein-polymer  with  different 
MV  of  -707  and  -3,200 

■j 

angstrom  ,  respectively. 

The  salient  features  are:  (1)  the  temporal  polarization  anisotropy  shown  in  Fig.  3(c)  exhibits  a 
faster  decay  than  the  profile  shown  in  Fig.  3  (d),  indicating  that  a  smaller  molecular  volume 
gives  faster  Brownian  rotation  as  described  by  Equation  (4);  (2)  The  initial  anisotropy  r0  shown 
in  Fig.  3  (c)  is  greater  than  that  shown  in  Fig.  3  (d). 

The  rotation  times  can  be  obtained  by  fitting  the  experimental  data  of  r(t)  with  Equation  (3)  and 
the  results  for  a  set  of  Fluorescein-polymer  conjugates  with  attached  chain  having  the  volume 
varying  from  -400  to  5,000  angstrom3  are  displayed  as  solid  lines  in  Fig.  4.  The  theoretical 
values  of  the  rotation  times  with  different  attached  polymer  length  were  calculated  using 
Equation  (4)  and  are  shown  as  the  dot-line  profile  in  Fig.  4  in  comparison  with  the  experimental 
values.  Since  the  fluorescein  dye  appears  as  approximated  oblate  symmetric  rotors  [10],  the 
volume  was  calculated  to  be  -  410  angstrom  ’  as  an  oblate  ellipsoid  with  semi-axes  of  7  and  2 
angstroms,  respectively  while  no  polymer  attached. 


Time-resolved  fluorescence  intensity  Time-resolved  fluorescence  intensity 


Time-resolved  polarization  anisotropy  Time-resolved  polarization  anisotropy 


Time  (Ps)  Time  (Ps) 

Fig.  3.  The  time  dependent  fluorescence  intensity  profdes  of 
parallel  (solid  line)  and  perpendicular  (dot  line)  components  of 
emission  from  conjugated  fluorescein-polymer  of  MV  at  (a)  -707 
and  (b)  -3,200  angstrom3.  The  temporal  profdes  of  the  polarization 
anisotropy  of  conjugated  fluorescein-polymer  of  MV  at  (c)  -707 
and  (d)  -3,200  angstrom3. 


9 


Comparison  of  experimental  and  calculated  rotation  time 


The  salient  features  of  Fig.  4  are:  (1)  the 
rotation  times  raise  with  the  volume  of  the 
attached  polymer  chain,  which  are  in  good 
agreement  with  experimental  values  when  the 
volume  of  conjugated  Fluorescein-polymers  is 
less  than  3,200  angstrom3;  (2)  The  rotation  time 
becomes  flat  when  the  volume  of  the  dye- 
polymer  increases  beyond  -3,200  angstrom3. 

This  indicates  that  for  large  size  molecules, 
where  the  volume  exceeds  a  critical  value,  the 
rotation  time  changes  very  slow  with  volume  of 
the  fluorescent  molecules.  The  incorrect 
prediction  of  Equation  (2)  may  be  because  the 
Fluorescein-polymer  can  no  longer  be 
estimated  as  a  simple  dimension  as  a  sphere  or 
an  oblate  since  it  is  difficult  to  control  the 
polymer  chain  extending  direction.  The 
irregular  shapes  of  molecules  result  in  the  different  rotation  rate  of  Dy  [6,  8,  9],  and  therefore  the 
depolarization  will  not  show  simple  exponential  decay  as  described  by  time -resolved 
polarization  anisotropy  r(t)  in  Equation  (3).  In  addition,  the  rapid  internal  motions  originated  by 
flexible  structure  of  the  conjugated  fluorescein-polymers  cause  the  “extra”  depolarization  other 
than  rotation  of  the  molecule  [8]. 


Fig.  4.  (a)  The  rotation  time  obtained  by  theoretical 
calculation  (dot  line)  and  by  fitting  the  experimental  data 
(solid  line)  of  the  fluorescein  dye-polymer  conjugate  as  a 
function  of  molecular  volume. 


(3)  Three  dimensional  location  of  cancerous  prostate  tissue  marked  with  a  receptor- 
targeted  contrast  agent  (Cytate)  using  near  infrared  scanning  imaging 


The  activities  reported  in  this  section 
partially  fulfilled  Task  2  described  in  the 
proposed  S.O.W  for  training  program  on 
biochemistry  knowledge  and  techniques  of 
receptor-targeted  contrast  agents. 

Cypate-Octreote  Peptide  Analogue 
Conjugate  (Cytate),  a  small  ICG- 
derivative  dye-peptide,  could  be  used  for 
effective  targeting  somatostatin  receptor- 
rich  tumor  cells  in  the  animal  model  [11] 
and  over-expressed  somatostatin  receptors 
on  prostate  cells  [12]  in  human  prostate 
tissue  [13]  studied  using  optical  imaging 
techniques.  To  enhance  the  imaging 
quality,  a  fast  algorithm,  namely  OPTICA 
(optical  imaging  using  independent 
component  analysis),  was  used  to 
characterize  the  intensity  spatial 
distribution  of  Cytate-stained  cancerous 
and  normal  prostate  tissues  embedded  in 


Fig.  5.  A  schematic  diagram  of  the  experimental  arrangement 
used  for  imaging  a  pair  of  Cytate-stained  cancerous  and  normal 
prostate  tissues  embedded  in  large  pieces  of  host  normal 
prostate  tissues.  L;  lense,  NF:  830  nm  narrow  band  filter. 


10 


un-stained  normal  prostate  tissue  [3].  OPTICA  extracts  the  spatial  distribution  of  light  “emitted” 
from  each  individual  target  by  treating  them  as  independent  virtue  sources  and  unmixes  them 
based  on  their  independence.  OPTICA  hence  highlights  the  difference  of  the  optical  properties 
between  different  targets  embedded  inside  the  turbid  medium.  OPTICA  can  be  used  to  determine 
the  existence  and  3D-location  of  inhomogeneities  hidden  in  the  homogeneous  host  medium  [3], 

Fig.  5  shows  a  schematic  diagram  of  the  experimental  setup.  The  sample  is  illuminated  by  a 
collimated  laser  diode  with  k=635nm  in  the  direction  close  to  the  normal  to  the  surface.  Two 
galvanometric  mirrors  are  used  to  scan  the  illuminating  beam  along  the  x-  and  y-directions  on 
the  front  surface  of  the  sample.  A  narrow  band-pass  filter  at  830  nm  was  placed  in  front  of  a 
CCD  camera  to  record  emission  images  based  on  the  emission  peak  of  Cytate  at  ~837nm  [13]. 
When  using  OPTICA,  the  leading  eigenvalues  provide  the  number  and  strength  of  the  targets.  In 
the  sample  preparation,  a  -  pair  of  small  pieces  of  cancerous  and  normal  prostate  tissues  were 
soaked  in  the  same  ~3.2  iiM  concentration  Cytate  solution  for  ~1 5  minutes,  and  then  put  into 
sodium  phosphate  buffer  to  wash  off  the  unbound  Cytate  [13].  Finally,  they  were  covered  by  a 
large  piece  of  normal  prostate  tissue  at  the  depth  of  ~3.6  mm  for  scanning  imaging. 


The  OPTICA-generated  independent  components  using  data  of  the  scanning  fluorescent  images 
of  Cytate-stained  cancerous  and  normal  prostate  tissues  embedded  in  large  normal  prostate  tissue 
are  shown  in  Figs.  6(a)  to  6(f).  Figs.  6  (a)  and  6  (b)  show  the  first  leading  IC  standing  for  Cytate- 
stained  cancerous  tissue  location;  6  (c)  and  6  (d)  show  the  second  IC  indicating  Cytate-stained 
normal  tissue  location;  6  (e)  and  6  (f)  show  residual  -  noise.  The  3D  locations  of  the  Cytate- 
stained  cancerous  prostate  tissues  were  obtained  by  numerical  matching  the  target  to  the  surface 
of  the  medium  until  matching  the  retrieved  IC  [3]. 


plane  for  the  sample  consisting  of  the  Cytate-stained  cancerous  and  normal  prostate  tissues 
embedded  in  un-stained  normal  prostate  tissue,  (a)  and  (b)  show  the  first  leading  IC  indicating  for 
Cytate-stained  cancerous  prostate  tissue;  (c)  and  (d)  show  the  second  leading  IC  standing  for 
Cytate-stained  normal  prostate  tissue;  and  (e)  and  (f)  show  the  residual  (noise)  component. 
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OPTIC A-generated  results  are  summarized  in  Table  II  under  different  conditions  using 
backscattering  polarized  imaging  and  OPTICA,  which  lists  the  OPTIC A-determined  positions  of 
the  different  objects  in  comparison  with  their  known  3D  locations  and  the  contrast  enhancement 
between  cancerous  and  normal  tissues. 


Table  II:  Comparison  of  the  known  and  OPTICA  determined  positions  of  embedded 
objects  and  the  contrast  enhancement  between  cancerous  and  normal  tissues. 


object 

covered  tissue 

known  (x,  y,  z) 
position  (mm) 

OPTICA-generated  (at, 
y,  z)  position  (mm) 

Ic/I„ 

stained  cancerous  tissue 

normal  prostate 

(17.8, 14.6,  2.5) 

(17.9, 14.4,  2.3) 

~  3.9 

stained  cancerous  tissue 

normal  prostate 

(18,  24,  3.5) 

(17.6,  23.8,  3.6) 

~  2.3 

(4)  Systematic  studies  of  fractal  dimension  parameters,  and  absorption  and  scattering 
coefficients  for  cancerous  and  normal  prostate  tissues 


A  systematic  investigation  of  the  absorption  coefficient  (jua),  the  scattering  coefficient  (/./,),  the 
reduced  scattering  coefficient  and  the  anisotropy  factor  (g)  of  tissue  is  necessary  for 

accurate  modeling  of  light  transport  in  tissues.  The  differences  in  the  values  of  jua,  jus,  /r' and  g 
between  normal  and  cancerous  tissues  may  provide  a  basis  for  biomedical  imaging  and  potential 
diagnostics  using  optical  techniques  [4].  In  this  study,  we  studied  jua,  /us,  f/s  and  g  of  cancerous 

and  nonnal  prostate  tissues  in  the  750  -  860  nm  spectral  range  using  extinction  measurements  on 
thin  samples  and  diffuse  reflection  measurements  on  thick  samples.  Light  extinction  by  a  thin 
specimen  (thickness  ((  mean  scattering  length,  ls )  within  the  spectral  range  investigated  here  is 
dominated  by  light  scattering.  Light  transport  through  a  thick  tissue  slice  is  usually  modeled  as  a 
diffusion  process,  and  diffuse  reflectance  measurements  provide  information  for  assessment  of 
jua,  and  n's  =  jus  (1  -  g)  ,  where  g  is  the  anisotropy  factor  [4]. 


Since  jua  (( jus  in  the  750  -  860  nm  spectral  range  for  prostate  tissue,  the  extinction  coefficient,  //, 
=  /us  +  /Lia  ~  jLis,  and  the  optical  density  (O.D.)  spectrum  of  a  thin  prostate  tissue  sample  is 
dominated  by  scattering.  The  scattering  coefficient  ( //v)  is  given  by: 


As  =M,=  2-303 


O.D. 

L 


(5) 


where  L  is  the  light  path  length  in  the  tissue. 

The  backscattered  diffuse  intensity  at  position  ( x ,  y,  0)  on  the  surface  of  the  thick  sample 
modeled  by  diffusion  approximation  is  given  by  [14]: 

I(x,y)  =  /,  C°S3^(1  +  Heff  A) exP(~/V A)  + 1‘  C°l  2Ze  (1  +  74//  A) exP ( - 74// Pi )  (6) 

^nP\  Q7Tp2 

where px  =  sin p)2  +  y2  +  f; cos2 /? ,  p2  =  yj(x-l, sin/?)2  +y2  +  (2ze  +/, cos /?)2  and peff  =  pp~JTt  ; 

lt  =  ( // )"'  is  the  transport  mean  free  path;  ze  is  the  extrapolation  length  that  depends  on  the 
mismatch  of  the  index  of  refraction  at  the  interface  [15].  The  extrapolation  length  is  computed  to 
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be  ze=l.69lt.  Equation  (6)  is  used  to  fit  the  measured  diffuse  reflectance  along  a  line  in  the  x 
direction  on  the  tissue  surface  to  obtain  jua  and  /,,  which  is  1  /«' .  The  anisotropy  factor  g  is  then 

obtained  using  g  =  1  -  (ju's/jus  )  . 

As  many  biological  tissues  have  fractal-like  organization  and  are  statistically  self-similar,  various 
tissue  light  scattering  models  based  on  fractal  geometry  have  been  proposed.  The  distribution  of 
the  number  density  of  particles  of  diameter  dl  is  given  by  [16,  17]: 


n(di)  = 


nd]  /  6 


^Ldf’ 

n 


(V) 


where  Df  is  the  fractal  volume  dimension,  rj0  is  a  scale-dependent  constant,  and  n(dt) Ad  gives 
the  number  density  of  the  particles  with  diameters  between  dt  and  <7(  +  Ad  .  The  value  of  rj0  can 
be  determined  from  the  total  volume  fraction  of  all  particles,  Tv,  using  equations 
7],  =  //(,  [  d*  Dfddi  =>  rj0  =  T  /  f  d]  Dfdd:  [16],  where  dmm  and  dmax,  are  diameters  of  the 

Jdmin  Jdn Jin 

smallest  and  largest  scattering  particles,  respectively. 

In  such  a  fractal  soft  tissue  model,  the  correlation  function  R(r )  of  the  random  fluctuation 
of  the  refractive  index  is  a  weighted  average  of  fluctuations  over  a  range  of  correlation  length  /, 
that  depends  on  the  fractal  dimension  Df  and  the  cutoff  length  or  the  maximum  correlation 
length,  dmax  [17].  Light  scattering  property  of  the  tissue  is  fully  determined  by  the  power 
spectrum  (the  Fourier  transform)  of  R(r).  This  model  has  been  shown  to  be  in  good  agreement 
with  NIR  light  scattering  experimental  results  in  biological  tissues  and  cell  suspensions  [17]. 


The  amplitude  scattering  function  S(0)  at  the  scattering  angle  $  is  given  [17]: 

|5(^)|2  =  fW“ - £^°kDI  'X  —dx 

I  'I  ihu  241  +  2(1-  cos  0)x2]2 


(8) 


2  4  2 

where  4=4n0  8m  ,  k=2nno//.,  no  is  the  index  of  refraction  of  the  surrounding  medium,  and  8m  is 
the  amplitude  of  the  relative  refractive  index  fluctuation.  The  anisotropy  factor  (g),  the  mean 
cosine  of  the  scattering  angle,  was  then  found  as  [17]: 


g  = 


(1  +  cos2 


<9)cos^S(<9)|2  /,  (1  +  cos2  9)S{6)2 

- ? — ! - —  \dQ. - ~ - 

2k2  /  J  2k2 


(9) 


where  £2  is  the  solid  angle  and  dD.=2Ti£\n6d6  for  a  scattering  angle  6  [17].  The  reduced  scattering 
coefficient  in  the  fractal  soft  tissue  model  follows  a  power  law: 


t  Rh—Df 
|IsOC/l 


(10) 


This  provides  a  simple  procedure  to  estimate  the  fractal  dimension  Df  from  the  wavelength 
dependence  of  ju's  [17]. 


Using  equations  (8)  and  (9),  and  the  known  value  of  dmin,  the  variation  of  g  as  a  function  of  dmax 
can  be  plotted  in  a  g  ~  dmax  curve.  Using  the  g  ~  dmax  curve,  the  maximum  correlation  length  dmax 
can  be  obtained  from  the  detennined  anisotropy  factor  g. 
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The  wavelength  dependence  of  jus,  jua,  ju's  and  g  in  the  750-860  mn  spectral  range  is  shown  in  Fig. 

7(a),  Fig.  7(b),  Fig.  7(c),  and  Fig.  7(d),  respectively.  Over  the  entire  wavelength  range  the  value 
of  jus  for  normal  prostate  tissue  is  larger  than  that  for  cancerous  prostate  tissue,  indicating  that 
normal  tissue  scatters  light  more  strongly  than  cancerous  tissue.  The  gradual  decrease  in  the 
value  of  jus  with  wavelength  is  consistent  with  Mie  theory  prediction. 


Wavelength  (nm)  Wavelength  (mn) 


Fig.  7  The  wavelength  dependence  of  (a)  the  scattering  coefficient  (//s)  of  cancerous  (solid  square) 
and  normal  (solid  circle)  prostate  tissues  extracted  from  the  measured  optical  density  (O.D)  data  in 
the  spectral  range  of  750  nm  to  860  nm  using  equation  (5)  in  the  text;  (b)  the  absorption  coefficient 
( jUa )  of  cancerous  (solid  square)  and  normal  (solid  circle)  prostate  tissues.  The  absorption  spectrum 
of  water  with  2  cm  thickness  in  the  spectral  range  of  750nm  -  860nm  is  shown  at  upper  left  location 
as  a  reference;  (c)  the  reduced  scattering  coefficients  (/r')  of  cancerous  (solid  square)  and  normal 

(solid  circle)  prostate  tissues;  and  (d)  the  scattering  anisotropy  factor  g  of  cancerous  (solid  square) 
and  normal  (solid  circle)  prostate  tissues  in  the  spectral  range  of  750nm  -  860nm. 

Fig.  7(b)  shows  that  the  absorption  coefficient  of  normal  prostate  tissue  is  higher  than  that  of  the 
cancerous  tissue  in  the  750  -  860  nm  spectral  range.  Comparison  with  the  absorption  spectrum  of 
water  in  the  inset  shows  that  the  absorption  spectra  of  the  tissues  follow  that  of  water  closely, 
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and  the  behavior  is  more  prominent  for  the  absorption  spectrum  of  the  normal  tissue.  This  result 
indicates  that  water  is  the  dominant  absorptive  tissue  constituent  in  this  spectral  range. 

Reduced  scattering  coefficient  (Fig.  7(c))  of  normal  prostate  tissue  decreases  and  that  of 
cancerous  prostate  tissue  increases  monotonically  in  the  750-860  mn  range,  while  the  behavior  is 
reversed  for  anisotropy  factor  g.  The  overall  magnitude  of  ju's  is  significantly  higher  for 
cancerous  prostate  tissue  than  for  normal  prostate  tissue.  The  observed  lower  values  of  g  for 
cancerous  tissue  than  that  for  normal  tissue  (Fig  7  (d))  follow  readily  from  the  behavior  of  ju's 
because  of  the  simple  algebraic  relation  between  the  two  parameters. 

To  study  how  the  absorption  and  scattering  spectra  reflect  micro-structural  properties  of 
cancerous  and  normal  prostate  tissues,  we  analyzed  the  scattering  data  of  these  two  types  of 
tissues  using  the  fractal  soft  tissue  model  [17].  The  following  parameters  were  used  in  our 
simulation:  the  total  volume  fraction  of  all  scatterers  Tv=  0.2  [16],  the  diameter  of  smallest 
particle  dmin  =  50  mn  [16],  the  index  of  refraction  of  the  cytoplasm  of  prostate  tissue  (host 
medium)  no  =  1-35  [18],  and  the  index  of  refraction  of  prostate  cellular  nuclei  (scattering 
particles)  ti\  =  1.46  [16,  18]. 


Fig.  8  (a)  The  anisotropy  factor  g  of  the  cancerous  and  normal  prostate  tissues  as  a  function  of  cutoff 
diameter  dmax  at  800  nm,  and  (b)  distribution  of  the  number  density  of  scatterers  with  different  diameters, 
evaluated  by  fractal  tissue  model  using  equation  (3).  The  curves  representing  the  behavior  of  cancerous 
and  normal  prostate  tissues  are  displayed  by  solid  and  dash  lines,  respectively.  The  cancerous  prostate 
tissue  (solid  line)  has  lower  weight  of  smaller  particles  and  higher  weight  of  larger  particles  than  the 
normal  prostate  tissue  (dash  line). 

The  values  of  fractal  dimension,  D/,  were  obtained  by  fitting  the  data  shown  in  Fig.  7  (c)  using 
Equation  (10),  and  the  results  show  that  Df= 1.9±0.04  and  4.4±0.06  for  the  cancerous  and  normal 
prostate  tissues,  respectively. 

The  g  values  of  the  cancerous  and  normal  prostate  tissues  at  800  nm  as  a  function  of  cutoff 
diameter  dmax  were  obtained  using  equations  (8)  and  (9),  and  the  results  are  displayed  as  solid  and 
dash  lines,  respectively,  in  Fig.  8  (a).  Using  the  previously  determined  values  of  g  shown  in  Fig. 
7(d),  the  cutoff  diameters  of  the  cancerous  and  normal  prostate  tissues  were  found  to  be  1.42  ± 
0.06  pm  for  the  cancerous  tissue,  and  6.48  ±  0.08  pm  for  the  normal  tissue. 
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The  distribution  of  the  number  density  as  a  function  of  the  diameter  d  of  scattering  particles 
(scatterers)  in  cancerous  and  normal  prostate  tissues  was  then  calculated  using  Equation  (7),  and 
the  results  are  shown  as  Fig.  8(b).  The  salient  feature  of  Fig.  8(b)  is  that  the  nonnal  tissue  has 
much  more  of  smaller  scatterers  than  the  cancerous  tissue  in  the  d  (  ~0.6  pm  range  while  the 
cancerous  tissue  has  much  more  of  larger  scatterers  than  the  normal  tissue  in  the  d  )  -0.6  pm 
range. 

The  fractal  parameters  for  the  normal  and  cancerous  prostate  tissues  are  summarized  in  Table  III. 
Table  III:  Df  and  dmax  for  cancerous  and  normal  prostate  tissues 


Tissue  type 

Df 

dmax  (jAHl) 

nucleolar  diameter 

nuclei  diameter 

Cancer 

1.9F0.04 

1.42  ±0.06 

1.38  -  1.96  pm  [19] 

8.4-10.31  pm  [19] 

Normal 

4.4±0.06 

6.48  ±0.08 

1.10-1.24  pm  [19] 

6.3-6.72  pm  [19] 

As  described  above,  the  fractal  soft  tissue  model  was  used  to  extract  the  values  of  fractal 
dimension  ( Dj )  from  the  wavelength  dependence  of/i'  and  were  found  to  be  1.9  and  4.4  for  the 
cancerous  and  normal  prostate  tissues,  respectively.  The  values  are  markedly  different.  Similar 
remarkable  difference  was  observed  in  the  values  for  cutoff  diameter,  dmax:  1.42  ±  0.06  pm  for 
the  cancerous  tissue,  and  6.48  ±  0.08  pm  for  the  normal  tissue.  Since  the  cutoff  diameter  may  be 
considered  to  be  an  upper  limit  on  the  size  of  scatterers  in  tissue  specimens  this  result  may  be 
quite  significant.  Further  probing  revealed  that  the  observed  value  of  the  cutoff  diameter  of  the 
cancerous  tissue  (1.42  ±  0.06  pm)  falls  within  the  range  of  (1.38  -  1.96  pm)  for  the  nucleolar 
diameter  of  cancerous  prostate  cell  [19],  and  the  cutoff  diameter  for  normal  tissue  (6.48  ±  0.08 
pm)  closely  matches  the  nuclear  diameter  of  nonnal  prostate  cell  (6.3  -  6.72  pm)  [19].  These 
close  agreements  point  to  the  intriguing  implication  that  in  nonnal  prostate  cell  the  largest 
scatterers  are  the  nuclei,  while  in  cancerous  prostate  cells  nucleoli  are  the  largest  scatterers.  - 

(5)  Chemotherapeutic  effects  on  malignant  cells  evaluated  by  native  fluorescence  spectra 

The  activities  reported  in  this  section  are  related  to  Task  1  of  the  proposed  training  program  on 
biological  knowledge  and  technique. 

The  aim  of  this  present  research  is  to  detennine  if  the  native  fluorescence  spectroscopy  is  effective 
enough  to  detect  changes  of  tissue  compositions  (fluorophores)  related  to  the  treatment  of 
chemotherapy  for  malignant  cell  lines  using  Nonnegative  Matrix  Factorization  (NMF)  analysis. 

The  human  breast  cells  were  purchased  from  American  Type  of  Culture  Collection  (ATCC), 
Rockville,  MD.  The  cell  lines  used  in  the  study  were:  malignant  human  breast  cell  lines:  ATCC 
HTB22  (Adenocarinoma  Pleural  Effusion)  and  ATCC  HTB126  (Dutal  Carcinoma).  The  cells  were 
cultured  with  90%  modified  Dulbeco’s  minimal  essential  medium  (Eagle)  (GIBO)  and  10%  of  the 
fetal  bovine  serum  (GIBCO)  in  a  humidified  atmosphere  of  4%  CO2  at  37°C.  Half  of  the  cells  were 
treated  with  -10  6  M  retinoic  acid  (RA)  and  the  remainder  served  as  controls  while  the  cell  lines  are 
cultured.  The  measured  cell  concentration  was  estimated  as  -105  -  1 10  x  106 cells/ml  [5]. 

The  fluorescence  emission  spectra  of  cells  were  measured.  The  excitation  wavelengths  of  300nm 
and  340  nm  were  selected  based  on  our  previous  spectral  study  of  the  main  fluorophores  in  breast 
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cells  [20].  Combined  with  NMF,  the  native  fluorescence  measurements  were  used  to  distinguish 
the  RA-treated  breast  cancer  cells  from  the  untreated  cell.  The  main  advantage  of  NMF  is  that  the 
NMF-recovered  spectral  data  and  concentrations  of  constituents  are  positive  values.  This  makes 
the  calculated  results  represent  more  physically  and/or  biologically  meaning  because  the  spectra 
and  contents  of  the  fluorophores  should  have  non-negative  values. 

Fourteen  pairs  of  the  RA-treated  and  untreated  breast  cancer  cell  samples  were  investigated 
using  fluorescence  spectroscopy  with  excitation  wavelength  of  300  nm.  The  average 
fluorescence  spectral  profiles  with  standard  deviation  error  bars  at  key  wavelengths  for  the  RA- 
untreated  and  treated  cancerous  breast  cells  are  shown  in  Figs.  9(a)  and  9(b),  respectively.  Each 
spectral  profile  was  normalized  to  unit  value  of  1  (i.e.,  the  sum  of  squares  of  the  elements  in  each 
emission  spectra  data  was  set  as  1)  before  taking  average  and  calculation. 


Fig.  9  Average  fluorescence  spectra  of  (a)  untreated  and  (b)  Retinoic  Acid  (RA)-treated 
cancerous  breast  cells  obtained  with  the  excitation  of  300  nm. 


The  main  emission  peaks  of  both  untreated  and  treated  cancerous  breast  cells  were  found  at 
337nm.  The  major  difference  of  the  fluorescence  profiles  between  the  untreated  and  treated  cells 
is  that  there  exists  a  much  lower  local  peak  at  ~  460nm  for  the  untreated  cells  while  no  such  a 
peak  for  the  treated  cells.  A  tiny  local  peak  at  529  nm  may  be  observed  in  the  emission  of  the 
untreated  cells,  which  can  be  more  clearly  seen  by  enlarging  the  fluorescence  profile  at  the 
spectral  range  of  700  nm  to  900  nm  as  shown  in  insert  of  Fig.  9  (a). 

To  investigate  the  relative  content  changes  of  Principal  Components  (PCs)  and  compare  with 
their  corresponding  fluorophores,  NMF  was  used  to  extract  spectra  of  PCs  (fluorophores)  of  the 
1st  component  (tryptophan),  the  2nd  (NADH)  and  the  3ld  (flavin)  for  the  cells,  which  are  shown  as 
dash,  dot  and  dash-dot  lines  in  Fig.  10(a),  respectively.  The  measured  spectra  of  individual 
tryptophan,  NADH  and  flavin  in  solutions  are  plotted  as  thicker,  thick  and  thin  solid  lines, 
respectively,  as  references.  The  two  groups  of  spectra  (the  extracted  and  the  measured)  for  each 
fluorophore  show  reasonable  agreement,  which  demonstrates  that  NMF  model  accounts  for  the 
major  spectroscopic  feature  observed,  and  indicates  that  the  measurements  are  reasonable. 

To  investigate  the  relative  content  changes  of  PCs,  the  contents  of  tryptophan,  NADH  and  flavin 
in  two  types  of  cells  were  extracted  from  the  measured  total  fluorescence  spectra  using  the  NMF 
analysis.  Fig.  10  (b)  shows  the  relative  content  of  the  1st  PC  (tryptophan)  vs.  the  2nd  PC  (NADH); 
10  (c)  displays  the  3ld  PC  (flavin)  vs.  the  1st  PC  (tryptophan)  and  10  (d)  exhibits  the  3ld  PC 
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(flavin)  vs.  the  2nd  PC  (NADH)  of  the  untreated  (solid  circles)  and  treated  (solid  squares) 
cancerous  breast  cells.  The  most  salient  feature  of  Fig.  10  (b)  is  that  all  data  points  for  the 
untreated  cells  are  located  on  the  upper  side  over  the  data  points  for  the  treated  cells,  indicating 
that  the  relative  contribution  of  tryptophan  is  higher  in  the  untreated  cancerous  cells  in 
comparison  with  the  treated.  Fig.  10  (c)  shows  again  that  the  relative  fluorescence  contribution  of 
tryptophan  is  higher  in  the  untreated  cancerous  cells  in  comparison  with  the  treated.  Fig.  10  (d) 
shows  that  most  data  points  for  the  untreated  cells  are  located  in  the  up-right  side  in  comparison 
with  the  data  points  for  the  treated  cells,  indicating  that  the  relative  contents  of  NADFI  and  flavin 
are  higher  in  the  untreated  cancerous  cells  in  comparison  with  the  treated. 
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Fig.  10  (a)  Comparison  of  the  extracted  spectra  of  the  three  PCs:  the  1st  PC  -  tryptophan  (dash  line), 
2nd  PC  -  NADFI  (dot  line)  and  3ld  PC  -  flavin  (dash-dot  line)  for  the  untreated  and  treated  breast 
cancerous  cells,  and  the  measured  spectra  of  individual  tryptophan  (thicker  line),  NADFI  (thick  line) 
and  flavin  (thin  line)  in  solution;  (b)  Fractional  content  of  the  1st  PC  -  tryptophan  vs.  that  of  the  2nd  PC 
-  NADFI;  (c)  Fractional  content  of  the  3rd  PC  -  flavin  vs.  that  of  the  1st  PC  -  tryptophan;  (d)  Fractional 
content  of  the  3rd  PC  -  flavin  vs.  that  of  the  2nd  PC  -  NADFI  by  analyzing  emission  spectra  obtained 
with  the  excitation  of  -  300  nm  using  NMF  method.  The  data  for  the  untreated  and  treated  cells  are 
displayed  with  the  solid  circles  and  solid  squares  in  10  (b)  -  10  (d),  respectively. 


The  following  biomedical  studies  may  help  to  understand  these  changes.  Brown  et  al  studied 
chemotherapy  effects  for  breast  cancer  [21].  After  chemotherapy,  the  decreased  plasma  amounts 
of  tryptophan  in  plasma  were  observed  [21]  since  increased  tryptophan  levels  are  somehow 
related  to  increased  tumor  cell  proliferation  [21].  NADH  is  one  of  the  most  important  coenzyme 
for  the  production  of  Adenosine-Tri -Phosphate  (ATP)  [22].  The  results  of  decrease  of  NADH  in 
the  RA-treated  cancerous  cells  may  reflect  a  reduction  in  the  production  of  ATP  in  these  treated 
cells.  Since  the  decrease  of  tryptophan  and  NADH  somehow  reflect  the  less  activity  of  cells  [21, 
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22],  it  can  be  hypothesized  that  the  RA-treatment  results  in  cancerous  cell  death  or  reduction  of 
the  activity  of  cancer  cells. 


Emission  spectra  of  nine  pairs  of  the  RA-treated  and  untreated  breast  cancer  cell  samples  were 
measured  with  an  alternate  excitation  of  340  nm.  The  average  fluorescence  spectral  profiles  of 
the  untreated  (solid  line)  and  treated  (dash  line)  cells  are  displayed  in  Figs.  11(a),  which  shows 
that  the  emission  profiles  of  the  untreated  and  treated  cells  were  mainly  contributed  from  NADH 
corresponding  to  a  peak  at  ~  459nm,  and  flavin  corresponding  to  a  peak  of  ~  530nm.  The  major 
difference  of  the  profiles  between  the  untreated  and  treated  cells  is  that  the  shoulder  peak  at  ~ 
530  nm  for  the  treated  cells  is  little  stronger  than  that  for  the  untreated  cells.  Fig.  1 1  (b)  displays 
the  relative  content  of  the  flavin  vs.  NADH  of  the  untreated  (solid  circle)  and  treated  (solid 
square)  cancerous  breast  cells  analyzed  using  the  NMF  analysis.  The  obvious  feature  of  Fig. 
11(b)  is  that  the  relative  contribution  of  NADH  is  higher  in  the  untreated  cancerous  cells  in 
comparison  with  the  treated,  which  also  reflects  the  decrease  of  the  relative  content  of  NADH  in 
the  treated  cancerous  breast  cells  in  comparison  with  the  untreated  cells. 


Wavelength(nm) 


concentration  of  NADH 


Fig.  11  (a)  Average  fluorescence  spectra  of  the  untreated  (solid)  and  RA-treated  (dash) 
cancerous  breast  cells  obtained  with  the  excitation  wavelength  of  340  nm.(b)  Relative  content 
of  the  NADH  vs.  flavin,  obtained  from  the  measured  fluorescence  spectroscopy  of  the  breast 
cells  with  340  nm  excitation  using  the  NMF  analysis. 


This  study  shows  detection  of  the  chemotherapeutic  effects  on  malignant  breast  cells  by 
measuring  the  nature  fluorescence  spectra  of  cells  and  extracting  the  fluorescence  spectral 
features  of  principal  components  using  an  NMF  method.  The  presence  and  contents  of  the 
detectable  fluorophores  in  breast  malignant  cells  such  as  tryptophan,  NADH  and  flavin  were 
analyzed  as  the  Principal  Components  for  the  NMF.  We  demonstrated  that  fluorescence 
intensities  of  tryptophan,  NADH  and  flavin  decreased  in  the  RA-treated  cells  in  comparison  with 
those  in  the  untreated  cells.  This  work  shows  the  change  of  relative  contents  of  tryptophan, 
NADH  and  flavin  studied  using  native  fluorescence  spectroscopy  with  the  NMF  analysis  may 
present  potential  criteria  for  evaluation  of  chemotherapeutic  effects  on  malignant  breast  cells. 

(6)  Specific  personnel  supported  by  this  grant: 
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postdoctoral  fellow 


Dr.  Yang  Pu,  PI, 

Key  Accomplishments 

•  The  training  program  on  biological  knowledge  and  technique  has  been  accomplished 
through  laboratory  rotations  and  attending  group  meetings  and  seminars  at  MSKCC, 
including  small  animal  training  sessions  ( RARC  Orientation,  RARC  Basic  Mouse,  RARC 
Rodent  Survival  Surgery),  training  on  cell  culture,  flow  cytometry,  fluorescence 
microscopic  imaging  techniques,  bioluminescence  assay  application,  cells  viability 
monitoring,  and  other  key  techniques  used  in  cancer  research. 

•  The  training  program  on  time-resolved  polarization  and  steady-state  spectroscopic 
techniques  has  been  accomplished  through  hands-on  training  at  the  mentor’s  laboratory 
in  the  Institute  for  the  Ultrafast  Spectroscopy  and  Lasers  (IUSL)  of  CUNY.  Time- 
resolved  polarization  and  steady-state  spectroscopic  measurements  of  dye-polymer 
conjugate,  target  contrast  agents  in  solution  and  in  scattering  medium  have  been 
perfonned. 

•  The  time -resolve  fluorescence  depolarization  caused  by  the  Brownian  rotation  of  the  dipoles 
was  investigated  experimentally  and  theoretically.  For  simple  and  stable  molecules,  the 
effect  of  Brownian  rotation  is  the  main  cause  of  the  depolarization.  We  show  this  effect  over 
a  limited  range  of  the  molecular  volume  of  Fluorescein  dye  attached  with  different  polymer 
chains. 

•  A  receptor  target  contrast  agent,  Cytate,  was  demonstrated  as  a  prostate  cancer  receptor- 
targeted  contrast  agent  using  NIR  scanning  imaging  technique.  The  results  show  the 
preferential  uptake  of  Cytate  by  cancerous  prostate  tissue  in  comparison  with  normal 
prostate  tissue.  The  results  also  indicate  that  OPTICA  can  be  used  to  enhance  the  contrast  of 
emission  intensities  between  cancerous  and  normal  tissue  areas.  These  techniques  (NIR 
imaging  and  OPTICA)  can  also  be  applied  to  obtain  the  three-dimensional  locations  of 
targets  (Cytate-stained  cancerous  tissues)  in  large  host  normal  prostate  tissue. 

•  A  systematic  study  was  performed  on  cancerous  and  normal  prostate  tissues  to  obtain  the 
scattering  coefficient  (jus),  the  absorption  coefficient  (jua),  the  reduced  scattering  coefficient 
( n's ),  and  the  anisotropy  factor  (g).  The  values  of  fractal  dimension  (D/)  of  cancerous  and 
normal  prostate  tissues  were  obtained  by  fitting  the  wavelength  dependence  of/r' .  The 

cutoff  diameter  dmax  as  a  function  of  g  were  investigated  using  the  fractal  soft  tissue  model 
and  Mie  theory.  Results  show  that  dmax  of  the  nonnal  tissue  is  larger  than  that  of  the 
cancerous  tissue.  The  cutoff  diameter  dmax  is  observed  to  agree  with  the  nuclear  size  for  the 
normal  tissue  and  the  nucleolar  size  for  the  cancerous  tissue. 

•  The  chemotherapeutic  effects  on  malignant  cells  were  evaluated  by  measuring  the  native 
fluorescence  spectra  of  cells  and  extracting  the  fluorescence  spectral  features  of  principal 
components  using  an  NMF  method.  The  presence  and  contents  of  the  detectable 
fluorophores  in  breast  malignant  cells  such  as  tryptophan,  NADH  and  flavin  were  analyzed 
as  the  Principal  Components  for  the  NMF.  We  demonstrated  that  fluorescence  intensities  of 
tryptophan,  NADH  and  flavin  decreased  in  the  RA-treated  cells  in  comparison  with  those  in 
the  untreated  cells.  This  work  shows  the  change  of  relative  contents  of  tryptophan,  NADH 
and  flavin  studied  using  native  fluorescence  spectroscopy  with  the  NMF  analysis  may 
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present  potential  criteria  for  evaluation  of  chemotherapeutic  effects  on  malignant  breast 
cells. 

Reportable  Outcomes 

•  Three  papers  in  peer-reviewed  journals  and  five  papers  in  conference  proceeding  resulted 
from  the  research  supported  by  and  related  to  this  grant  have  been  published  in  the  first  year. 

•  Eight  conference  and  symposium  presentations  resulted  from  the  research  supported  by  this 
grant  have  been  made. 

The  following  lists  the  publications  and  presentations: 

(1)  Paper  Publications  in  peer-reviewed  journals: 

1.  Yang  Pu,  Wubao  Wang  Mohammad  AL-Rubaiee,  Swapan  Kumar  Gayen,  and  Min  Xu, 
“Determination  of  optical  coefficients  and  fractal  dimensional  parameters  of  cancerous  and 
nonnal  prostate  tissues”,  accepted  by  Appl.  Spectrosc.,  to  be  published  on  July  issue  66  (7), 
(2012). 

2.  E.  Jeyasingh,  Yang  Pu,  W.  B.  Wang,  C.  H.  Liu,  and  R.  R.  Alfano,  “Stokes  Shift 
Spectroscopy  Pilot  Study  for  Cancerous  and  Normal  Prostate  Tissues”,  accepted  by  Appl. 
Opt.,  to  be  published  on  June  issue  5J_  (2012). 

3.  Y.  Pu,  W.  B.  Wang,  M.  Xu,  G.  C.  Tang,  Y.  Budansky,  M.  Sharanov,  S.  Achilefu,  J.  A. 
Eastham  and  R.  R.  Alfano,  “Near  Infrared  Photonic  Finger  Imager  for  Prostate  Cancer 
Screening”,  Technol.  Cancer  Res.  Treat.,  ( TCRT ),  f0,  selected  as  a  cover  page  paper,  507- 
517(2011). 

(2)  Paper  Publications  in  conference  proceedings: 

1.  Yang  Pu,  Guichen  Tang,  Wubao  Wang,  H.  E.  Savage,  S.  P.  Schantz,  and  Robert  R.  Alfano, 
“Chemotherapeutic  Effects  on  Breast  Malignant  Cells  Evaluated  by  Native  Fluorescence 
Spectroscopy”,  in  Biomedical  Optics  (BIOMED)/  Digital  Holography  and  Three- 
Dimensional  Imaging  (DH),  JM3A45-1-3  (2012). 

2.  Y.  Pu,  W.  B.  Wang,  R.  B.  Dorshow,  R.  R.  Alfano,  “Picosecond  polarization  spectroscopy  of 
fluorescein  attached  to  different  molecular  volume  polymer  influenced  by  rotational  motion” 
in  Organic  Photonic  Materials  and  Devices  XIV,  edited  by  Christopher  Tabor,  Francois 
Kajzar,  Toshikuni  Kaino,  Yasuhiro  Koike,  Proceedings  of  SPIE  Vol.  8258  (SPIE, 
Bellingham,  WA  2012)  825818. 

3.  Y.  Pu,  W.  B.  Wang,  G.  C.  Tang,  Y.  Budansky,  M.  Sharonov,  M.  Xu,  S.  Achilefu,  J.A. 
Eastham,  R.  R.  Alfano,  “Screening  prostate  cancer  using  a  portable  near  infrared  scanning 
imaging  unit  with  an  optical  fiber-based  rectal  probe”  in  Optical  Biopsy  X,  edited  by  Robert 
R.  Alfano,  Stavros  G.  Demos,  Proceedings  of  SPIE  Vol.  8220  (SPIE,  Bellingham,  WA  2012) 
822002. 

4.  Y.  Pu,  W.  B.  Wang,  M.  Alrubaiee,  S.  K.  Gayen,  M.  Xu,  “Systematic  studies  of  fractal 
dimension  parameters,  absorption  and  scattering  coefficients  for  cancerous  and  nonnal 
prostate  tissues”  in  Photonic  Therapeutics  and  Diagnostics  VIII,  edited  by  Hyun  Wook 
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Kang,  Bodo  E.  Knudsen,  Proceedings  of  SPIE  Vol.  8207  (SPIE,  Bellingham,  WA  2012) 
8207 1H. 

5.  C.-H.  Liu,  X.  H.  Ni,  Y.  Pu,  Y.  L.  Yang,  F.  Zhou,  R.  Zuzolo,  W.  B.  Wang,  V.  Masilamani,  A. 
Rizwan,  R.  R.  Alfano,  “Optical  spectroscopic  characteristics  of  lactate  and  mitochondrion  as 
new  biomarkers  in  cancer  diagnosis:  understanding  Warburg  effect”  in  Optical  Biopsy  X, 
edited  by  Robert  R.  Alfano,  Stavros  G.  Demos,  Proceedings  of  SPIE  Vol.  8220  (SPIE, 
Bellingham,  WA  2012)  82200Y. 

(3)  Presentations  and  Abstracts: 

1 .  R.  V.  Simoes,  E.  Ackerstaff,  N.  Kruchevsky,  Y.  Pu,  G.  Sukenich,  and  J.  Koutcher,  “Dynamic 
metabolic  signature  of  metastatic  and  non-metastatic  breast  cancer  cells”  ISMRM  2012 
Annual  Meeting,  5638,  May  5-11,  (2012). 

2.  Yang  Pu,  Guichen  Tang,  Wubao  Wang,  H.  E.  Savage,  S.  P.  Schantz,  and  Robert  R.  Alfano, 
“Chemotherapeutic  Effects  on  Breast  Malignant  Cells  Evaluated  by  Native  Fluorescence 
Spectroscopy”,  presented  in  Biomedical  Optics  ( BIOMED )/  Digital  Holography  and  Three- 
Dimensional  Imaging  (DH)  OSA  Optics  &  Photonics  Congress,  held  in  Miami,  Florida,  on 
April  29  -  May  02,  2012. 

3.  Y.  Pu,  W.  B.  Wang,  R.  B.  Dorshow,  R.  R.  Alfano,  “Picosecond  polarization  spectroscopy  of 
fluorescein  attached  to  different  molecular  volume  polymer  influenced  by  rotational  motion” 
presented  at  the  2012  OPTO,  in  SPIE  Photonics  West  2012,  Conference  8258:  Organic 
Photonic  Materials  and  Devices  XIV,  held  in  San  Francisco,  CA,  on  January  25,  2012.  In  the 
Technical  Program,  PW12O-OE104-3. 

4.  Y.  Pu,  W.  B.  Wang,  M.  Alrubaiee,  S.  K.  Gayen,  M.  Xu,  “Systematic  studies  of  fractal 
dimension  parameters,  absorption  and  scattering  coefficients  for  cancerous  and  nonnal 
prostate  tissues”  presented  at  the  2012  International  Biomedical  Optics  Symposium  (BIOS), 
SPIE  Photonics  West  2012,  Conference  BO101:  Photonic  Therapeutics  and  Diagnostics  VIII, 
held  in  San  Francisco,  CA,  on  January  21,  2012.  In  the  Technical  Program,  PW12B-BO101- 
53. 

5.  Y.  Pu,  W.  B.  Wang,  G.  C.  Tang,  Y.  Budansky,  M.  Sharonov,  M.  Xu,  S.  Achilefu,  J.  A. 
Eastham,  R.  R.  Alfano,  “Screening  prostate  cancer  using  a  portable  near  infrared  scanning 
imaging  unit  with  an  optical  fiber-based  rectal  probe”  presented  at  the  2012  International 
Biomedical  Optics  Symposium  (BIOS),  SPIE  Photonics  West  2012,  Conference  BO207: 
Optical  Biopsy  X,  held  in  San  Francisco,  CA,  on  January  24-25,  2012.  In  the  Technical 
Program,  PW 1 2B-BO207- 1 . 

6.  Yang  Pu,  “Prostate  cancer  detection  using  receptor-targeted  contrast  agents  and  a  near 
infrared  scanning  imaging  unit  with  fiber-based  rectal  probe”,  presented  at  “December 
Career  Day”,  the  NYAS/PepsiCo  Day  in  NY  December  12,  2011,  In  the  Technical  Program, 
112. 

7.  Y.  Pu,  W.  B.  Wang,  S.  Achilefu,  M.  Xu,  and  R.  R.  Alfano,  “Detecting  prostate  cancer  in 
three  dimensional  position  using  a  receptor-targeted  contrast  agent  and  near  infrared 
scanning  imaging”,  presented  as  the  NIH  Workshop  2011  7th  Inter-Institute  Workshop  on 
Optical  Diagnostic  and  Biophotonic  Methods  from  Bench  to  Bedside,  held  in  Natcher 
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Conference  Center,  National  Institutes  of  Health  Bethesda,  Maryland,  ob  September,  15-16, 
201 1,  In  the  Technical  Program,  NIH  100-68 

8.  Yang  Pu,  W.  B.  Wang,  and  R.  R.  Alfano  “Detecting  Cancerous  Prostate  Tissues  in  3D 
position  after  Enhanced  by  a  Receptor-targeted  Contrast  Agent  using  Near  Infrared  Scanning 
Imaging  and  Independent  Component  Analysis”  in  The  City  College  of  New  York  2011 
Structure  Biology  Faculty  Retreat,  held  in  Concourse  Level  of  The  Graduate  Center  of 
CUNY,  on  June  6,2011. 

Conclusions: 

1.  Successfully  accomplished  Task  1  of  planned  training  objectives  through  laboratory  rotations 
and  attending  group  meetings  and  seminars  at  MSKCC.  The  biological  and  biomedical 
training  include  small  animal  training  sessions  (RARC  Orientation,  RARC  Basic  Mouse, 
RARC  Rodent  Survival  Surgery),  training  on  cell  culture,  flow  cytometry,  fluorescence 
microscopic  imaging  techniques,  bioluminescence  assay  application,  cells  viability 
monitoring,  and  other  key  techniques  used  in  cancer  research. 

2.  Successfully  accomplished  Task  3  of  planned  training  objectives  through  hands-on 
laboratory  training  at  mentor’s  Lab  in  the  Institute  for  the  Ultrafast  Spectroscopy  and  Lasers 
(IUSL)  of  CCNY.  The  biomedical  optics  training  includes  time-resolved  polarization  and 
steady-state  spectroscopic  measurements  of  dye-polymer  conjugate,  target  contrast  agents  in 
solution  and  in  scattering  medium,  the  theory  of  time-resolve  fluorescence  depolarization 
caused  by  the  Brownian  rotation  of  the  dipoles. 

3.  Started  Task  2  of  planned  training  objectives  through  hands-on  laboratory  training  at 
Washington  University  School  of  Medicine  at  St  Louis  (WUSM).  The  biochemistry 
knowledge  and  techniques  training  include  cellular  and  molecular  assays,  cell  labeling,  target 
contrast  agent  synthesizing.  (Because  conflicting  schedule,  Task  3  was  arranged  to  be 
performed  before  Task  2.  However,  task  2  was  partially  achieved  by  remote  mentoring  and 
Lab-training  was  started  on  April  1 1th,  2012) 

4.  Following  the  training  program  of  tasks  1  to  3,  the  following  research  programs  have  been 
carried  out:  (1)  the  time -resolve  fluorescence  depolarization  caused  by  the  Brownian  rotation 
of  the  dipoles  was  investigated  experimentally  and  theoretically;  (2)  a  receptor  target  contrast 
agent,  Cytate,  was  demonstrated  as  a  prostate  cancer  receptor-targeted  contrast  agent  using 
NIR  scanning  imaging  technique,  and  NIR  imaging  and  an  algorithm  named  OPTICA  were 
applied  to  obtain  the  three-dimensional  locations  of  Cytate-stained  cancerous  tissues  in  large 
host  nonnal  prostate  tissue;  (3)  a  systematic  study  was  performed  on  cancerous  and  normal 
prostate  tissues  to  obtain  /us,  jua,  ju's ,  and  g  of  cancerous  and  nonnal  prostate  tissues.  The 

values  of  fractal  dimension  (Dj)  and  the  cutoff  diameter  dmax  of  cancerous  and  nonnal 
prostate  tissues  were  investigated  using  the  fractal  soft  tissue  model  and  Mie  theory;  and  (4) 
the  chemotherapeutic  effects  on  malignant  cells  were  evaluated  by  measuring  the  native 
fluorescence  spectra  of  cells  and  extracting  the  fluorescence  spectral  features  of  principal 
components  using  an  NMF  method.  The  fluorescence  intensities  of  tryptophan,  NADH  and 
flavin  were  observed  to  be  decreased  in  the  RA-treated  cells  in  comparison  with  those  in  the 
untreated  cells. 
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Plan  for  the  Coming  Year 

The  research  tasks  for  the  second  year  (04/15/2012  -  04/14/2013)  are  described  by  Tasks  2,  4,  5 
of  the  Statement  of  Work.  In  the  coming  year,  we  will: 

(1)  Continue  the  Pi’s  training  on  basic  techniques  for  cellular  assay  synthesizing  in  the  lab  of 
Professor  Alchilefu;  study  optical  properties  of  different  contrast  agents  to  choose  derisible 
spectra  for  synthesis;  and  synthesize  2-3  kinds  of  objective  receptor  target  contrast  agents,  which 
will  be  used  to  target  corresponding  over-expressed  receptors  in  prostate  cancer  cells  under 
direction  of  Professor  Alchilefu. 

(2)  Submicron  resolution  microscopic  study  of  targeting  contrast  agents  conjugated  with  prostate 
cancer  cells  using  Confocal  microscopy;  training  of  operations  of  AlphaA/SOM  under  direction  of 
Professor  Alfano;  exploiting  the  specificity  and  affinity  of  synthesized  contrast  agents  with 
prostate  cancer  cells  monitoring  using  Confocal  microscope;  Submicron  resolution  microscopic 
study  of  targeting  contrast  agents  conjugated  with  prostate  cancer  cells;  control  experiments  with 
nonnal  prostate  cells;  develop  new  method  to  image  submicron  resolution  of  receptor  expression 
on  cellular  and  molecular  level. 

(3)  Time-resolved  polarization  and  steady-state  spectroscopic  study  of  targeting  contrast  agents 
to  prostate  cancer  cells  to  investigate  the  affinity  between  target  contrast  agents  and  prostate 
cancer  cells;  research  on  different  kinds  of  over-expressed  receptors  of  prostate  cancer  cell  lines; 
time-resolved  polarization  and  steady-state  spectroscopic  measurements  of  synthesized  targeting 
contrast  agents;  time-resolved  polarization  and  steady-state  spectroscopic  measurements  of 
prostate  cancer  cells;  time-resolved  polarization  and  steady-state  spectroscopic  measurements  of 
t  contrast  agents  conjugated  with  prostate  cancer  cells;  study  conjugation  of  targeting  contrast 
agents  with  prostate  cancer  cells  based  on  the  time-resolved  polarization  spectroscopic 
measurements  and  empirical  theory  of  rotational  dynamics  and  fluorescence  polarization 
anisotropies  in  scattering  medium;  and  Control  experiments  with  nonnal  prostate  cells. 
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A  abstract  of  presentation  entitled  “Dynamic  metabolic  signature  of  metastatic  and  nonmetastatic 
breast  cancer  cells”  by  R.  V.  Simoes,  E.  Ackerstaff,  N.  Kruchevsky,  Y.  Pu,  G.  Sukenich,  and  J. 
Koutcher,  presented  in  ISMRM  2012  Annual  Meeting,  5638,  May  5-11,  (2012). 
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DYNAMIC  METABOLIC  SIGNATURES  OF  METASTATIC  AND  NON-METASTATIC  BREAST  CANCER  CELLS 

Rui  Vasco  Simoes1,  Ellen  Ackerstaff1,  Natalia  Kruchevsky1,  Yang  Pu1  George  Sukenick1,  and  Jason  A  Koutcher1 
1 Memorial  Sloan-Kettering  Cancer  Center,  New  York,  NY,  United  States,  :The  City  College  of  New  York,  New  York,  NY,  United  States 


INTRODUCTION:  Distinguishing  metastatic  from  non-metastatic  clinical  breast  cancer,  non-in vasively,  would  significantly  improve 
evaluation  of  patient  prognosis.  Recent  studies  with  two  isogenic  murine  breast  cancer  lines,  metastatic  4T 1  and  non-metastatic  67NR, 
revealed  differences  at  the  level  of  (a)  LDH-A  expression,  in  response  to  changing  oxygen  tensions  in  vitro,  and  (b)  lactate 
accumulation,  during  tumor  growth  in  vivo  (1).  Understanding  how  each  cell  line  metabolically  adapts  to  dynamic  microenvironmental 
changes,  such  as  oxygen  availability  vs.  hypoxia,  and  major  nutrient  supply  vs.  deprivation  (glucose  and  glutamine  (2,3)),  could  provide 
unique  non-invasive  biomarkers  of  metastatic  and  non-metastatic  breast  cancer  phenotypes,  useful  in  developing  new  therapies. 

PURPOSE:  To  study  4T1  and  67NR  breast  cancer  lines  in  an  NMR-compatible  cell  perfusion  system  and  investigate  differences  in 
their  dynamic  metabolic  responses  to  (a)  oxygen  availability  vs.  hypoxia,  and  (b)  glutamine  (Gin)  supply  vs.  Gin  deprivation. 


METHODS:  4T1  and  67NR  cells  were  initially  studied  in  vitro  to  determine  the  response  of  cell  growth  to  modified  media  composition 
(DME,  supplemented  with  1%  P/S  and  10%  FBS:  DMECOmpi):  1)  DMEcompi,  contains  25  mM  glucose  (Glc)  and  6  mM  Gin;  2)  DME  with  2 
mM  Gin;  3)  DME  without  Gin;  4)  DME  without  Glc.  Cells  were  incubated  for  48  h  in  each  medium,  harvested,  and  counted  (ViaCount 
Assay,  Guava  Technologies).  Three  independent  experiments,  with  duplicate  samples  per  experiment,  were  performed  for  each  cell 
line.  Slight  modifications  were  introduced  in  our  cell  perfusion  system  described  previously  (4)  to  accommodate  its  use  in  a  500  MHz 
vertical-bore  Bruker  magnet.  For  each  study,  3.0  xIO6  cells  were  seeded  on  microcarriers  (Plastic  Plus,  SoloHill)  and  cultured  for  3 
days.  Cells  were  then  loaded  into  the  perfusion  system:  a  custom-made  10  mm  screw-cap  tube  where  culture  medium  and  gas  were 
circulated  and  the  temperature  controlled  (37  °C)  by  a  heating  water  jacket.  The  MR  experiments  consisted  in  repeatedly  acquiring  13C, 
1H,  and  31P  spectra  (standard  Bruker  BBO  probe)  while  challenging  the  perfused  cells  with  different  metabolic  perturbations.  These 
were  carried  out  in  six  steps  (I  -  VI,  5  h  each):  I,  DMEcompi  and  carbogen  gas  (95%  O2  /  5%  CO2);  II,  100%  13C(1)-Glc  labeled  medium, 
without  Gin,  and  carbogen  gas;  III,  same  medium  as  in  II  but  during  hypoxia  (95%  N2  /  5%  CO2);  IV,  still  hypoxia  but  replacing  the 
medium  with  100%  13C(1)-Glc  and  2  mM  Gin;  V,  same  medium  as  in  IV  but  shifting  the  gas  back  to>  carbogen;  VI,  same  as  in  I.  Each 
experiment  was  repeated  3  times  for  each  cell  line.  Spectral  analysis  of  peak  areas  we  carried  out  with  AMARES  (jMRUI  v4.0). 
Statistical  analyses  were  carried  out  with  the  paired  t-Test  (SPSS  15.0)  to  compare  4T1  vs.  67NR  data  (significance  level:  p<0.05). 


RESULTS:  As  shown  in  Fig.  1,  the  in  vitro  cell  growth  is  dependent  on  the  availability  of  both  Glc  and  Gin,  in  4T1  and  67NR  lines:  cell 
growth  drops  by  85-95%  in  the  absence  of  either  Gin  or  Glc,  and  remains  at  35-45%  when  Gin  alone  is  decreased  by  2/3;  no  significant 
differences  were  observed  in  growth  between  the  two  cell  lines  in  response  to  each  nutrient  deprivation.  As  for  dynamic  cell  perfusion 
studies,  31P-monitored  ATP  levels  were  kept  within  ±27%  of  their  initial  concentration  during  the  experiments  (l-VI),  in  both  cell  lines, 
although  showing  significant  differences  during  stages  IV-VI  (data  not  shown),  ^-visible  total  choline  compounds  (3.21  ppm)  and 
mobile  lipids  (1.29  ppm,  mixed  contribution  with  lactate)  showed  consistent  stress-response  changes  at  each  stage,  similar  in  both  cell 
lines  (data  not  shown).  13C-labeling  studies  revealed  that67MR  cells  are  significantly  more  glycolytic  than  4T1  cells  in  the  experimental 
conditions  tested:  in  67NR  cells,  the  rate  of  glucose-derived  lactate  (Lac)  production  is  2-fold  higher  than  in  4T1  cells.  The  glycolytic  Lac 
synthesis  during  hypoxia  is  Gin-dependent  in  4T1  cells  (2-fold  higher  rate  in  presence  of  glutamine:  vs.  deprivation),  but  not  in  67NR 
cells.  Supply  of  Gin  vs.  deprivation,  generates  significantly  different  dynamic  responses  between  4T1  and  67NR  cells,  as  far  as  turnover 
of  glycolytic  intermediates,  e.g.  glyceraldehyde-3-phosphate,  alanine  synthesis,  and  TCA  cycle  activity.  In  aerobic  conditions  and 
presence  of  Gin,  4T1  cells  incorporate  more  Glc  into  the  TCA  cycle  (3-fold  higher  rate  of  glutamate  synthesis)  than  67NR  cells, 
whereas  the  latter  show  higher  glycolytic  activity,  i.e.  2.4-fold  higher  rate  of  Lac  synthesis,  and  also  25-fold  faster  synthesis  of  alanine. 


Figure  1  -  Cell  growth  (A)  and  cell  perfusion 
studies  (B-E)  of  4T1  (red)  and  67NR  (green) 
lines.  A,  Stress  conditions,  from  left  to  right: 
DMEcompi;  Gin  deprivation;  2mM  Gin;  Glc 
deprivation.  B-E,  13C-NMR  monitoring:  I  (0-5h), 
Gin  deprivation;  II  (5-1  Oh),  hypoxia;  III  (1 0-1 5h), 
2mM  Gin;  IV  (15-20h),  oxygen  supply.  Chemical 
Shifts  (ppm):  glyceraldehyde-3-P  C3,  65.98; 
lactate  C3,  21.18;  alanine  C3,  17.30;  glutamate 
C4,  34.59.  Time-points  normalized  to  the  initial 
value  of  each  curve:  *,  p<0.05.  Slopes  (synthesis/ 
consumption  rates):  #,  r2>0.8;  ##,  r>0.9. 
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CONCLUSIONS:  Although  no  significant  differences  in  cell  growth  were  observed  in  vitro  between  4T1  and  67NR  cells,  when  exposed 
to  defined  metabolite  stress  conditions  each  cell  line  adopts  different  metabolic  strategies  to  cope  with  said  changes,  depending  on  the 
level  of  oxygen  availability.  Mostly,  non-metastatic  67NR  cells  revealed  higher  glycolytic  activity  than  metastatic  4T1  cells.  Unlike  in 
67NR  cells,  glycolysis  is  Gin-dependent  in  the  more  aggressive  cell  line  (4T1)  during  hypoxia,  dropping  when  oxygen  becomes 
available  to  favor  TCA  cycle  activity.  Ongoing  13C(3)-Gln  studies  should  provide  further  insight  into  these  metabolic  differences. 
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Optical  extinction  and  diffuse  reflection  spectra  of  cancerous  and  normal 
prostate  tissues  in  the  750  to  860  nm  spectral  range  were  measured. 
Optical  extinction  measurements  using  thin  ex  vivo  prostate  tissue  samples 
were  used  to  determine  the  scattering  coefficient  (ps),  while  diffuse 
reflection  measurements  using  thick  prostate  tissue  samples  were  used  to 
extract  the  absorption  coefficient  (pa)  and  the  reduced  scattering 
coefficient  (p').  The  anisotropy  factor  (g)  was  obtained  using  the  extracted 
values  of  ps  and  p'.  The  values  of  fractal  dimension  (Dr)  of  cancerous  and 
normal  prostate  tissues  were  obtained  by  fitting  to  the  wavelength 
dependence  of  p'.  The  number  of  scattering  particles  contributing  to  ps  as 
a  function  of  particle  size  and  the  cutoff  diameter  </nlax  as  a  function  of  g 
were  investigated  using  the  fractal  soft  tissue  model  and  Mie  theory. 
Results  show  that  rfmax  of  the  normal  tissue  is  larger  than  that  of  the 
cancerous  tissue.  The  cutoff  diameter  £/nlax  is  observed  to  agree  with  the 
nuclear  size  for  the  normal  tissues  and  the  nucleolar  size  for  the  cancerous 
tissues.  Transmission  spectral  polarization  imaging  measurements  were 
performed  that  could  distinguish  the  cancerous  prostate  tissue  samples 
from  the  normal  tissue  samples  based  on  the  differences  of  their 
absorption  and  scattering  parameters. 

Index  Headings:  Prostate  tissue;  Prostate  cancer  detection;  Biomedical 
imaging;  Absorption  coefficient;  Scattering  coefficient;  Reduced  scatter¬ 
ing  coefficient;  Anisotropy  factor;  Near-infrared  imaging;  NIR  spectros¬ 
copy;  Diffusion;  Mie  theory;  Fractal  tissue  model;  Transmission  imaging. 


INTRODUCTION 

Prostate  cancer  is  the  most  common  visceral  malignancy  in 
American  men.  Hence,  there  exists  a  need  for  development  of 
affordable,  noninvasive  modalities  for  early  detection  and 
diagnosis  of  prostate  tumors.  Near-infrared  (NIR)  optical 
imaging  and  spectroscopy  offer  the  potential  for  noninvasive 
detection  and  diagnosis  based  on  differences  in  pathological 
architecture  and  optical  characteristics  between  cancerous  and 
normal  prostate  tissues.  The  disruption  of  the  normal 
architecture  of  tissue  that  accompanies  cancer  and  its 
progression  may  be  modeled  by  fractal  geometry.1  A 
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combination  of  optical  spectroscopic  investigations  and  fractal 
geometric  analysis  presents  an  intriguing  possibility  for 
development  of  a  useful  modality  for  prostate  cancer  detection 
and  diagnosis. 

Variations  in  extracellular2  and  nuclear3  structures  and 
chromophore  compositions  cause  differences  in  optical  scat¬ 
tering  and  absorption  properties  of  cancerous  and  normal 
prostate  tissues,  which  in  turn  provide  optical  contrast,  and  may 
serve  as  native  markers  for  prostate  cancer  detection.  A  survey 
of  the  literature  suggests  that  the  change  of  water  content 
within  living  cells  might  be  the  primary  factor  in  carcinogen¬ 
esis.4-6  Tissues  with  high  water  content  will  usually  be  hypo- 
dense  in  relation  to  normal  tissue.5  The  dominant  chromophore 
in  biological  tissue,  which  absorbs  light  in  the  infrared  range,  is 
water.4,5 

A  systematic  investigation  of  the  absorption  coefficient  (pa), 
the  scattering  coefficient  (ps),  the  reduced  scattering  coefficient 
(p'),  and  the  anisotropy  factor  (g)  of  tissue  is  necessary  for 
accurate  modeling  of  light  transport  in  tissues.7  The  differences 
in  the  values  of  pa,  ps,  p',  and  g  between  normal  and  cancerous 
tissues  may  provide  a  basis  for  biomedical  imaging  and 
potential  diagnostics  using  optical  techniques.7,8 

Several  authors  have  suggested  fractal  analysis  as  a  useful 
tool  in  the  study  of  cancer.9-14  Baish  and  Jain  elucidated  the 
role  of  fractal  geometry  in  describing  the  pathological 
architecture  of  tumors  and  its  potential  for  providing  insight 
into  tumor  growth  mechanism.9  Bizzarri  et  al.  emphasized  that 
the  appropriate  level  of  observation  for  the  study  of 
carcinogenesis  is  tissue  and  suggested  that  fractal  measures 
may  provide  reliable  information  about  the  system  complex¬ 
ity.10  The  fractal  approach  was  used  to  model  soft  tissues  and, 
along  with  the  Mie  scattering  theory,  provided  the  basis  for 
numerical  evaluation  of  the  phase  functions,  the  optical 
scattering  coefficient,  the  reduced  scattering  coefficient,  the 
backscattering  coefficient,  and  the  anisotropy  factor. 11-14 
However,  there  is  a  paucity  of  systematic  and  detailed 
investigations  exploring  the  correlation  between  optical 
properties  and  fractal  characteristics  of  normal  and  cancerous 
tissues  involving  different  types  of  cancers  and  cancers  in 
different  stages  of  progression. 
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Fig.  1.  (a)  Schematic  diagram  of  the  diffuse  reflectance  imaging  experimental  setup;  (b)  NIR  spectral  polarization  imaging  arrangement  (CCD:  charge-coupled 

device,  M:  mirror;  L:  lens,  Pi:  polarizer,  P2:  analyzer,  and  WBF:  wide-band  filter). 


In  this  paper,  we  report  on  the  initial  results  of  a  study  of 
optical  properties  and  fractal  parameters  of  ex  vivo  normal  and 
cancerous  prostate  tissues  and  how  these  parameters  vary 
between  these  two  types  of  tissues.  We  obtained  pa,  ps,  p^,  and 
g  of  cancerous  and  normal  prostate  tissues  in  the  750-860  nm 
spectral  range  using  extinction  measurements  on  thin  samples 
and  diffuse  reflection  measurements  on  thick  samples.  Light 
extinction  by  a  thin  specimen  (thickness  <<  mean  scattering 
length,  /s)  within  the  spectral  range  investigated  here  is 
dominated  by  light  scattering.  Therefore,  measurements  on 
thin  tissue  may  approximate  ps  in  tissue  and  present  a  good 
starting  estimation.  Light  transport  through  a  thick  tissue  slice 
is  usually  modeled  as  a  diffusion  process,  and  diffuse 
reflectance  measurements  provide  information  for  assessment 
of  pa,  and  p;  =  ps(l  -  g),  where  g  is  the  anisotropy  factor.7,8,15 

EXPERIMENTAL  METHODS  AND  SAMPLES 

Prostate  tissue  specimens  were  obtained  from  the  Co¬ 
operation  Human  Tissue  Network  (CHTN),  the  National 
Disease  Research  Interchange  (NDRI),  and  Hackensack 
University  Medical  Center  (HUMC)  under  an  Institutional 
Review  Board  (IRB)  approval  at  the  City  College  of  New  York 
(CCNY).  For  every  cancerous  tissue  sample,  a  normal  tissue 
specimen  from  the  same  patient  was  obtained  for  comparison. 
CHTN,  NDRI,  and  HUMC  provided  pathology  reports  for  the 
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samples.  The  samples  were  received  on  dry  ice  and  were 
neither  chemically  treated  nor  frozen  before  spectroscopic  and 
imaging  measurements.  The  time  elapsed  between  tissue 
resection  and  measurements  varied  from  sample  to  sample 
due  to  different  sample  sources,  but  all  were  within  30  hours. 

Six  pairs  of  different  prostate  tissue  specimens  obtained  from 
six  patients  of  age  between  49  and  64  years  were  investigated. 
The  Gleason  scores  of  the  six  cancerous  tissue  specimens 
varied  from  5  to  8  as  indicated  in  the  pathology  reports.  The 
cancerous  samples  were  carefully  checked  to  determine  the 
hard  parts  to  locate  the  cancerous  regions  before  making  any 
measurements.  Assessment  using  hardness  is  acknowledged  to 
be  a  simple  way  to  locate  regions  of  malignancy.6 

For  extinction  spectral  measurements,  the  cancerous  and 
normal  prostate  tissues  were  cut  to  a  thickness  of  —100  pm  and 
an  area  of  —2  cm  X  3  cm.  The  tissue  samples  were  then  put 
into  100  pm  thick  quartz  cells  and  slightly  compressed  to 
ensure  same  thickness  for  the  extinction  measurements.  The 
extinction  spectra  (optical  density,  O.D.,  vs.  wavelength)  of 
normal  and  cancerous  prostate  tissues  in  the  400  to  2400  nm 
spectral  range  were  measured  using  a  Perkin-Elmer  Lambda-9 
UV/VIS/NIR  Spectrophotometer. 

For  diffuse  reflectance  measurements,  the  cancerous  and 
normal  prostate  tissue  samples  were  cut  into  —2  cm  X  2  cm  X 
1.1  cm  pieces.  The  schematic  diagram  of  the  diffusive 
reflectance  experimental  setup  is  shown  in  Fig.  la.  A 
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Ti :  sapphire  laser  (Spectra-Physics  Tsunami)  beam  tunable  from 
750  to  860  nm  was  used  to  illuminate  the  sample.  The  laser  was 
operated  in  continuous-wave  (cw)  mode  and  had  a  spectral  full 
width  at  half-maximum  bandwidth  less  than  1  nm  The  beam 
passed  through  a  40-cm  focal  length  converging  lens  and 
illuminated  the  sample  at  an  incidence  angle  of  a  =  72.4°. 
Backscattered  light  in  the  direction  normal  to  the  surface  of  the 
sample  was  recorded  by  a  charge-coupled  device  (CCD)  camera. 
Measurements  made  in  10  nm  steps  over  the  750  to  860  nm 
spectral  range  resulted  in  12  sets  of  backscattering  data. 

The  effective  area  recorded  for  the  backscattered  light 
intensity  distribution  used  for  analyzing  is  ~2  mm  X  5  mm. 
Light  reflectance  predominantly  probes  the  volume  of  the 
sample  around  the  point  of  incidence.  The  effective  area 
recorded  for  the  backscattered  light  intensity  distribution  and 
used  for  analysis  is  on  the  order  of  ~2  mm  X  5  mm.  While  the 
light  penetrates  deep  into  the  tissue,  the  major  contribution  to 
backscattered  signal  at  the  detection  angle  comes  mainly  from 
a  depth  of  approximately  1  mm.  The  incident  beam  was  aimed 
at  the  cancer  site.  In  the  case  of  a  tumor  of  size  smaller  than  the 
probed  volume,  the  fitted  optical  parameters  should  be 
regarded  as  the  average  property  of  the  probed  volume 
(including  the  tumor  and  its  surrounding  environment).  The 
cancer  sites  examined  in  the  diffuse  reflectance  measurements 
were  of  sizes  comparable  to  the  effective  area  mentioned 
above,  and  hence  the  probed  optical  properties  approximate 
those  of  the  cancerous  areas. 

A  spectral  polarization  imaging  setup  used  for  transmission 
imaging  measurements  is  schematically  shown  in  Fig.  lb.  A 
white  light  beam  collimated  to  a  diameter  of  ~2  cm  was  passed 
through  wide  band  filters  (WBF)  with  a  typical  bandwidth  of 
60  nm  centered  at  700  nm  and  800  nm  to  select  the  desirable 
spectral  range  for  imaging.  A  polarizer  (Pi)  was  located  in  the 
incident  beam  path  before  the  sample  to  obtain  a  linearly 
polarized  incident  beam  of  light.  The  second  polarizer  (P2)  was 
positioned  in  front  of  the  CCD  camera  for  selecting  the 
detection  polarization  direction  (parallel  or  perpendicular  to  the 
orientation  of  Pi).  A  cooled  silicon  CCD  camera  (Photometric 
CH350L),  which  has  spectral  response  in  the  400  to  1000  nm 
range  and  is  equipped  with  a  50-mm  focal-length  zoom  lens, 
was  used  to  record  the  images  in  transmission  geometry. 

THEORETICAL  FORMALISM 

Determination  of  pa,  ps,  p',  g  from  Optical  Spectroscopic 
Measurements.  Since  pa  <<  ps  in  the  750  to  860  nm  spectral 
range  for  prostate  tissue,  the  extinction  coefficient,  pt  =  ps  +  pa 
~  pSi  and  the  optical  density  (O.D.)  spectrum  of  a  thin  prostate 
tissue  sample  is  dominated  by  scattering.  The  scattering 
coefficient  (ps)  is  given  by: 

O.D. 

ps-m  =  2.303—  (1) 

where  L  is  the  light  path  length  in  the  tissue. 

Backscattered  light  intensities  measured  in  diffuse  reflection 
experiments  for  a  thick  prostate  tissue  sample  with  thickness  of 
~1.1  cm  were  used  to  extract  the  optical  absorption  and  the 
reduced  scattering  coefficients.  The  ordinate  of  (x,  y,  z)  is 
defined  in  Fig.  la.  In  the  case  of  an  obliquely  incident  pencil 
beam,  the  origin  of  the  system  co-ordinate  is  set  on  the  surface 
(z  =  0)  and  the  incident  direction  is  (sin  ot,  0,  —cos  a),  where  a 
is  the  incidence  angle  (a  =  72.4°  in  our  case).  The  Snell’s  law 


predicted  direction  of  the  diffraction  beam  within  the  tissue  is 
(sin  P,  0,  —cos  p),  where  sin  a  =  n  sin  p,  and  n  is  the  index  of 
refraction  of  prostate  tissue  (7?  =  1.35). 16  The  backscattered 
diffuse  intensity  at  position  (. x ,  y,  0)  on  the  surface  of  the 
sample  is  given  by:15 

/tcosP 

ffay)  =^3-(1  +  fleff  Pl)®xP(  fleff  Pi) 

/,cosP  +  2ze 

+ - ^3 - (!  +  Feffpzjexpi-fleffpz)  (2) 

where 

Pi  =  \j ( x  -  /,sinP)2  +  y  +  /t2cos2p, 


p2  =  \/(x-  /tsinP)2  +  y2  +  (2 ze  +  /,cosp)2 


and  peff  =  \/3pa/A;  U  =  (fis)  *s  the  transport  mean  free  path; 
ze  is  the  extrapolation  length  that  depends  on  the  mismatch  of 
the  index  of  refraction  at  the  interface.17  The  extrapolation 
length  is  computed  to  be  ze  =  1.69/t.  Equation  2  is  used  to  fit 
the  measured  diffuse  reflectance  along  a  line  in  the  x  direction 
on  the  tissue  surface  to  obtain  pa  and  /,.  The  anisotropy  factor  g 
is  then  obtained  using  g  =  1  —  (p'/ps). 

Fractal  Soft  Tissue  Model  for  Determination  of  Df  and 
dmax.«  As  many  biological  tissues  have  fractal-like  organization 
and  are  statistically  self-similar,  various  tissue  light  scattering 
models  based  on  fractal  geometry  have  been  proposed.  The 
distribution  of  the  number  density  of  particles  of  diameter  dt  is 
given  by1 1,12 


n(di ) 


riofif  Pf  _  6q0  Df 
ndf  /6  7t  ' 


(3) 


where  Df  is  the  fractal  volume  dimension,  r|o  is  a  scale- 
dependent  constant,  and  n(di)Ad  gives  the  number  density  of 
the  particles  with  diameters  between  <7,  and  d,  +  Ad.  The  value 
of  r|o  can  be  determined  from  the  total  volume  fraction  of  all 
particles,  Tv,  using  equations1112 

r  ^max  /  r  P^max 

Tv  =  %  /  dtD,ddi  =►  Bo  =  TJ  '\\dh 

d dmin  /  '  d:rnr, 


where  dmin  and  dmax  are  diameters  of  the  smallest  and  largest  ?1 
scattering  particles,  respectively. 

The  discrete  particle  model  may  not  be  appropriate  to 
describe  tissue  inhomogeneities,  as  on  a  microscopic  scale  the 
constituents  of  tissue  have  no  clear  boundaries  and  merge  into 
a  quasi-continuum  structure.  As  the  refractive  index  variation 
in  biological  tissue  is  weak,  tissue  is  better  modeled  as  a 
continuous  random  medium  where  light  scattering  is  not  due  to 
the  discontinuities  in  refractive  index  but  is  considered  to  be  a 
consequence  of  weak  random  fluctuations  of  the  dielectric 
permittivity. 

In  such  a  fractal  soft  tissue  model,  the  correlation  function  K(r) 
of  the  random  fluctuation  of  the  refractive  index  is  a  weighted 
average  of  fluctuations  over  a  range  of  correlation  length  /,  that 
depends  on  the  fractal  dimension  Df  and  the  cutoff  length  or  the 
maximum  correlation  length,  <7max.13  Light  scattering  property  of 
the  tissue  is  fully  determined  by  the  power  spectrum  (the  Fourier 
transfonn)  of  R(r).  This  model  has  been  shown  to  be  in  good 
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measured  optical  density  (O.D.)  data  in  the  spectral  range  of  750  nm  to  860  nm  using  Eq.  1  in  the  text;  (b)  the  absorption  coefficient  (pa)  of  cancerous  (solid  square) 
and  normal  (solid  circle)  prostate  tissues;  (c)  the  reduced  scattering  coefficients  (p')  of  cancerous  (solid  square)  and  normal  (solid  circle)  prostate  tissues.  The 
absorption  spectrum  of  water  with  2  cm  thickness  in  the  spectral  range  of  750  to  860  nm  is  shown  at  upper  left  as  a  reference;  and  (d)  the  scattering  anisotropy  factor 
g  of  cancerous  (solid  square)  and  normal  (solid  circle)  prostate  tissues  in  the  spectral  range  of  750  to  860  nm. 


agreement  with  NIR  light  scattering  experimental  results  in 
biological  tissues  and  cell  suspensions.13,14 

The  amplitude  scattering  function  5(9)  at  the  scattering  angle 
0  is  given:13 


|5(0)|2 


£2T|0£D| 


f  —  1  r6— D{ 


2n[l  +  2(1  —  cos  0)x2] 


■2l2 


dx 


(4) 


where  s2  =  4hq8ot2,  k  =  2ti«oA,  «o  is  the  index  of  refraction  of 
the  surrounding  medium,  and  § m  is  the  amplitude  of  the 
relative  refractive  index  fluctuation.  The  anisotropy  factor  (g), 
the  mean  cosine  of  the  scattering  angle,  was  then  found  as13 

(1 +cos20)cos0|5(0)|2  /  y  (1  +  cos20)|5(0)|2 

t£-  /  J dL  2 i? 

(5) 


where  Q  is  the  solid  angle  and  c/Q  =  27tsin0t/0  for  a  scattering 
angle  0.13  The  reduced  scattering  coefficient  in  the  fractal  soft 
tissue  model  follows  a  power  law: 

(6) 
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This  provides  a  simple  procedure  to  estimate  the  fractal 
dimension  Df  from  the  wavelength  dependence  of  p(. 13 

Using  Eqs.  4  and  5,  and  the  known  value  of  <7min,  the 
variation  of  g  as  a  function  of  dmax  can  be  plotted  in  a  g  ~  <7max 
curve.  Using  the  g  ~  dmax  curve,  the  maximum  correlation 
length  dmm  can  be  obtained  from  any  determined  anisotropy 
factor  g. 

EXPERIMENTAL  RESULTS  AND  ANALYSIS 

Determination  of  Optical  Coefficients  of  Cancerous  and 
Normal  Prostate  Tissues.  The  wavelength  dependence  of  ps, 
pa,  p(,  and  g  in  the  750  to  860  nm  spectral  range  is  shown  in 
Figs.  2a,  2b,  2c,  and  2d,  respectively.  Typical  values  of  the 
scattering  coefficients  (Fig.  2a)  of  the  normal  and  cancerous 
prostate  tissues  at  800  nm  were  found  to  be  psnor"'al  =  131.0  ± 
9.5  cm-1  and  pscancer  =  94.9  ±6.1  cm-1,  respectively.  Over 
the  entire  wavelength  range,  the  value  of  ps  for  nonnal  prostate 
tissue  is  larger  than  that  for  cancerous  prostate  tissue, 
indicating  that  normal  tissue  scatters  light  more  strongly  than 
cancerous  tissue.  The  gradual  decrease  in  the  value  of  ps  with 
wavelength  is  consistent  with  Mie  theory  prediction. 

As  Fig.  2b  shows,  the  absorption  coefficient  of  normal 
prostate  tissue  is  higher  than  that  of  the  cancerous  tissue  in  the 


?2 
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750  to  860  nm  spectral  range,  with  typical  values  of  p"ormal  = 
(1.73  ±  0.14)  X  10“ 2  cm- 1  and  p“ncer  =  (1.62  ±  0.13)X10"3 
cm-1,  respectively  at  800  nm.  Comparison  with  the  absorption 
spectrum  of  water  in  the  inset  shows  that  the  absorption  spectra 
of  the  tissues  follow  that  of  water  closely,  and  the  behavior  is 
more  prominent  for  the  absorption  spectrum  of  the  normal 
tissue.  This  result  indicates  that  water  is  the  dominant 
absorptive  tissue  constituent  in  this  spectral  range.  What  is 
even  more  important  is  that  water  content  is  much  higher  in 
nonnal  prostate  tissues  than  in  prostate  tumors,  an  observation 
consistent  with  known  relative  abundance  of  water  in  normal 
and  cancerous  tissues.4,5 

Reduced  scattering  coefficient  (Fig.  2c)  of  normal  prostate 
tissue  decreases  and  that  of  cancerous  prostate  tissues  increases 
monotonically  in  the  750  to  860  nm  range,  while  the  behavior 
is  reversed  for  anisotropy  factor  g.  The  overall  magnitude  of  p' 
is  significantly  higher  for  cancerous  prostate  tissue  than  for 
nonnal  prostate  tissue,  with  typical  values  of  (p')normai  =  4.13 
±  0.21  cm"1  and  (p')cancer  =  5.74  ±  0.19  cm"1  at  800  nm.  The 
observed  lower  values  of  g  for  cancerous  tissue  than  that  for 
nonnal  tissue  (Fig.  2d)  follow  readily  from  the  behavior  of  p' 
because  of  the  simple  algebraic  relation  between  the  two 
parameters.  Typical  values  of  g  at  800  nm  for  normal  and 
cancerous  prostate  tissues  are  estimated  to  be  0.968  ±  0.04  and 
0.939  ±  0.06.  It  is  intriguing  that  the  g  factor  of  normal  tissue 
is  always  larger  than  that  of  cancerous  tissues  in  this  spectral 
range.  Since  ps,  p',  and  g  depend  on  a  number  of  parameters, 
such  as  index  of  refraction  and  scattering  particles’  size  and 
distribution,  detailed  analysis  is  given  in  the  next  section  using 
the  fractal  soft  tissue  model  to  illustrate  the  origin  of  these 
observations. 

Numerical  Evaluation  of  Fractal  Dimension  (Df)  and 
Cutoff  Diameters  of  Scatterers  in  Cancerous  and  Normal 
Prostate  Tissues.  To  study  how  the  absorption  and  scattering 
spectra  reflect  micro-structural  properties  of  cancerous  and 
nonnal  prostate  tissues,  we  analyzed  the  scattering  data  of  these 
two  types  of  tissues  numerically  using  the  fractal  soft  tissue 
model  outlined  in  the  Theory  section.  The  following 
parameters  were  used  in  our  simulation:  the  total  volume 
fraction  of  all  scatterers  Tv  =  0.2, 11,12  the  diameter  of  the 
smallest  particle  <7min  =  50  nm, 11,18  the  index  of  refraction  of 
the  cytoplasm  of  prostate  tissue  (host  medium)  no  =  1.35, 16 
and  the  index  of  refraction  of  prostate  cellular  nuclei  (scattering 
particles)  m  =  1.46. 12-14,16 

The  values  of  fractal  dimension,  Df,  were  obtained  by  fitting 
the  data  shown  in  Fig.  2c  using  Eq.  6,  and  the  results  show  that 
Df  =  1.9  ±  0.04  and  4.4  ±  0.06  for  the  cancerous  and  normal 
prostate  tissues,  respectively. 

The  g  values  of  the  cancerous  and  normal  prostate  tissues  at 
800  nm  as  a  function  of  cutoff  diameter  <7max  were  obtained 
using  Eqs.  4  and  5,  and  the  results  are  displayed  as  solid  and 
dashed  lines,  respectively,  in  Fig.  3a.  Using  the  previously 
determined  values  of  g  in  Fig.  2d,  the  cutoff  diameters  of  the 
cancerous  and  normal  prostate  tissues  were  found  to  be  1 .42  ± 
0.06  pm  for  the  cancerous  tissue  and  6.48  ±  0.08  pm  for  the 
normal  tissue. 

The  distribution  of  the  number  density  as  a  function  of  the 
diameter  d  of  scattering  particles  (scatterers)  in  cancerous  and 
normal  prostate  tissues  was  then  calculated  using  Eq.  3,  and  the 
results  are  shown  as  Fig.  3b.  The  salient  feature  of  Fig.  3b  is 
that  the  normal  tissue  has  many  more  smaller  scatterers  than 
the  cancerous  tissue  in  the  d  <  —0.6  pm  range  while  the 


Fig.  3.  (a)  The  anisotropy  factor  g  of  the  cancerous  and  normal  prostate 

tissues  as  a  function  of  cutoff  diameter  dmax  at  800  nm,  and  (b)  distribution  of 
the  number  density  of  scatterers  with  different  diameters,  evaluated  by  fractal 
tissue  model  using  Eq.  3.  The  curves  representing  the  behavior  of  cancerous 
and  normal  prostate  tissues  are  displayed  by  solid  and  dashed  lines, 
respectively.  The  cancerous  prostate  tissue  (solid  line)  has  lower  weight  of 
smaller  particles  and  higher  weight  of  larger  particles  than  the  normal  prostate 
tissue  (dashed  line). 

cancerous  tissue  has  more  larger  scatterers  than  the  normal 
tissue  in  the  d  >  —0.6  pm  range. 

The  fractal  parameters  for  the  normal  and  cancerous  prostate 
tissues  are  summarized  in  Table  I. 

NEAR-INFRARED  OPTICAL  POLARIZATION 
IMAGING 

The  efficacy  of  near-infrared  (NIR)  optical  imaging  as  a 
detection  and  diagnostic  tool  for  prostate  cancer  will  depend  on 
the  contrast  between  normal  and  cancerous  prostate  tissues, 
which  arise  from  differences  in  intrinsic  optical  spectroscopic 
properties  when  no  external  agents  are  used.  The  differences  in 
optical  scattering  coefficient  and  anisotropy  factor  reported  in 
the  previous  section  prompted  us  to  explore  the  potential  of 
NIR  imaging  in  prostate  cancer  detection  using  a  model 
prostate  tissue  sample  consisting  of  a  piece  of  cancerous  tissue 
on  the  right  side  and  a  piece  of  normal  tissue  on  the  left  side. 
The  NIR  spectral  polarization  imaging  system  is  schematically 
shown  in  Fig.  lb.  The  sample  with  dimensions  of  —2  cm  X  3 

APPLIED  SPECTROSCOPY 


//Xinet/production/a/apls/liveJobs/apls-66-07/apls-66-07-09/layouts/apls-66-07-09.3d  ■  Wednesday,  25  April  2012  ■  9:46  am  ■  Allen  Press,  Inc.  ■  Page  5 


TABLE  I.  Df  and  dmax  for  cancerous  and  normal  prostate  tissues. 


Tissue  type 

Df 

rfmax  (pm) 

Cancer 

1.9  ±  0.04 

1.42  ±  0.06 

Normal 

4.4  ±  0.06 

6.48  ±  0.08 

cm  X  0.035  cm  was  held  between  two  glass  plates  under  slight 
compression  to  ensure  uniform  thickness.  Figures  4a  and  4b 
display  the  polarization-gated  images  of  the  sample  recorded 
with  transmitted  light  polarized  parallel  and  perpendicular  to 
polarization  of  the  incident  light,  respectively.  The  center 
wavelength  of  the  imaging  beam  was  800  nm.  It  can  be  seen 
from  Fig.  4  that  the  image  of  the  right  side  of  the  sample 
(consisting  primarily  of  cancerous  tissue)  is  brighter  than  that 
of  the  left  side  (consisting  mainly  of  normal  tissue).  This  can  be 
explained  by  the  differences  of  scattering  and  absorption 
coefficients  between  cancerous  and  normal  prostate  tissues. 
Since  both  the  scattering  and  absorption  coefficients  of  the 
cancerous  tissue  at  800  nm  are  smaller  than  the  corresponding 
coefficients  of  the  normal  tissue,  as  shown  in  Figs.  2a  and  2b, 
the  light  propagating  in  the  cancerous  tissue  undergoes  less 
scattering  and  absorption.  As  a  result,  the  light  transmitted 
through  the  cancerous  tissue  area  is  stronger  than  that  through 
the  normal  tissue. 

The  differences  between  the  intensities  transmitted  through  the 
cancerous  and  normal  prostate  tissue  areas  can  be  seen  more 
clearly  from  the  spatial  intensity  distributions  obtained  by 
integrating  the  image  intensity  over  a  horizontal  rectangular 
region  as  marked  with  the  dashed  boxes  in  Figs.  4a  and  4b.  The 
integrated  image  intensity  distributions  for  the  images  shown  in 
Figs.  4a  and  4b  are  displayed  in  Figs.  4c  and  4d,  respectively.  The 
average  image  intensities  of  the  cancerous  tissue  area  were  found 
to  be  1 .04  and  1 .7 1  times  stronger  than  that  of  the  normal  tissue 
area  for  the  parallel  and  perpendicular  images,  respectively. 

The  image  intensity  distributions  were  also  used  to  calculate 
the  degree  of  polarization,  defined  as  F  =  (/  -  I±)/(I\\  +  7j_), 
where  7||  and  /  are  intensities  of  transmitted  light  parallel  and 
perpendicular  to  the  polarization  of  the  incident  light  beam, 
respectively.  The  value  of  P  for  the  normal  and  cancerous 
prostate  tissues  at  800  nm  are  Pcancer  =  0.24  and  Fnorma|  = 
0.46.  Similar  higher  values  of  P  for  normal  tissue  than  for 
cancerous  tissue  are  observed  for  images  with  other  wave¬ 
lengths  in  the  760  to  850  nm  range.  This  is  expected  since  the  g 
factor  of  normal  prostate  tissue  is  larger  than  that  of  the 
cancerous  tissue,  as  shown  in  Fig.  2d.  Linearly  polarized  light 
propagating  in  a  more  strongly  forward-peaked  scattering 
medium  (characterized  by  a  larger  value  of  g)  is  less 
depolarized  in  one  scattering  event.  The  smaller  g  factor  for 
the  cancerous  tissue  results  in  more  depolarization  and  a 
smaller  value  of  P  in  comparison  with  that  in  the  normal 
tissue.19 

DISCUSSION 

In  this  paper  we  have  presented  the  initial  results  of  a  study 
to  explore  potential  correlation  between  optical  properties  and 
fractal  dimensional  parameters  in  normal  and  cancerous  human 
prostate  tissues.  The  extinction  spectra  of  cancerous  and 
normal  prostate  tissues  in  the  750  to  860  nm  spectral  range 
were  used  to  obtain  the  scattering  coefficient  (ps).  The 
absorption  coefficient  (pa)  and  the  reduced  scattering  coeffi- 
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Fig.  4.  Transmission  images  of  a  normal  (N)  and  cancerous  (C)  human 
prostate  tissue  sample,  obtained  using  800  nm  light  for  (a)  parallel  and  (b) 
perpendicular  polarization  configurations.  The  size  of  tissue  for  imaging  is  — 2 
cm  X  3  cm  X  0.035  cm  (length  X  height  X  depth),  (c,  d)  The  digitized  spatial 
intensity  distributions  of  images  (a)  and  (b),  respectively,  obtained  by 
integrating  the  image  intensity  over  the  horizontal  rectangular  region  as 
marked  with  the  dashed  boxes  in  (a)  and  (b). 


cient  (p')  were  acquired  by  fitting  the  diffuse  reflectance  of  an 
obliquely  incident  beam  using  the  diffusion  model  in  a  thick 
tissue  sample.  The  values  of  p'  were  found  to  be  greater  than  pa 
by  2  to  3  orders  of  magnitude  in  this  spectral  range,  and  the 
results  are  consistent  with  literature  values  for  prostate  tissue.7-8 
The  anisotropy  factor  g  was  calculated  from  ps  and  p'. 

A  fractal  soft  tissue  model  was  used  to  extract  the  values  of 
fractal  dimension  (Df)  from  the  wavelength  dependence  of  p' 
and  were  found  to  be  1 .9  and  4.4  for  the  cancerous  and  normal 
prostate  tissues,  respectively.  The  values  are  markedly 
different.  Similar  remarkable  difference  was  observed  in  the 
values  for  cutoff  diameter,  <7max:  1.42  ±  0.06  pm  for  the 
cancerous  tissue  and  6.48  ±  0.08  pm  for  the  normal  tissue. 
Since  the  cutoff  diameter  may  be  considered  to  be  an  upper 
limit  on  the  size  of  scatterers  in  tissue  specimens,  this  result 
may  be  quite  significant.  Further  probing  revealed  that  the 
observed  value  of  the  cutoff  diameter  of  the  cancerous  tissue 
(1.42  ±  0.06  pm)  falls  within  the  range  of  1.38-1.96  pm  for 
the  nucleolar  diameter  of  cancerous  prostate  cell,20  and  the 
cutoff  diameter  for  normal  tissue  (6.48  ±  0.08  pm)  closely 
matches  the  nuclear  diameter  of  normal  prostate  cell  (6.3-6.72 
pm).20  These  close  agreements  point  to  the  intriguing 
implication  that  in  normal  prostate  cells  the  largest  scatterers 
are  the  nuclei,  while  in  cancerous  prostate  cells  nucleoli  are  the 
largest  scatterers.  However,  such  a  hypothesis  raises  the 
question  why  it  is  that  the  nucleoli  represent  the  upper  limit 
to  the  size  of  scatterers  in  cancerous  cells,  while  nuclei 
determine  the  upper  limit  in  normal  cells.  In  this  context  it 
should  be  noted  that  the  diameter  of  the  normal  prostate 
nucleoli  he  in  the  1.10  to  1.24  pm  range,  and  that  of  cancerous 
nuclei  is  in  the  8.4  to  10.31  pm  range.20 
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A  quantitative  answer  to  this  question  is  beyond  the  scope  of 
this  paper.  However  some  qualitative  and  plausible  observa¬ 
tions  based  on  the  morphological  changes  associated  with 
cancer  and  consequent  changes  in  optical  spectroscopic 
characteristics  can  be  made  in  support  of  the  hypothesis.20-25 
Normal  cells,  in  general,  are  characterized  by  uniformity  in  size 
and  shape  of  nuclei  (round  or  oval)  and  have  an  ordered 
morphological  structure.21  The  cell  nucleus  is  dense  and 
compact  and  may  be  considered  as  an  entity  with  a  localized 
index  of  refraction  inhomogeneity.  So,  the  entire  nucleus  acts 
as  a  scattering  entity  that  happens  to  be  larger  than  any  other 
scatterer  within  the  normal  cell.  This  neat,  regular,  and  ordered 
structure  breaks  down  in  cancerous  cells.  Malignant  cells  are 
characterized  by  irregular  morphology,  including  changes  in 
the  shape  and  size  of  the  nucleus,  the  nuclear  envelope,  the 
nucleolus,  and  the  chromatin.21  Richards-Kortum  et  al.  found 
that  light  scattering  is  strongly  influenced  by  the  difference  in 
morphology  of  nucleus  and  nucleolus,  DNA  content,  and 
clumping  and  coarse  chromatins.22  Fischer  et  al.  show  that  two 
important  cell  structure  changes  in  human  prostate  cancer  are 
nucleolar  enlargement  and  nuclear  envelope  irregularity.23 
Varma  et  al.  investigated  the  morphologic  criteria  for  diagnosis 
of  prostate  adenocarcinoma  in  250  consecutive  cases  and 
observed  that  dysplastic  nuclei  have  irregular  size  and  shape, 
clumps  of  condensed  chromatin,  and  multiple  nucleoli.24  They 
suggest  that  the  prominent  nucleoli,  marginated  nucleoli 
(location  of  nucleoli  is  close  to  or  at  the  inner  nuclear 
membrane)  and  multiple  nucleoli  are  criteria  having  highest 
sensitivity  and  specificity  for  diagnosis  of  prostate  cancer.24  To 
sum  up,  because  of  high  density  of  malignant  prostate  cells, 
irregularity  of  nuclear  envelope,  enlargement  of  nucleoli, 
clumping  and  condensation  of  the  chromatin  and  marginated 
nucleoli,22-25  and  fragmentation  of  nuclei  and  nuclear 
degeneration,  it  seems  that  nuclei  no  longer  behave  as  the 
largest  compact  entities  in  cancerous  cells  when  it  comes  to 
light  scattering:  that  role  is  taken  up  by  nucleoli. 

Fractal  model  analysis  of  optical  spectroscopic  measure¬ 
ments  on  normal  prostate  tissues  and  prostate  tumors  at 
different  stages  of  growth  and  different  Gleason  grades 
corroborated  by  microscopic  and  histological  examinations  of 
the  samples  will  be  needed  to  put  the  hypothesis  presented 
above  to  stringent  test.  Affirmation  of  the  hypothesis  would 
imply  that  the  fractal  dimension  and  the  cutoff  diameter 
predicted  by  the  fractal  soft  tissue  model  may  be  used  as 
criteria  for  prostate  cancer  detection.  The  results  further 
indicate  that  the  NIR  spectral  polarization  imaging26  has  the 
potential  to  be  developed  as  a  modality  for  prostate  cancer 
detection. 
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Stokes  shift  spectroscopy  (S3)  is  an  emerging  approach  for  cancer  detection.  The  goal  of  this  paper  is  to 
evaluate  the  diagnostic  potential  of  the  S3  technique  for  the  detection  and  characterization  of  normal  and 
cancerous  prostate  tissues.  Pairs  of  cancerous  and  normal  prostate  tissue  samples  were  taken  from  each 
of  eight  patients.  Stokes  shift  spectra  were  measured  by  simultaneously  scanning  both  the  excitation  and 
emission  wavelengths  while  keeping  a  fixed  wavelength  interval  Ai  =  20  nm  between  them.  The  salient 
features  of  this  technique  are  the  highly  resolved  emission  peaks  and  significant  spectral  differences 
between  the  normal  and  cancerous  prostate  tissues,  as  observed  in  the  wavelength  region  of  250  to 
600  nm.  The  Stokes  shift  spectra  of  cancerous  and  normal  prostate  tissues  revealed  distinct  peaks  around 
300,  345,  440,  and  510  nm,  which  are  attributed  to  tryptophan,  collagen,  NADH,  and  flavin,  respectively. 
To  quantify  the  spectral  differences  between  the  normal  and  cancerous  prostate  tissues,  two  spectral 
ratios  were  computed.  The  findings  revealed  that  both  ratio  parameters  f?!  =^297/^345  an(f  -^2  = 
/307//345  were  excellent  diagnostic  ratio  parameters  giving  100%  specificity  and  100%  sensitivity  for 
distinguishing  cancerous  tissues  from  the  normal  tissue.  Our  results  demonstrate  that  S3  is  a  sensitive 
and  specific  technique  for  detecting  cancerous  prostate  tissue.  ©  2012  Optical  Society  of  America 


1.  Introduction 

Prostate  cancer  is  the  second  leading  cause  of  cancer 
death  among  American  men.  During  2010,  an  esti¬ 
mated  217,730  new  prostate  cancer  cases  were  re¬ 
ported  in  the  United  States,  and  it  was  predicted 
that  among  them,  32,050  men  might  die  [1].  Cur¬ 
rently,  prostate  cancer  diagnosis  is  based  on  a  pros¬ 
tate  specific  antigen  (PSA)  test  and  digital  rectal 
examination  (DRE).  If  either  of  these  two  tests  is  ab¬ 
normal,  a  biopsy  is  usually  performed  guided  by 
transrectal  ultrasound  (TRUS).  Although  a  TRUS- 
guided  biopsy  is  considered  the  gold  standard,  it  suf¬ 
fers  from  lack  of  sensitivity  and  specificity  and  leads 
to  a  significant  number  of  false  negatives,  which  then 
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often  leads  to  unnecessary  biopsies,  patient  trauma, 
and  the  time  needed  to  obtain  a  histopathological  di¬ 
agnosis  [2],  Consequently,  new  detection  technolo¬ 
gies  are  needed  that  can  overcome  the  current 
limitations  and  improve  patient  well-being. 

Optical  spectroscopy  offers  a  novel  diagnostic  ap¬ 
proach  for  tissue  and  can  be  considered  as  an  alterna¬ 
tive  technique  to  conventional  methods  because  of  its 
advantages,  such  as  minimal  invasiveness,  less  time, 
and  reproducibility.  For  more  than  three  decades, 
various  optical  spectroscopic  techniques  for  cancer 
diagnosis  have  been  widely  explored  as  potential  di¬ 
agnostic  tools  in  the  discrimination  of  normal  from 
abnormal  tissues.  In  1980  s,  Alfano’s  group  first 
reported  on  the  light-induced  fluorescence  differences 
from  malignant  and  nonmalignant  breast  and 
lung  tissues  [3],  Since  then,  some  groups  around 
the  world  have  reported  the  use  of  native  fluorescence 
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spectroscopy  [4,5]  and/or  diffuse  reflectance  spectro¬ 
scopy  [6,7]  for  diagnosis  of  cancers  in  different  organs. 
These  reported  studies  indicate  that  optical  spectro¬ 
scopy  has  the  potential  to  improve  the  screening  and 
early  detection  of  cancer,  the  optical  spectra  method, 
known  as  “optical  biopsy,”  as  termed  by  Alfano,  is  a 
potential  technique  by  which  the  medical  community 
can  diagnose  tissues  without  removing  them. 

Native  fluorescence  of  tissues,  also  called  auto¬ 
fluorescence,  is  believed  to  be  produced  by  several  en¬ 
dogenous  fluorophores,  such  as  tryptophan,  collagen, 
elastin,  a  reduced  form  of  nicotinamide  adenine 
dinucleotide  (NADH),  flavin  adenine  dinucleotide 
(FAD),  and  endogenous  porphyrins.  During  the  de¬ 
velopment  from  benign  prostatic  hyperplasia  to 
premalignant  (dysplastic)  and  malignant  stages, 
prostate  cells  undergo  proliferations  that  modify 
these  fluorophore  levels  [8].  Additionally,  the  connec¬ 
tive  tissue  frameworks  of  prostate  tissue  can  be  im¬ 
paired  during  cancer  development  [9].  Such  changes 
are  likely  to  alter  both  tissue  morphology  and  bio¬ 
chemistry,  which  could  be  detected  using  tissue  fluor¬ 
escence  spectroscopy  [10].  Thus,  by  measuring  the 
fluorescence  signal  from  tissue,  changes  in  tissue 
structure  and  composition  and  information  concern¬ 
ing  the  pathological  state  of  the  tissue  could  be 
obtained  [11,12]. 

Cells/tissues  contain  several  key  fluorophores  with 
broad  and  partly  overlapped  excitation  and  emission 
spectra  [R3] .  Further,  each  fluorophore  has  a  distinct 
absorption  peak  wavelength  and  a  characteristic 
“fingerprint”  peak  emission  wavelength  [14].  The 
emission  spectra  at  one  or  more  excitation  wave¬ 
lengths  or  excitation  spectra  corresponding  to  one 
or  more  emission  wavelengths  have  been  used  for 
diagnostic  purposes  [15].  The  conventional  laser- 
induced  fluorescence  spectroscopic  method  for  cancer 
detection  has  limited  applicability  since  most  tissue 
fluorescence  spectra  have  a  series  of  overlapping 
bands  from  different  fluorophores,  and  it  often  can¬ 
not  be  resolved  satisfactorily  by  a  single-excitation 
wavelength.  To  overcome  this  problem,  multiple-ex¬ 
citation  wavelengths  are  sequentially  used  to  gener¬ 
ate  an  excitation-emission  matrix  (EEM)  in  order  to 
identify  the  optimal  excitation  wavelengths  at  which 
tissue  classification  is  enhanced  and  to  determine 
the  origin  of  the  measured  fluorescence  signal  in  a 
more  reliable  manner  [5].  However,  EEM  measure¬ 
ment  requires  a  large  number  of  fluorescence  emis¬ 
sion  scans  at  sequential  excitation  wavelengths  at 
small  wavelength  intervals,  and  this  is  time-  con¬ 
suming.  In  order  to  reduce  time  consumption,  a  two- 
dimensional  detector,  such  as  a  CCD  camera,  is  used 
in  conjunction  with  a  spectrograph. 

Alfano  and  Yang  were  the  first  to  introduce  the 
physical  concept  of  Stokes  shift  spectroscopy  (S3) 
as  a  rapid  method  for  diagnosing  diseased  tissue  that 
depends  on  both  absorption  and  emission  properties 
of  the  fluorophores  [16],  In  S3,  the  excitation  wave¬ 
length  Aex  and  the  emission  wavelength  lem  are 
scanned  synchronously  with  a  constant  wavelength 


interval  A  A  =  /<,m— /ex  between  the  excitation  and 
emission.  Earlier,  this  technique  was  called  synchro¬ 
nous  fluorescence  spectroscopy,  and  it  begins  with  a 
multicomponent  analysis  to  obtain  spectral  emission 
peaks  of  tissue  samples,  and  it  enhances  the  selectiv¬ 
ity  in  the  assay  of  complex  systems  [17].  Recently,  S3 
was  used  as  a  potential  tool  for  the  diagnostic  pur¬ 
poses.  Ebenezar  et  al.  investigated  the  potential 
use  of  S3  as  a  diagnostic  tool  for  the  discrimination 
between  normal  and  abnormal  cervical  and  breast 
tissues  [18,19].  Masilamani  et  al.  studied  the  use 
of  S3  in  the  discrimination  of  benign  and  malignant 
prostate  tissues  [20].  However,  to  the  best  of  our 
knowledge,  no  report  is  available  using  S3  technique 
for  the  comparison  and  discrimination  of  cancerous 
and  normal  prostate  tissues  from  the  same  patients. 
Therefore,  this  paper  is  designed  to  demonstrate  the 
potential  use  of  the  S3  technique  to  detect  and  char¬ 
acterize  the  normal  and  cancerous  prostate  tissues 
from  the  same  patients,  and  to  determine  the  optimal 
offset  wavelength  A k  for  distinguishing  cancerous 
tissue  from  normal  prostate  tissue. 

2.  Materials  and  methods 

A.  Tissue  Samples 

A  total  of  sixteen  fresh  human  prostate  tissue  sam¬ 
ples  were  collected  from  eight  patients  from  the  Co¬ 
operation  Human  Tissue  Network  (CHTN)  under  the 
IRB  approval  of  The  City  College  of  New  York.  Each  Q 
patient  provided  a  pair  of  normal  and  cancerous  sam¬ 
ples.  The  spectroscopic  data  were  classified  into  two 
groups:  (1)  normal  and  (2)  cancerous  tissues  accord¬ 
ing  to  their  post-spectroscopy  studies.  The  optical 
biopsy  of  each  tissue  sample  was  compared  to  the  his- 
topathology  analyses. 

B.  S3  Measurements 

The  Stokes  shift  (SS)  spectra  of  ex  vivo  prostate  tis¬ 
sues  were  recorded  using  spectrofluorometer  (LS  50, 
PerkinElmer).  The  prostate  tissue  samples  of  solid 
chunks  of  approximately  1.5  mm  x  1  mm  x  0.3  mm 
(length  x  width  x  thickness)  were  placed  in  a  quartz 
cuvette  with  the  epithelium  toward  the  face  of  the 
cuvette  and  the  beam.  The  excitation  light  of  ap¬ 
proximately  2  mm  x  4  mm  was  incident  perpendicu¬ 
lar  to  the  face  of  the  tissue  epithelium  surface,  and 
the  emitted  fluorescence  light  was  collected  at  a 
90°  angle  to  the  excitation  light.  During  data  acqui¬ 
sition,  the  excitation  and  emission  monochromators 
had  fixed  band  passes  of  3  nm  each,  and  the  wave¬ 
length  increment  was  set  at  0.5  nm.  The  variation 
in  the  excitation  light  source  intensity  as  a  function 
of  wavelength  was  taken  into  account  during  SS 
spectra  measurements.  This  was  done  by  detecting 
the  fluorescence  signal  (S)  by  the  PMT  as  well  as  re-  @ 
cording  the  reference  excitation  intensity  (R)  by  a 
photodiode  and  taking  their  ratio  S/R  to  serve  as 
the  final  S3  signal  to  account  for  the  wavelength- 
dependent  light  intensity. 
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C.  Optimization  of  Offset  Wavelength  A/l 

The  SS  is  dependent  upon  the  vibrational  and  elec¬ 
tronic  interaction  and  the  polarity  of  the  host  envir¬ 
onment  surrounding  the  emitting  organic  molecules. 
The  wavelength  interval  of  the  SS  spectrum  can  be 
simply  compressed  or  expanded  just  by  decreasing 
or  increasing  the  selected  A A  parameter  [19] .  During 
S3  acquisition,  the  excitation  and  emission  mono¬ 
chromators  were  scanned  simultaneously  at  the 
same  speed  with  a  constant  wavelength  interval  A/ 
between  them.  In  this  paper,  we  investigated  various 
AA  value,  ranging  from  10  to  50  nm,  in  increments  of 
10  nm.  Among  the  five  AA  alues,  a  A  A  =  20  nm  re¬ 
vealed  four  key  identifiable  fluorophores  with  a  good 
signal-to-  noise  ratio.  Therefore,  the  value  of  A  A  = 
20  nm  was  selected  as  the  favored  offset  wavelength 
interval  for  the  entire  series  of  investigation.  The  S3 
of  normal  and  cancerous  prostrate  tissues  were 
recorded  in  wavelengths  of  250  to  600  nm,  where 
AA  =  20  nm,  and  with  a  scanning  speed  of 
5  nm/sec.  Because  the  S3  measurements  involved 
the  simultaneous  scanning  of  both  excitation  and 
emission  monochromators,  significant  contribution 
from  the  excitation  light  source  was  expected.  By  di¬ 
viding  each  of  the  SS  spectra  by  the  lamp  spectrum 
measured  under  the  same  instrumental  parameters, 
the  contribution  from  the  excitation  light  source  to 
the  SS  spectra  was  counted. 

D.  Statistical  Analysis 

The  measured  spectral  data  were  analyzed  statisti¬ 
cally  to  discriminate  cancerous  tissue  from  normal 
prostate  tissue.  Data  were  initially  preprocessed  by 
normalizing  each  spectral  profile  with  respect  to  the 
maximum  fluorescence  intensity  of  spectrum.  From 
the  normalized  spectrum,  two  intensity  ratios  were 
calculated  at  different  wavelengths  corresponding 
to  the  characteristic  spectral  features  of  different 
groups  of  the  prostate  tissues  studied.  The  mean  and 
standard  deviation  values  were  calculated  for  each 
group  of  prostate  tissue,  and  their  statistical  signifi¬ 
cance  was  verified  using  an  unpaired  Student’s  t-test 
to  calculate  the  level  of  physical  significance  p  with 
95%  confidence  interval  using  origin  statistical  soft¬ 
ware  (OriginPro7).  The  corresponding  tissues,  which 
yielded  significant  differences  (p  <  0.005)  in  the  ratio 
values,  were  selected  for  further  histopathological 
analyses. 

The  classification  analysis  was  determined  on  the 
basis  of  discrimination  cutoff  values.  Discrimination 
cutoff  lines  are  drawn  between  the  normal  (mean  of 
all  normal  tissues)  and  cancerous  (mean  of  all  cancer 
tissues),  at  values  that  correspond  to  the  mean  ratio 
value  of  the  respective  groups.  The  sensitivity  and 
specificity  in  discriminating  each  of  these  categories 
were  determined  on  the  basis  of  the  cutoff  values 
identified  by  ratio  method  by  validation  with  the  gold 
standard,  namely,  histopathological  results  of  biopsy 
specimens  taken  from  SS  emission  measurement 
sites.  Similar  ratio  intensity  variable  methods  have 
been  adopted  by  earlier  researchers  to  discriminate 


normal  from  different  tissue  categories  of  bladder 
tissues,  and  the  sensitivity  and  specificity  were  de¬ 
termined  on  the  basis  of  the  discrimination  cutoff 
lines  [5], 

3.  Results 

The  normalized  mean  S3  of  cancerous  (solid)  and  nor¬ 
mal  (dash)  prostate  tissues  in  the  wavelength  region 
between  250  and  600  nm  is  shown  in  Fig.  1.  The  nor¬ 
mal  tissues  show  a  small  hump  around  of  307  nm, 
with  a  single  primary  band  centered  at  345  nm,  a 
small  shoulder  around  375  nm,  and  secondary  peaks 
observed  at  405,  440  and  510  nm.  Although  cancer¬ 
ous  tissue  is  found  to  exhibit  two  primary,  highly  re¬ 
solved  bands  at  297  and  345  nm,  with  secondary 
peaks  observed  at  440  and  510  nm,  and  the  spectral 
band  at  297  nm  was  higher  than  that  at  345  nm. 
From  this  normalized  spectrum,  it  was  also  observed 
that  all  the  spectral  bands  are  highly  resolved,  and 
one  could  discern  the  spectral  shift  of  primary  emis¬ 
sion  bands  between  normal  and  cancerous  tissues. 
For  example,  the  normalized  spectra  of  normal 
tissues  show  an  emission  band  at  307  nm,  which 
was  blueshifted  to  297  nm  for  cancerous  tissues.  In 
addition,  a  spectral  dip  or  valley  was  also  observed 
around  425  nm  for  normal  and  cancerous  tissues. 
It  was  noted  that  the  valley  at  425  nm  is  deeper 
in  the  case  of  cancerous  than  normal  tissues. 

A.  S3  Measurements  of  Standard  Fluorophores 

To  investigate  the  origin  of  the  observed  spectral 
peaks  from  normal  and  cancerous  prostate  tissues 
and  to  assign  them  to  various  endogenous  fluoro¬ 
phores,  the  S3  of  commercially  obtained  fluorophores 
(Sigma  Chemicals,  USA)  of  tryptophan,  collagen, 
NADH,  and  FAD  were  measured  with  the  same 
experimental  parameters  maintained  as  for  the  pros¬ 
tate  tissues.  The  SS  spectra  of  the  commercial-grade 
fluorophores  were  corrected  for  the  variation  in  the 
lamp  intensity  and  normalized  with  respect  to  the 
corresponding  peak  intensity.  Figure  2  shows  the 
normalized  SS  spectra  of  the  standard  fluorophores, 


Fig.  1.  (Color  online)  Normalized  mean  Stokes  shift  spectra  of 
cancerous  (solid)  and  normal  (dash)  prostate  tissues. 
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Fig.  2.  (Color  online)  Normalized  Stokes  shift  spectra  of  standard 
fluorophores  of  tryptophan  (solid),  elastin  (dash),  collagen  (dot), 
NADH  (dash  dot),  and  FAD  (dash  dot  dot). 


such  as  tryptophan  (solid),  elastin  (dash),  collagen 
(dot),  NADH  (dash  dot)  and  FAD  (dash  dot  dot),  with 
their  narrow  emission  peaks  observed  at  311,  347, 
353,  397  and  533  nm,  respectively.  The  SS  spectra 
of  elastin  also  showed  secondary  peaks  at  285  nm, 
with  a  small  shoulder  around  414  nm.  Further,  the 
fine  structures  observed  in  the  wavelength  region  be¬ 
tween  440  and  500  nm  may  be  due  to  the  stray  light 
contribution  from  the  xenon  lamp  passing  through 
the  exit  slit  of  the  excitation  monochromator.  A  sharp 
spectral  feature  was  observed  at  467  nm  due  to  the 
contribution  from  the  xenon  lamp  source. 

B.  Results  of  Statistical  Analysis 

In  order  to  validate  the  diagnostic  utility  of  the  ob¬ 
served  spectral  signatures  between  normal  and  can¬ 
cerous  tissues,  two  ratio  variables  were  introduced 
for  SS  spectra,  with  Ai  =  20  nm.  The  mean  and 
standard  deviation  of  the  SS  spectral  ratios  of  nor¬ 
mal  and  cancerous  prostate  tissues  were  calculated, 
and  the  results  are  listed  in  Table  1  along  with  p  va¬ 
lues.  The  statistical  significance  was  computed 
through  an  unpaired  Student’s  t-test  to  determine 
the  level  of  significance  p.  The  p  values  for  both 
the  ratio  variables  are  less  than  0.005,  indicating  a 
high  statistical  significance.  From  Table  1,  for  both 
ratio  variables  Rx  =  Amy/Tus  and  R2  =  /307//345, 
the  mean  value  of  the  normal  tissue  was  lower  than 
that  of  the  cancerous  tissue.  In  the  SS  spectrum,  it  is 
the  relative  intensities  between  the  maxima  that  are 
of  importance  rather  than  the  numerical  intensity 
value  of  each  maximum.  Table  1  shows  the  computed 
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Fig.  3.  Scatter  plot  of  intensity  ratio  of  /297A345  f°r  normal  (•) 
and  cancerous  (■)  prostate  tissues. 


ratios  of  the  fluorescence  peak  intensities  at  297  to 
345  nm  (Ri),  and  307  to  345  nm  (R2),  correlating  with 
their  histologic  findings.  For  example,  the  mean  ratio 
value  of  Ri  in  normal  tissues  was  0.36  ±  0.05,  which 
is  significantly  different  from  the  mean  value  of 
1.24  ±  0.09  in  cancerous  tissues  (unpaired  Student’s 
t-test,  p  <  0.005).  By  selecting  the  I2^/I^  of  ~0.80 
as  the  demarcation  ratio  level  (between  the  normal 
and  cancerous  tissues),  it  is  found  to  yield  100%  spe¬ 
cificity  and  100%  sensitivity  for  differentiating  can¬ 
cerous  from  normal  tissues  (Fig.  3).  Similarly,  the 
mean  ratio  value  of  R2  for  normal  tissue  was 
0.57  ±  0.05,  which  is  also  significantly  different  from 
the  mean  value  of  0.92  ±  0.05  for  cancerous  prostate 
tissues  (unpaired  Student’s  t-test,  p  <  0.005).  By  se¬ 
lecting  /307//345  of  ~0.74  as  a  cutoff  for  tissue  demar¬ 
cation,  100%  specificity  and  100%  sensitivity  were 
achieved  (Fig.  4). 

4.  Discussion 

The  potential  use  of  the  S3  technique  in  detection 
and  characterization  of  normal  and  cancerous  pros¬ 
tate  tissues  was  explored,  and  the  spectral  data  were 
analyzed  by  a  simple  statistical  method,  which  vali¬ 
dated  the  diagnostic  potentiality  of  the  present  meth¬ 
od.  The  comparison  of  S3  spectra  of  normal  and 
cancerous  tissues  with  those  corresponding  to  com¬ 
mercial-grade  standard  fluorophores  suggests  that 
the  peaks  observed  around  297,  307,  345,  440,  and 
510  nm  may  be  primarily  attributed  to  tryptophan, 
collagen  and/or  elastin,  NADH,  and  FAD  respec¬ 
tively.  Similar  peaks  of  S3  were  reported  earlier  in 
the  discrimination  of  normal  and  different  pathologi¬ 
cal  studies  of  cervical  tissues  [18].  In  this  paper,  the 


Table  1.  Mean  ±SD  of  the  Stokes  Shift  Spectral  Ratios  of  Normal  and  Cancerous  Prostate  Tissues 


Spectral  Intensity 


Tissue  types 

At  297  nm 

At  307  nm 

At  345  nm 

Intensity  ratio  R1  =  I 297/^345 

Intensity  ratio  R2  =  I307N345 

Normal 

0.36  ±  0.05 

0.57  ±  0.05 

0.99  ±  0.00 

0.36  ±  0.05 

0.57  ±  0.05 

Cancer  oua 

0.98  ±  0.01 

0.73  ±  0.03 

0.79  ±  0.05 

1.24  ±  0.09 

0.92  ±  0.05 

P  value 

<0.005 

<0.005 

<0.005 

<0.005 

<0.005 
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Fig.  4.  Scatter  plot  of  intensity  ratio  of  J307/J345  for  normal  (•) 
and  cancerous  (■)  prostate  tissues. 


wavelengths  of  the  spectral  peaks  observed  for  pros¬ 
tate  tissues  were  not  exactly  the  same  as  those  of 
the  standard  fluorophores.  This  observed  spectral 
shift  of  the  peak  positions  in  the  spectra  of  prostate 
tissues  may  be  due  to  the  different  microenviron¬ 
ments  of  these  endogenous  fluorophores  in  tissues. 
In  particular,  the  band  position  of  standard  fluoro- 
phore  of  NADH  at  397  nm  does  not  exactly  match 
our  prostate  tissue  band  position,  which  has  been 
observed  around  440  nm.  Earlier  reported  studies 
have  shown  that  exact  determination  of  NADH  emis¬ 
sion  is  difficult  because  this  compound  significantly 
changes  the  spectral  shape  and  intensity  after  under¬ 
going  oxidation  after  resection  [13] .  In  addition,  the 
fluorescence  data  of  NADH  compound  obtained  from 
ex  vivo  tissue  may  not  accurately  represent  the  in 
vivo  or  living  situation  [21].  Thus  oxidation  of  tissue 
fluorophores  will  influence  the  fluorescence  proper¬ 
ties  of  the  tissues,  which  leads  to  change  in  spectral 
peak  positions  and  the  line  width  and  fluorescence 
intensities  [22,23]. 

Distinct  salient  differences  in  the  spectral  signa¬ 
tures  of  tryptophan,  collagen  and/or  elastin,  NADH, 
and  FAD  were  observed  from  normal  and  different 
pathological  conditions  of  prostate  tissues,  as  shown 
in  Fig.  1.  The  normal  tissue  exhibits  a  small  peak 
around  307  nm,  most  likely  due  to  tryptophan,  an 
aromatic  amino  acid.  This  peak  was  significantly 
higher,  and  it  shifted  to  297  nm  for  cancerous  tissues. 
Table  1  also  shows  that  the  intensity  of  the  trypto¬ 
phan  emission  band  was  higher  in  the  cancerous 
tissues  (0.98  ±  0.01  at  297  nm  and  0.73  ±  0.03  at 
307  nm)  than  in  normal  tissues  (0.36  ±  0.05  at 
297  nm  and  0.57  ±  0.05  and  307  nm).  Significantly, 
the  tryptophan  emission  band  was  lower  (0.57  ±  0.05 
at  307  nm)  for  normal  tissues,  and  it  was  noticed  that 
the  same  band  was  blueshifted  with  increased  emis¬ 
sion  intensity  (0.98  ±  0.01  at  307  nm)  for  cancerous 
tissues.  The  striking  observation  is  the  increase  in 
the  tryptophan  fluorescence  signals  in  cancerous  tis¬ 
sues  relative  to  its  normal  counterpart.  This  increase 
is  probably  due  to  hyperactivity  or  epithelium  gland¬ 
ular  proliferation  resulting  in  higher  transient 


concentration  of  proteins  in  cells  and/or  increased 
thickness  of  the  epithelium  [9,24].  The  observed  in¬ 
crease  in  tryptophan  fluorescence  correlates  with  the 
earlier  reported  data  on  cervical  epithelial  cancer 
cells,  non-melanoma,  and  bladder  cancerous  tissues 
[4,5,25].  In  addition  to  the  increase  in  tryptophan 
fluorescence  intensity,  the  peak  emission  wavelength 
was  blueshifted  to  297  nm  for  cancerous  tissues 
in  comparison  with  the  normal  tissues.  This  observed 
blueshift  may  be  attributed  to  changes  in  the  local 
microenvironment  and  conformational  changes  of 
the  protein-bound  tryptophan,  relative  to  the  free 
amino  acid  [26].  Similar  spectral  blueshifts  have 
been  noted  in  human  laryngeal  epithelial  cells  [27], 
carcinogen-transformed  human  bronchial  epithelial 
cells  [28],  and  SS  spectra  of  breast  tissues  [19]. 
Further,  it  is  suggested  that  sensitivity  of  tryptophan 
fluorescence  and  its  peak  emission  due  to  polarity  and 
mobility  of  environment  makes  tryptophan  fluores¬ 
cence  useful  for  monitoring  structural  changes  in  a 
protein  [29].  Based  on  the  results,  it  is  concluded  that 
changes  in  intrinsic  tryptophan  fluorescence  and  its 
peak  emission  wavelength  from  cancerous  prostate 
tissues  can  be  used  as  one  of  the  diagnostic  criteria 
for  classifying  early  neoplastic  changes. 

The  structural  proteins  collagen  and/or  elastin 
have  been  noted  in  the  S3  spectrum.  As  we  know, 
collagen  and/or  elastin  is  a  major  component  of  extra¬ 
cellular  matrix  and  dominates  the  stoma  fluores¬ 
cence  associated  with  crosslinks.  From  Fig.  1,  the  pri¬ 
mary  prominent  band  was  observed  for  the  normal 
tissue  at  345  nm,  corresponding  to  structural  protein 
collagen  and/or  elastin,  with  a  small  shoulder  around 
375  nm,  and  secondary  peak  was  observed  at  405  nm 
due  to  elastin.  These  bands  are  drastically  decreased 
in  the  cases  of  cancerous  prostate  tissues.  Table  1 
also  shows  that  the  intensity  of  collagen  emission 
band  (0.99  ±  0.00  at  345  nm)  was  high  in  normal  tis¬ 
sues  and  low  (0.79  ±  0.05  at  345  nm)  in  cancerous  tis¬ 
sues.  In  normal  prostate  tissues,  the  collagen  and 
elastin  network  is  dense,  with  larger  numbers  of 
fibers  compared  to  non-uniform  of  size  and  shape  of 
the  epithelial  cells,  and  loss  and  disintegration  of  col¬ 
lagen  and  elastin  fibers  in  the  cancerous  prostate  tis¬ 
sue  was  expected  [24,30,31].  The  significant  decrease 
of  collagen  and/or  elastin  fluorescence  in  cancerous 
prostate  tissues  may  be  due  to  structural  changes  in 
the  epithelium,  and  sub-epithelial  collagen  network 
may  be  due  to  nuclear  enlargement  and  because 
epithelial  cells  are  crowded  with  increased  cell 
density,  loss  of  cellular  maturation,  and  overall  thick¬ 
ening  of  the  epithelium  [81032].  Thus,  increased 
thickness  of  the  epithelium  in  malignant  prostate 
glands  attenuates  both  excitation  light  and  emission 
fluorescence  light  from  the  vascular  stroma,  leading 
to  reduced  penetration  depth  of  the  excitation  light 
as  well  as  the  re-absorption  of  the  emitted  fluores¬ 
cence.  This  may  present  a  potential  criterion  for  pros¬ 
tate  cancer  detection.  This  criterion  agreement  with 
reported  studies  has  shown  that  collagen  fluores¬ 
cence  is  decreased  in  cancerous  prostate  tissues  [8], 
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The  characteristic  hemoglobin  absorption  features 
were  clearly  found  as  a  dip  at  420  nm  for  normal  and 
cancerous  prostate  tissues.  The  hemoglobin  absorp¬ 
tion  is  deeper  in  the  case  of  cancerous  tissues  than 
for  normal  tissues  due  to  possible  increase  of  hemo¬ 
globin  content  from  angiogenesis  [32].  Most  cancer¬ 
ous  tissue  is  known  to  exhibit  increased  vasculature 
due  to  angiogenic  developments  with  recurrence  in 
hypoxia,  and,  hence,  there  is  increased  blood  content 
in  superficial  tissue  layers,  which  leads  to  increased 
hemoglobin  absorption  when  compared  to  normal  tis¬ 
sues  [33],  This  hypothesis  is  also  consistent  with  the 
findings  of  earlier  reports  on  S3  of  cervical  [18]  and 
breast  tissues  [19]  and  other  studies  by  Demos  et  al. 
[14]  and  Volynskaya  et  al.  [34]. 

In  addition  to  this,  in  Fig.  1,  the  normal  tissue  ex¬ 
hibits  a  broad  peak  around  440  nm  and  has  a  small 
sharp  peak  at  510  nm  at  a  longer  wavelength.  How¬ 
ever,  cancerous  prostate  tissues  also  reveal  similar 
spectral  bands  but  with  a  less  fluorescence  intensity 
compared  with  the  normal  tissues.  The  observed 
spectral  bands  at  440  and  510  nm  may  be  attributed 
to  the  coenzyme  of  NADH  and  FAD.  The  relative  de¬ 
crease  of  NADH  and  FAD  fluorescence  in  malignant 
prostate  tissues  may  be  due  to  changes  in  tissue  mor¬ 
phology  (thickening  of  epithelium),  the  decrease  of 
fluorophores  (NADH  and  FAD)  quantum  yields, 
and/or  oxidation  of  tissue  fluorophores  [5,22], 

The  present  result  shows  that  the  overall  S3  signa¬ 
ture  provides  information  about  tissue  architecture 
(epithelial  thickness),  changes  in  absorption  due 
to  hemoglobin,  and  concentration  of  fluorescing  mo¬ 
lecules,  which  can  be  correlated  to  histological 
changes.  Therefore,  alteration  in  tissue  architecture 
that  inhibits  excitation  photons  from  reaching  the 
native  fluorophores  and/or  the  emission  photons 
from  the  fluorophores  from  escaping  from  the  tissue 
could  be  detected  by  the  S3.  Additionally,  as  the  he¬ 
moglobin  concentration  increased,  the  fluorescence 
intensity  decreased  and  the  fluorescence  profile  was 
altered  significantly.  Thus,  changes  in  the  concentra¬ 
tion  of  the  key  native  fluorophores  could  alter  the 
emitted  fluorescence  in  a  single  scan,  which  is  man¬ 
ifested  in  our  study. 

Many  studies  have  been  reported  that  the  fluores¬ 
cence  spectral  intensity  ratio  could  be  used  as  a  po¬ 
tential  tool  for  tissue  diagnosis  [5,18-20].  In  this 
paper,  a  simple  classification  method  was  adopted 
based  on  the  ratio  values  estimated  from  the  inten¬ 
sities  at  different  wavelengths  characterizing  the 
peaks  of  different  groups  under  study.  To  quantify 
the  observed  spectral  results  and  to  estimate  the  di¬ 
agnostic  potentiality  of  the  present  technique,  two 
intensity  ratios,  such  as  I297//345  f°r  normal  tissues 
and  /307/I345  for  cancerous  tissues,  were  selected. 
Each  emission  wavelength  used  in  these  ratios  repre¬ 
sents  a  specific  fingerprint  of  one  fluorophore.  The 
component  identified  as  /30 7  represents  the  peak 
emission  wavelength  of  tryptophan  for  normal  tis¬ 
sues,  whereas  I2 07  represents  the  peak  emission 
wavelength  of  tryptophan  for  cancerous  tissues,  and 


/345  represents  the  peak  emission  wavelength  of  col¬ 
lagen  and/or  elastin  from  normal  and  cancerous 
tissues.  Both  intensity  ratio  values  increase  signifi¬ 
cantly  for  cancerous  tissues  in  comparison  with  nor¬ 
mal  tissues.  From  these  ratio  parameters,  we  found 
variations  in  the  relative  distribution  of  tryptophan, 
collagen  and/or  elastin  corresponding  to  the  different 
histopathology  of  the  tissues. 

In  this  paper,  we  performed  the  classification  ana¬ 
lysis  for  two  ratio  parameters  C/297/^345  and 
/307//345)  in  order  to  find  out  the  specificity  and  sen¬ 
sitivity  of  the  present  technique  to  discriminate  be¬ 
tween  two  groups.  Figure  3  shows  the  scatter  plot  of 
the  spectral  intensity  ratio  of  /297/-Z345  from  sixteen 
samples  from  eight  patients  categorized  as  normal 
prostate  tissues  ( n  =  8)  and  cancerous  prostate  tis¬ 
sues  (n  =  8).  The  classification  specificity  and  sensi¬ 
tivity  in  discriminating  each  of  these  categories 
were  determined  on  the  basis  of  discrimination  cut¬ 
off  values,  given  in  the  scatter  plot  (Fig.  3).  For  ex¬ 
ample,  the  cutoff  line  discriminating  the  normal 
from  the  cancerous  was  drawn  at  0.80,  correspond¬ 
ing  to  the  mean  ratio  of  normal  (0.36)  and  the  mean 
of  the  cancerous  (1.24)  group.  For  the  /297/-Z345  ratio, 
by  selecting  a  cutoff  at  the  mean  (0.80)  of  the  normal 
and  cancerous  tissues  yield,  100%  specificity  and 
100%  sensitivity  was  obtained  for  discrimination 
of  the  normal  from  the  cancerous.  Similarly,  we  per¬ 
formed  other  ratio  parameter  /307  //3 45  for  the  same 
classification  analysis.  Figure  4  shows  the  scatter 
plot  of  the  spectral  intensity  /3o7/I345  ratio,  the  cut- 
offline  drawn  at  0.74  discriminates  the  normal  from 
cancerous,  with  the  specificity  and  sensitivity  of 
100%  and  100%,  respectively. 

Statistical  analysis  showed  that  both  intensity  ra¬ 
tios  yielded  100%  diagnostic  sensitivity  and  specifi¬ 
city.  Since  the  increased  value  of  the  fluorescence 
ratios  of/297/-f345  and /307//345  for  cancerous  prostate 
tissues  compared  to  normal  tissue  reflects  an  in¬ 
crease  in  the  tryptophan  fluorescence  contribution 
and  a  decrease  in  the  collagen  and/or  elastin  fluores¬ 
cence  contribution.  Among  the  all  the  fluorophores, 
the  finding  of  increased  band  of  tryptophan  fluores¬ 
cence  and  decreased  fluorescence  of  collagen  and/or 
elastin  band  from  prostate  cancer  could  be  particu¬ 
larly  useful  for  detecting  prostate  cancer. 

The  major  advantage  of  the  S3  technique  is  that  all 
the  fluorophores  in  complex  structures,  such  as  tis¬ 
sue,  are  excited  under  optimal  conditions,  such  as 
wavelength  of  excitation,  slit  widths,  scan  speed, 
and  scan  range,  which  is  practically  impossible  in 
the  case  of  conventional  fluorescence  measurements. 
In  the  case  of  conventional  fluorescence  measure¬ 
ments,  optimal  conditions  are  chosen  based  on  the 
fluorophores  of  interest.  Though,  EEM  offers  great 
resolution,  it  takes  some  time  to  carry  out,  even  with 
multichannel  detectors.  Also,  the  large  amount  of 
data  of  EEM  could  significantly  slow  down  the  data 
processing  procedures  and  clearly  is  impractical  in  a 
clinical  setting.  On  the  other  hand,  S3  represents  the 
diagonal  scan  over  the  entire  EEM,  thereby  reducing 
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the  time  of  data  acquisition  without  any  alterations 
in  the  details  over  the  entire  spectral  range  of  the  dif¬ 
ferent  fluorophores  present  in  the  different  tissues. 
The  S3  technique  reduces  the  time  of  data  acquisi¬ 
tion  and  can  be  identified  in  a  specific  spectral  range 
with  narrower  spectral  band  attributed  by  each 
fluorophore  than  the  conventional  fluorescence  spec¬ 
trum.  S3  improves  the  diagnostic  capability  of  optical 
spectroscopy  for  prostate  cancer  diagnosis  because 
all  the  key  fluorophores,  such  as  tryptophan,  collagen 
and/or  elastin,  hemoglobin  absorption,  NADH,  and 
FAD,  can  be  obtained  in  a  single  scan  [16] ,  and  they 
are  reported  as  cancer  markers  [8-11] .  On  the  basis 
of  above  factors,  S3  may  also  be  considered  as  a  tech¬ 
nique  alternative  or/and  complementary  to  the  exist¬ 
ing  conventional  methods  for  tissue  diagnosis. 

5.  Conclusion 

This  paper  suggests  that  S3  holds  the  potential  for 
differentiating  between  diseased  and  normal  pros¬ 
tate  tissues  in  vitro.  Results  of  the  current  study  de¬ 
monstrate  that  spectral  changes  due  to  the  changes 
of  contents  of  tryptophan,  collagen  and/or  elastin, 
NADH,  and  FAD  have  good  diagnostic  potential. 
Two  ratio  parameters  (i?i  =  and  R2  =  I  urn  / 

1 345)  were  empirically  selected  and  evaluated  by  an 
unpaired  Student’s  t-test,  which  showed  useful  diag¬ 
nostic  information  with  excellent  discrimination  of 
mean  ratio  values  between  normal  and  cancerous 
prostate  groups.  To  improve  on  the  diagnostic  cap¬ 
ability  of  S3,  a  larger  number  of  samples  of  normal 
and  different  pathological  prostate  tissues  would 
need  to  be  validated  in  blinded  manner.  Such  inves¬ 
tigations  with  a  large  group  of  tissue  biopsies  will  be 
useful  for  the  development  of  a  statistical  database 
and  a  user-friendly  diagnostic  algorithm  that  could 
facilitate  a  real-time  in  vivo  clinical  diagnosis  of  pros¬ 
tate  cancer.  Currently,  we  are  working  on  the  above 
aspects  to  prove  S3  as  a  useful  and  novel  technique  in 
the  diagnostic  clinical  oncology  for  different  tissue 
types. 
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A  paper  entitled  “Near  Infrared  Photonic  Finger  Imager  for  Prostate  Cancer  Screening”,  by  Y. 
Pu,  W.  B.  Wang,  M.  Xu,  G.  C.  Tang,  Y.  Budansky,  M.  Sharanov,  S.  Achilefu,  J.  A.  Eastham  and 
R.  R.  Alfano,  published  in  Technol.  Cancer  Res.  Treat.,  ( TCRT ),  J_0,  507-517  (2011). 
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A  portable  rectal  near  infrared  (NIR)  scanning  polarization  imaging  unit  with  an  optical 
fiber-based  rectal  probe,  designated  as  a  Photonic  Finger  (PF),  was  designed,  developed, 
built  and  tested.  PF  was  used  to  image  and  locate  the  three  dimensional  (3D)  positions  of 
abnormal  prostate  tissue  embedded  inside  normal  prostate  tissue.  An  inverse  image  recon¬ 
struction  algorithm,  namely  Optical  Tomography  using  Independent  Component  Analysis 
(OPTICA)  was  developed  to  unmix  the  signal  from  targets  (cancerous  tissue)  embedded  in 
a  turbid  media  (normal  tissue)  in  the  backscattering  imaging  geometry.  The  Photonic  Finger 
combined  with  OPTICA  was  ex  vivo  tested  to  characterize  different  target(s)  inside  different 
tissue  medium,  including  cancerous  prostate  tissue  embedded  inside  large  pieces  of  normal 
tissue.  This  new  developed  instrument,  Photonic  Finger,  may  provide  an  alternative  imaging 
technique,  which  is  accurate,  of  high  spatial  resolution  and  non-or-less  invasive  for  prostate 
cancers  screening. 

Key  words:  Photonic  finger  (PF);  Backscattering;  Scanning  optical  polarization  imaging; 
Optical  fiber-based  rectal  probe;  Three-dimensional  localization;  Pluman  prostate  tissue; 
Receptor-target  contrast  agents;  Independent  component  analysis  (ICA);  Photon  propagation 
model. 

Introduction 

In  year  2010,  217,730  new  cases  of  prostate  cancer  were  diagnosed,  and  approxi¬ 
mately  32,050  men  died  from  prostate  cancer  in  the  U.S.A.  (1).  The  common 
screening  tests  for  prostate  cancer  diagnosis  are  digital  rectal  examination 
(DRE),  measurement  of  the  serum  tumor  marker  namely  prostate  specific  antigen 
(PSA),  and  the  transrectal  ultrasound  (TRUS)  imaging  (2).  A  value  of  PSA  over 
4.0ng/ml  is  the  commonly  used  threshold  for  further  diagnostic  evaluation  (3). 
Although  PSA  test  appears  to  have  acceptable  sensitivity  for  late  stage  cancers 
(4)  and  disease  with  histopathologic  features  associated  with  tumor  progression 
of  a  large  volume,  poorly  differentiated  cells,  and  extracapsular  penetration  (5), 
its  accuracy  is  limited  as  low  as  28%-35%  (6).  During  the  DRE,  a  doctor  inserts 
a  lubricated,  gloved  finger  into  the  patient’s  rectum  to  feel  for  lumps,  enlarge¬ 
ments,  or  hard  areas  of  prostate  that  might  indicate  prostate  cancer.  DRE  has 
a  reported  sensitivity  of  18%-22%  (7,  8).  TRUS  is  no  longer  considered  as  a 
first-line  screening  test  for  prostate  cancer  (9)  because  of  its  poor  spatial  resolu¬ 
tion  and  contrast,  but  it  does  play  a  role  in  mapping  the  locations  of  the  biopsy 
sampling  (9).  The  confirmation  of  prostate  cancer  finally  needs  a  needle  biopsy 
of  the  prostate.  In  the  biopsy,  a  number  [  12  to  18]  of  cores  of  prostate  tissue  are 


Abbreviations:  Photonic  Finger  (PF);  Independent  Component  Analysis  (ICA);  Optical  Tomography 
using  Independent  Component  Analysis  (OPTICA). 
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randomly  taken  from  whole  region  of  the  prostate  using  a 
thin  needle  with  the  help  of  TRUS  to  map  the  locations  of  the 
sampling  (10). 

The  early  detection  and  treatment  of  prostate  cancers  can 
reduce  mortality  (11).  Conventional  oncology  imaging  meth¬ 
ods  for  prostate  cancer  diagnosis  still  depend  on  bulk  physical 
properties  of  cancer  tissue  and  are  not  effective  for  early-stage 
primary  tumors  (12).  Since  PSA  and  DRE  have  limited  accu¬ 
racy,  TRUS  has  poor  contrast  between  normal  and  abnormal 
tissue  regions,  and  needle  biopsy  is  invasive  and  may  cause 
damages  of  prostate,  it  is  highly  desirable  to  develop  a  better 
method  which  is  accurate,  of  higher  spatial  resolution  and 
non-or-less  invasive  for  prostate  cancer  screening. 

Optical  imaging  technique  using  near  infrared  (NIR)  light 
provides  an  attractive  noninvasive  approach  for  screening 
human  diseases.  The  “tissue  optical  window”  in  the  NIR 
range,  which  corresponds  to  lower  absorption  from  major 
tissue  chromophores  such  as  water,  oxygenated  and  deoxygen- 
ated  hemoglobin,  allows  light  to  penetrate  deep  into  the  tissue, 
up  to  several  centimeters  (12).  The  other  main  advantages  of 
NIR  optical  approaches  are  its  inherent  low-cost,  the  ability  to 
monitor  multiple  independent  optical  reporters  simultaneously 
in  vivo  based  upon  wavelength,  the  absence  of  radioactive 
intermediates,  and  the  relative  simplicity  of  the  imaging  hard¬ 
ware  as  compared  with  Magnetic  resonance  imaging  (MRI) 
and  Positron  emission  tomography  (PET)  equipment.  These 
advantages  make  optical  imaging  unmatched  by  any  other  in 
vivo  imaging  techniques.  Over  the  past  decade,  Indocyanine 
Green  (ICG,  also  called  cardio-green),  a  clinically  approved 
NIR  dye  by  FDA,  has  been  investigated  as  a  contrast  agent 
for  optical  detection  of  tumors.  However,  ICG  is  not  designed 
to  specifically  target  cancer  cells.  The  investigations  of  recep¬ 
tor  expression  in  normal  and  cancer  tissues  suggest  that  small 
peptide-dye  conjugates  can  be  used  to  target  over-expressed 
receptors  on  tumors  to  enhance  specificity  (13).  Biological 
studies  have  indicated  that  somatostatin  receptors  (SSTR)  are 
over-expressed  in  human  prostate  tumor  (14).  The  previous 
investigation  showed  that  a  small  ICG-derivative  dye-peptide, 
namely  Cypate-Octreote  Peptide  Analogue  Conjugate  (Cytate) 
could  be  used  for  effectively  targeting  somatostatin  receptor- 
rich  tumor  in  the  animal  model  because  of  the  high  affinity  of 
Cytate  for  the  somatostatin  receptors  (13).  This  motivates  us 
to  apply  Cytate  in  optical  scanning  imaging  in  human  prostate 
tissue  for  cancer  detection.  The  major  disadvantage  of  optical 
imaging  approaches  is  that  strong  scattering  by  biological  tis¬ 
sue  causes  most  photons  diffused.  The  very  few  percentage  of 
ballistic  and  snake  photons  make  direct  imaging  practical  only 
in  surface  and  subsurface  layers  of  tissue  (15,  16).  Scientists 
have  to  explore  optical  tomography  methods  and/or  inverse 
image  reconstruction  approaches  to  locate  the  three  dimension 
(3D)  positions  of  abnormal  tissue  or  recover  the  3D  spatial  dis¬ 
tribution  information  of  optical  parameters  of  the  tissue  (16). 


In  this  research,  a  portable  rectal  NIR  scanning  polarization 
imaging  unit  with  an  optical  fiber-based  rectal  probe,  named 
Photonic  Finger  (PF),  was  developed.  The  transrectal  opti¬ 
cal  imaging  approach  was  used  to  locate  the  3D  positions 
of  abnormal  prostate  sites  hidden  in  normal  prostate  tissue 
based  on  differences  of  optical  parameters  between  cancer¬ 
ous  and  normal  prostate  tissues.  The  difference  can  also  be 
enhanced  using  an  extrinsic  chromophore  or  fluorophore  such 
as  a  receptor-targeted  contrast  agent.  The  rectal  NIR  scan¬ 
ning  polarization  imaging  and  3D  inverse  location  technique 
has  potential  to  address  the  critical  issues  of  the  conventional 
screening  methods  for  prostate  cancer  detection.  The  scanning 
polarization  imaging  acquires  2D  images  by  sequentially  scan¬ 
ning  a  polarized  illuminating  light  beam  at  different  areas  of  a 
prostate  gland  through  rectum,  and  recording  the  distribution 
of  light  intensity  backscattered  from  the  prostate  using  a  CCD 
camera.  An  Independent  Component  Analysis  (ICA)-based 
inverse  3D  location  reconstruction  algorithm  was  improved 
specifically  for  the  application  of  backscattering  configuration 
and  used  to  locate  the  3D  positions  of  foreign  inhomogene¬ 
ities  from  the  recorded  array  of  the  2D  images.  This  research 
provides  a  noninvasive  optical  imaging  technique  for  detect¬ 
ing  and  3D  locating  cancerous  sites  in  prostate.  Therefore, 
PF  may  introduce  a  new  criteria/indicator  for  prostate  can¬ 
cer  screening  in  addition  to  the  conventional  examinations  to 
enhance  the  accuracy  of  prostate  cancer  detection. 

Experimental  Setup  and  Methods 

Design  and  Construction  of  Photonic  Finger 

The  Photonic  Finger  is  a  portable  rectal  NIR  scanning  polar¬ 
ization  imaging  unit  with  an  optical  fiber-based  rectal  probe. 
The  unit  was  designed  to  be  capable  of  recording  sets  of  2D 
images  of  the  prostate  by  scanning  the  illuminating  beam 
on  the  prostate  through  rectum  walls.  The  major  optical  and 
electronic  components  of  the  unit  include  laser  diodes  and 
their  power  supplies,  miniature  scanning  Galvanometric  mir¬ 
rors  with  their  electronic  control  boards  (servo  driver  circuit 
boards)  and  the  Lab  VIEW  control  software,  and  illumination 
and  imaging  coherent  optical  fiber-bundles. 

A  photograph  of  the  portable  rectal  NIR  scanning  polariza¬ 
tion  imaging  unit  is  shown  in  Figure  1(A).  Three  diode  lasers 
emitting  at  635 nm,  750nm,  and  980nm,  respectively,  are 
alternatively  used  as  light  sources.  These  wavelengths  were 
selected  to  probe  the  native  molecules  such  as  Hb,  Hb02 
and  H20  in  tissues.  The  output  beam  from  a  laser  diode  is 
directed  to  the  scanning  galvanometric  mirror  system  (Thor- 
Lab  GVSM002  with  Dual  Axis  Galvo  Mirrors)  after  passing 
through  two  pinholes  and  a  polarizer  (PJ.  The  beam  can  be 
scanned  by  two  miniature  galvanometric  mirrors  in  the x-  and 
v-directions,  respectively.  The  beam  output  from  the  galva¬ 
nometric  mirrors  is  focused  using  a  microscopy  objective 
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Figure  1:  (A)  A  photograph  of  the  Photonic  Finger:  a  portable  rectal  NIR  scanning  polarization  imaging  unit.  The  power  supplies  and  computer  system  are 

not  shown  in  the  picture.  (B)  The  schematic  diagram  of  the  portable  rectal  NIR  scanning  polarization  imaging  unit,  where  L-laser,  M-mirror,  A-aperture, 
P-polarizer,  GV-galvanometric  scanning  mirror,  BPF-band  pass  filter,  MOL-microscopy  objective  lens. 


lens  on  a  coherent  fiber-bundle  (Mytiad  Fiber  Imaging, 
20-0826  Fiberscope  Assy)  used  for  illumination.  The  posi¬ 
tion  of  the  microscopy  objective  lens  can  be  adjusted  in  the 
x-,  y-  and  z-  three  directions  for  the  beam-fiber  coupling. 
The  output  beam  from  the  illumination  fiber  is  directed 
to  a  small  reflection  prism  located  inside  the  rectal  probe. 
The  beam  reflected  from  the  prism  is  used  to  illuminate  a 
prostate  sample.  This  illumination  beam  with  a  diameter 
of  ~1  mm  can  be  scanned  in  the  x-  and  y-directions  on  the 
prostate  sample. 


The  light  backscattered  (or  emitted)  from  a  prostate  sample 
is  first  passing  through  the  rectal  wall  and  reflected  from 
the  prism  inside  the  probe  head.  The  diameter  of  the  probe 
is  ~2cm  and  the  length  is  ~12cm.  The  beam  is  then  col¬ 
lected  by  a  lens  into  another  coherent  fiber-bundle  (Mytiad 
Fiber  Imaging,  20-0826  Fiberscope  Assy)  used  for  imag¬ 
ing.  The  diameter  of  a  single  fiber  in  the  bundle  is  ~3.2 
pm  each  with  numerical  aperture  N.A.  =  0.4.  The  image 
information  formed  in  the  optical  fiber  bundle  was  sent  to  a 
CCD  camera  through  the  coherent  imaging  fiber  bundle  and 
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coupling  lens  for  recording  2D  images  of  the  prostate  sample. 
The  coupling  loss  of  the  system  is  -10%  and  the  transverse 
resolution  is  100  Jim  associated  with  Air  Force  resolution 
target  bar  chart  (AFBC)  at  group  3.  A  polarizer  (P2)  and 
a  band  pass  filter  are  placed  in  front  of  the  CCD  camera 
to  record  the  2D  images  at  different  polarization  configura¬ 
tion  and  wavelengths.  A  cross-polarization  image  for  each 
scanned  illumination  position  of  the  laser  beam  is  recorded 
when  the  polarization  direction  of  P2  is  perpendicular  to 
that  of  P!  to  suppress  the  contribution  of  light  scattered  (or 
emitted)  from  the  surface  and  sub-surfaces  to  the  images  of 
the  prostate  sample.  For  the  backscattering  light  imaging,  a 
narrow  band  filter  corresponding  to  the  illumination  wave¬ 
length  is  used  in  front  of  the  CCD  camera  to  ensure  that  the 
recorded  images  are  formed  only  by  the  light  backscattered 
from  the  prostate  sample.  For  the  tissue  emission  and/or 
contrast  agent  emission  light  imaging,  a  long  pass  filter  is 
used  in  front  of  the  CCD  camera  to  ensure  that  the  recorded 
images  are  formed  only  by  the  light  emitted  from  the  pros¬ 
tate  sample.  When  the  illumination  light  beam  is  scanned  in 
the  x-y  plane  of  the  sample  with  n  X  n  points,  an  array  of  n1 
2D  images  will  be  recorded. 

The  key  part  of  the  scanning  imaging  unit  is  the  scan  of  the 
illumination  beam  on  the  surface  of  the  prostate  sample  with 
adjustable  scanning  parameters  such  as  scanning  area,  step, 
speed  and  the  number  of  scanning  points.  Figure  1(B)  sche¬ 
matically  shows  the  layout  of  the  major  parts  and  the  con¬ 
trol  boards  for  the  PF  unit.  This  scanning  system  consists  of 
(1)  the  Dual  Axis  Galvo  mirrors  (ThorLab  GVS  002),  (2)  the 
two  servo  driver  circuit  boards  and  their  power  supply  (Thor¬ 
Lab,  GPS011),  (3)  the  drive  unit  for  sending  voltage  output 
to  servo  circuit  boards  (National  Instruments,  NI  DAQmx 
USB-9263  USB  DAQ  -  data  acquisition),  and  (4)  a  PC  with 
an  installed  Lab  VIEW  software  (Lab  VIEW  Signal  Express 
2009)  to  power  and  send  a  command  to  DAQ  USB  9263  to 
generate  desired  output  voltage  through  a  USB  connection.  In 
the  scanning  system,  the  drive  unit  NI-9263  powered  by  USB 
interface  of  the  PC,  is  used  to  generate  an  analog  voltage, 
which  can  be  varied  from  — 10Vto+10V  using  the  Lab- 
VIEW  software.  This  analog  voltage  is  sent  to  the  two  servo 
driver  boards  to  drive  the  rotations  of  the  Galvo  Mirrors.  For 
example,  IV  input  to  the  servo  board  can  make  the  mirror 
rotate  1°.  One  servo  board  controls  the  X-axis,  and  another 
is  for  Y-axis.  Both  servo  driver  boards  are  powered  by  the 
power  supply  of  GPS011.  The  beam  scanning,  the  image 
acquiring  and  recording,  and  the  synchronization  of  scan¬ 
ning  and  imaging  are  controlled  by  the  Graphical  User  Inter¬ 
face  (GUI)  software  developed  using  LabVIEW.  Figure  2 
shows  a  screenshot  of  the  GUI  after  completing  an  n  X  n 
scan.  All  parameters  of  the  scanning  imaging  unit  (the  posi¬ 
tion  of  the  original,  scanning  steps,  step  size,  exposure  time, 
and  waiting  time  between  two  adjacent  imaging  acquiring), 
can  be  adjusted  through  the  GUI. 


Test  Samples 

The  Photonic  Finger  scanning  imaging  unit  was  tested  for 
different  types  of  tissue  samples.  The  first  sample  used  for 
scanning  imaging  was  a  black  rubber  absorber  with  the  size 
of  -2.7  mm  X -3.0 mm  X -1.5  mm  embedded  in  chicken 
breast  tissue  at  different  depth.  The  fresh  chicken  breast  tis¬ 
sues  were  purchased  at  the  local  super  market.  The  black 
rubber  was  covered  by  the  chicken  breast  tissues  at  two  dif¬ 
ferent  depths  of  -5.1  mm  and  -6.6 mm  with  a  lateral  dimen¬ 
sion  38  mm  X  29  mm  for  scanning  imaging  measurements. 
The  thickness  of  the  tissue  behind  the  object  is  -8  mm.  A 
transport  length  lt=  1.1mm  and  an  absorption  coefficient 
pa  =  0.007  mm  were  taken  for  chicken  breast  tissue  for 
imaging  analysis  (17). 

The  second  study  was  performed  to  distinguish  in  vitro 
cancerous  prostate  tissue  from  surrounding  normal  prostate 
tissue.  Human  prostate  tissues  were  obtained  from  the  Co¬ 
operation  Human  Tissue  Network  (CHTN)  and  the  National 
Disease  Research  Interchange  (NDRI)  under  the  approval  of 
the  Institutional  Review  Board  (IRB)  at  CCNY.  The  cancer¬ 
ous  and  normal  prostate  samples  were  diagnosed  by  patholo¬ 
gist  before  the  optical  imaging  experiments.  Samples  were 
neither  chemically  treated  nor  were  frozen  prior  to  the  experi¬ 
ments.  The  time  elapsed  between  tissue  resection  and  taking 
the  scanning  imaging  measurements  may  vary  for  different 
sample  sources.  The  longest  elapsed  time  is  about  30  hours. 
The  sample  consists  of  a  small  piece  of  cancerous  prostate  tis¬ 
sue  (4  mm  X  4  mm  X  1.5  mm)  embedded  inside  a  large  piece 
of  normal  prostate  tissue  at  the  depth  of  z  =  3.0mm  from 
the  front  surface.  The  thickness  of  the  whole  tissue  sample 
is  10  mm.  For  imaging  analysis,  the  absorption  and  reduced 
scattering  coefficients  for  normal  prostate  tissue  (the  host 
medium)  were  taken  as  pa  ~  0.026mm-1  and p'  ~  0.53  mm-1. 
The  cancerous  prostate  tissue  absorbs  much  less  light  with 
pa  ~  0.0025  mm-1  and  scatters  less  with  p'  ~  0.44mm-1  at 
the  probing  wavelength  of  635  nm  estimated  from  our  earlier 
spectroscopic  investigations  of  ex  vivo  prostate  tissues  (18). 
The  piece  of  cancerous  prostate  tissue  embedded  behaves 
predominantly  as  an  absorption  inhomogeneity  in  homoge¬ 
neous  normal  prostate  tissue.  Here,  we  simply  assume  that 
the  optical  parameters  of  normal  tissue  are  macroscopically 
homogeneous  throughout  the  tissue  volume.  Although  this  is 
not  a  rigorous  description  of  tissue,  for  many  cases  of  inter¬ 
est,  it  will  be  sufficient.  OPTICA  is  applied  here  to  extract 
different  independent  components  by  treating  cancerous  tis¬ 
sue  as  inhomogeneity  embedded  inside  a  “homogeneous” 
normal  tissues. 

The  third  study  was  designed  to  characterize  cancerous  prostate 
tissue  after  enhancing  imaging  contrast  by  extrinsic  fluores¬ 
cent  marker,  a  receptor-targeted  contrast  agent,  namely  Cytate. 
Cy  tate  was  synthesized  by  Achilefu’ s  group  at  the  W ashington 
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Figure  2:  A  screenshot  of  the  Graphical  User  Interface  (GUI)  after  completing  annXn  scan.  The  GUI  was  developed  using  Lab  VIEW  for  operation  of  the 
Photonic  Finger.  All  parameters  of  the  scanning  imaging  (the  position  of  the  original,  scanning  steps,  step  size,  exposure  time,  and  waiting  time  between  two 
adjacent  imaging  acquiring,  e.g .)  can  be  adjusted  through  the  GUI. 


University  School  of  Medicine.  This  contrast  agent  is  mainly 
composed  of  ICG  and  the  somatostatin  receptor  ligand, 
which  delivers  ICG  to  the  corresponding  somatostatin  recep¬ 
tors  over-expressed  in  the  tumor  (13,  14).  The  prostate  tis¬ 
sue  samples  for  the  scanning  imaging  study  were  prepared 
following  this  protocol:  (1)  samples  (cancerous  and  normal 
prostate  tissues)  were  cut  into  very  tiny  size,  which  is  less 
than  ~lmm  (in  length,  width  and  thickness)  pieces.  Each 
pair  of  cancerous  tissue  and  the  corresponding  normal  tissue 
used  as  a  control  sample  was  obtained  from  same  patient;  (2) 
both  small  pieces  of  prostate  tissue  samples  were  soaked  in  a 
same  Cytate  solution  with  a  concentration  of  ~3.2  x  1 0  6  M 
for  -lOminutes;  (3)  Cytate-stained  tissue  samples  were  put 
into  sodium  phosphate  buffer  (Sigma-Aldrich)  to  wash  off 
and  consequently  reduce  the  amount  of  un-bound  Cytate  (19); 
and  (4)  the  Cytate-stained  small  pieces  of  cancerous  and 
normal  prostate  tissues  were  covered  by  a  large  piece  of  nor¬ 
mal  prostate  tissue  at  two  different  depths  of  ~2.5  mm  and 
~3.7  mm  for  scanning  imaging  measurements. 

The  orthogonal  co-ordinate  system  to  locate  3D  positions  of 
the  target  in  the  scanning  polarized  image  measurements  is 
schematically  shown  in  Figure  3.  The  samples  to  be  tested 
were  embedded  inside  large  slices  of  host  tissue  at  differ¬ 
ent  depths.  The  illuminating  laser  beam  was  scanned  along 
the  x-  and  y-directions  on  the  front  surface  of  the  sample. 
The  2D  images  formed  by  the  light  backscattered  from  the 


sample  in  the  normal  direction  are  recorded.  In  the  first  and 
second  studies,  only  one  target  (black  rubber  or  cancerous 
prostate  tissue)  is  involved  in  each  case  and  the  sample  was 
illuminated  by  a  collimated  laser  beam  with  X  =  635  nm. 
The  elastic  backscattering  from  the  samples  was  collected 
through  a  narrow  band  filter  (NF).  In  the  third  study  for 
contrast  agents’  fluorescence,  the  Cytate-stained  tiny  can¬ 
cerous  and  normal  prostate  tissue  samples  were  embedded 
inside  large  pieces  of  normal  prostate  tissue.  The  sample 
was  illuminated  by  a  laser  beam  with  X  =  635  nm  in  the 
direction  close  to  the  normal  to  the  surface.  A  narrow  band¬ 
pass  filter  at  830  nm  was  placed  in  front  of  a  CCD  camera 
to  record  emission  images  based  on  the  emission  peak  of 
Cytate  at  -837  nm  (20). 

Algorithm  of  OPTICA 

The  theory  of  Optical  Tomography  using  Independent  Com¬ 
ponent  Analysis  (OPTICA)  method  was  described  in  else¬ 
where  (21).  ICA,  the  core  of  OPTICA,  is  a  solution  to  the  blind 
source  separation  problems,  which  is  a  class  of  the  problems 
that  no  precise  knowledge  is  available  on  neither  the  mixing 
channels  nor  the  sources.  Typically  the  observations  are  the 
output  of  a  set  of  sensors,  where  each  sensor  receives  a  set  of 
mixed  source  signals  (21).  The  algorithm  of  OPTICA  is  based 
on  that  ICA  of  the  perturbations  in  the  spatial  intensity  dis¬ 
tributions  provides  the  corresponding  independent  intensity 
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Figure  3:  The  schematical  diagram  that  displays  our  calculated  position  of 
objects  in  an  orthogonal  co-ordinate  system  to  locate  3D  positions  of  the 
target  using  the  scanning  polarized  image  unit. 

distributions  (IID)  on  the  detector  planes  (21).  The  measured 
array  image  is  given  by  (16): 

I  (P,b  Ps )  =  foGo  (rd,  rs)  -  I0Az}G0(rd\  p.  z)djujp',  z') 

G0  ip',  z\  rs)  d1 2 3 4 5 6p  I '  ] 

where  pd  covers  the  whole  2D  array,  ps  is  the  scanning  posi¬ 
tions  of  incident  light  source  with  the  intensity  of  I0  at  the 
Z  =  0  plane,  G0  denotes  the  exact  Green’s  function  for  light 
propagation  in  the  host  medium.  It  is  assumed  that  one  tiny 
object  (with  absorptive  perturbation  dpa  of  volume  AV)  is 
located  at  r  =  ( p ,  z')  with  the  extension  Az  «  1  along  the 
axial  direction  and  that  dpa  is  constant  within  AV.  The  speed 
of  light  is  set  to  be  unity.  Scattering  targets  can  be  treated  in 
a  similar  fashion. 

The  background  image  could  be  approximately  generated 
using  an  “average”  image  of  all  acquired  images  by  shifting 
the  scanning  position  for  each  image  ps  to  the  origin  0,  which 
is  given  by  (16): 

T(pd,0)  =  IBG0(rd,0)-h(pd,0)  [2] 

where  the  term  h  stands  for  the  error  term.  This  “averaged” 
background  image  is  called  “dirty”  background  image  since 
the  object(s)  with  dps  and  Sps  is  included  in  the  averaging 
calculation.  To  generate  a  “clean”  background  image,  the 


shifted  images  which  are  minimally  perturbed  by  the  embed¬ 
ded  objects  were  selected  and  averaged  (16): 

Tc=T(pd,0)+h(pd,0)  [3] 


where  the  error  term  is  chosen  by: 

%-/?’0)VlA fU-ZJ.O)  [4] 

™  B  p,eB 

where  B  denotes  the  perimeter  of  the  scanning  grid  which 
contains  NB  scanning  positions.  After  h  and  the  clean  image 
of  the  host  medium  have  been  obtained,  the  difference  images 
are  calculated  more  accurately  (16)  since  the  cleaning  imag¬ 
ing  suppresses  the  contributions  from  embedded  objects 
(cancerous  tissues)  to  simulate  a  more  “homogeneous”  back¬ 
ground. 

The  difference  images  were  then  generated  by  subtracting 
the  “clean”  background  image  from  the  recorded  2D  images. 
Based  on  the  difference  images,  ICA  was  used  to  unmix  the 
signal  arising  from  individual  targets  and  the  independent 
intensity  distribution  due  to  the  target(s).  Each  target  is  asso¬ 
ciated  with  one  independent  component  (IC),  which  consists 
of  the  projection  of  the  Green’s  functions,  G0(rs,  rt)  and  G0 
( rd ,  r,),  on  the  source  or  detector  plane,  respectively  (21).  The 
position  of  the  target  is  obtained  by  numerically  marching 
the  target(s)  to  the  surface  until  matching  the  retrieved  inde¬ 
pendent  component,  incorporating  both  the  beam  profile  and 
the  surface  property  of  the  sample  (14).  The  3D  locations  of 
a  black  rubber  hidden  inside  chicken  breast  tissue  medium, 
and  3D  locations  of  a  small  piece  of  cancerous  prostate  tis¬ 
sue  embedded  in  the  host  normal  prostate  tissue  were  used 
to  verify  our  new  OPTICA  algorithm  using  the  recorded  sets 
of  2D  backscattering  images.  The  procedures  of  OPTICA 
for  obtaining  3D  locations  of  objects  in  scattering  media  are 
listed  in  Table  I. 


Table  I 

Steps  of  OPTICA  to  locate  3D  position  of  target(s)  inside  scattering 

medium  at  backscattering  configuration. 

1 .  Record  a  set  of  images  with  embedded  inhomogeneities  at  different 
illuminated  scanning  positions; 

2.  Generate  a  “clean”  background  image  from  the  “dirty”  background  image 
produced  by  selecting  and  averaging  the  minimally  perturbed  images; 

3.  Calculate  the  difference  images  between  the  recorded  images  and  the 
“clean”  background  image; 

4.  Do  ICA  analysis  of  the  set  of  difference  images  to  generate  the 
Independent  Components  (ICs); 

5.  Obtain  the  contribution  of  each  target  to  intensity  distribution; 

6.  Estimate  3D  location  of  the  target  relative  to  boundaries  by 
numerically  marching  the  target  to  the  surface  until  matching  the 
retrieved  IC. 
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Experimental  Results 

To  illustrate  the  procedures  and  improvement  of  OPTICA  using 
PF  particularly  developed  for  backscattering  geometry,  a  set  of 
scanning  optical  polarized  images  of  a  small  piece  of  cancer¬ 
ous  prostate  tissue  embedded  in  normal  prostate  tissue  were 
acquired  by  the  Photonic  Finger.  Figure  4  shows  the  typical 
scanning  images  chosen  from  the  total  8X8  polarized  images 
recorded  by  scanning  the  incident  beam  on  the  x-y  plane  of  the 
sample.  The  images  shown  on  the  four  comers  [images  (1,  1), 
(1,  8),  (8,  1),  (8,  8)]  are  chosen  by  the  x-y  positions  of  the  inci¬ 
dent  light  far  from  the  embedded  object  while  other  images  in 
the  central  area  [images  (3,  5),  (3,  6),  (4,  5),  (4,  6)]  are  chosen 
by  the  x-y  positions  of  the  incident  light  close  to  the  embedded 
cancerous  tissue.  There  is  no  big  difference  of  the  light  intensity 
distribution  pattern  among  these  recorded  images  except  the 
position  at  the  maximum  intensity.  This  is  because  of  the  sharp 
peak  of  light  intensity  scattered  from  the  surface  and  subsurface 
in  the  backscattered  direction  (16),  which  suppresses  the  pertur¬ 
bation  caused  by  dfi„  and/or  dp'  from  the  embedded  cancerous 
tissue.  To  overcome  this  difficulty,  a  “clean”  background  image 
(CBI)  needs  to  be  synthesized.  In  “clean”  background  image 
synthesis,  the  incident  light  position  of  each  of  8  X  8  images  is 


set  to  origin.  All  array  images  are  then  shifted  to  the  origin  and 
the  size  of  each  image  is  cropped  at  the  boundary  while  incident 
light  distribution  reaches  the  noise  level.  The  images  minimally 
perturbed  by  the  embedded  targets  are  selected  for  synthesizing 
the  “clean”  image. 

After  the  “clean”  background  image  is  obtained,  the  perturba¬ 
tion  2D  (x-y)  images  are  generated  by  extracting  the  “clean” 
background  image  from  the  recorded  images.  In  the  pertur¬ 
bation  2D  images,  the  perturbation  caused  by  the  embedded 
cancerous  tissue  was  highlighted.  ICA  is  then  performed 
upon  the  perturbation  images  to  recognize  leading  indepen¬ 
dent  components  (ICs).  Each  target  is  associated  with  one 
independent  component  (IC).  As  an  example,  Figure  5  shows 
the  OPTICA-generated  independent  intensity  distributions  of 
a  cancerous  prostate  tissue  embedded  in  normal  prostate  tis¬ 
sue.  Figure  5(A)  is  for  the  leading  independent  component, 
and  Figure  5(B)  is  for  residual  (noise)  component.  It  is  clear 
that  the  existence  of  target  (the  cancerous  prostate  tissue) 
can  be  discerned  from  Figure  5(A).  The  x-  and  y-locations  of 
the  cancer  tissue  can  be  obtained  from  the  position  of  maxi¬ 
mum  intensity  of  Figure  5(A).  The  signal  strength  of  Figure 
5(A)  is  much  stronger  than  that  of  Figure  5(B),  indicating 


Figure  4:  A  set  of  64  (8  X  8)  images  recorded  by  scanning  the  incident  beam  on  the  x-y  plane  of  the  sample  (a  small  piece  of  cancerous  prostate  tissue 
hidden  in  the  host  normal  tissue). 
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Figure  5:  OPTIC  A-generated  intensity  distributions  of  the  independent  components  on  the  detector  plane  for  the  sample  consisting  of  a  small  piece  of  can¬ 
cerous  prostate  tissue  embedded  in  normal  prostate  tissues.  (A)  is  the  leading  IC;  and  (B)  is  the  residual  (noise)  component. 


independent  intensity  distribution  in  Figure  5(A)  generated 
by  OPTICA  is  the  leading  IC.  The  OPTICA-generated  inten¬ 
sity  distributions  of  the  leading  independent  component(s) 
on  the  detector  plane  can  be  used  to  locate  the  z-position  of 
the  cancerous  tissue  by  numerical  matching  the  target  to  the 
surface  of  the  medium  until  matching  the  retrieved  IC.  (14). 
The  result  of  z  =  3.1  mm  obtained  from  OPTICA  analysis  is 
in  good  agreement  with  the  actual  depth  of  -3  mm. 

The  OPTICA-generated  independent  components  using  data 
of  the  scanning  fluorescent  images  of  Cytate- stained  cancer¬ 
ous  and  normal  prostate  tissues  embedded  in  large  normal 
prostate  tissue  are  shown  in  Figures  6(A)  to  6(F).  Figure  (A) 
and  (B)  show  the  first  leading  IC  standing  for  Cytate-stained 
cancerous  tissue  location;  (C)  and  (D)  show  the  second  IC 
indicating  Cytate-stained  normal  tissue  location;  (E)  and  (F) 
show  residual  -  noise.  The  3D  locations  of  the  Cytate-stained 
cancerous  prostate  tissues  were  obtained  using  similar  meth¬ 
ods  described  above. 

OPTICA-generated  results  are  summarized  in  Table  II  under 
different  conditions  using  backscattering  polarized  imaging 
and  OPTICA,  which  lists  the  OPTICA-determined  positions 


of  different  objects  in  the  different  experiments  mentioned 
above  in  comparison  with  their  known  3D  locations. 

Another  salient  feature  shown  in  Figure  6  is  the  higher 
emission  intensity  of  the  stained  cancerous  tissue  com¬ 
pared  to  the  stained  normal  tissue.  This  is  attributed  to  the 
preferential  uptake  of  Cytate  by  cancerous  prostate  tissue 
(12,  17-19).  Using  the  data  shown  in  Figure  6,  the  ratio 
of  IJI„  after  OPTICA  analysis  was  found  to  be  ~3.9  and 
~2.3,  for  the  samples  with  depths  of  2.5  mm  and  3.5  mm, 
respectively. 

Discussion 

Prostate  cancer  is  classified  as  an  adenocarcinoma,  or  glan¬ 
dular  cancer  which  is  developed  from  epithelial  cells  (22). 
About  70%  of  prostate  cancer  arises  in  the  peripheral  zone 
(22).  The  mean  mucosa  thickness  of  830  ±  60  pm  and  the 
mean  rectal  wall  thickness  of  2.57  ±  0.15  mm  (22,  23)  allows 
the  NIR  scanning  laser  beam  to  penetrate  the  rectal  wall  to 
reach  prostate  for  optical  imaging  through  rectum  (12,  24). 
By  employing  scanning  NIR  imaging,  Photonic  Finger  can 
be  used  to  characterize  cancerous  prostate  tissue  up  to  ~3  mm 


Table  II 


Comparison  of  the  known  and  OPTICA  determined  positions  of  embedded  objects. 


Object 

Covered  tissue 

Known  ( x ,  y,  z) 
position  (mm) 

OPTICA-generated 
(x,  y,  z )  position  (mm) 

Cancerous  prostate  tissue 

Normal  prostate  tissue 

(13.5,  10.1,  3.0) 

(13.8,  10.3,  3.1) 

Stained  cancerous  tissue 

Normal  prostate  tissue 

(17.8,  14.6,  2.5) 

(17.9,  14.4,  2.3) 

Stained  cancerous  tissue 

Normal  prostate  tissue 

(18,  24,  3.5) 

(17.6,  23.8,  3.6) 

Black  rubber 

Chicken  tissue 

(15.6,  10.4,5.1) 

(15.2,  10.2,  4.7) 

Black  rubber 

Chicken  tissue 

(15.6,  10.4,  6.6) 

(15.3,  10.8,  6.2) 
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Figure  6:  OPTICA-generated  intensity  distributions  of  the  independent  components  on  the  detector  plane  for  the  sample  consisting  of  the  Cytate-stained 
cancerous  and  normal  prostate  tissues  embedded  in  normal  prostate  tissue.  (A)  and  (B)  show  the  Hrst  leading  IC  indicating  for  Cytate-stained  cancerous  prostate 
tissue;  (C)  and  (D)  show  the  second  leading  IC  standing  for  Cytate-stained  normal  prostate  tissue;  and  (E)  and  (F)  show  the  residual  (noise)  component. 


without  contrast  agents.  With  the  enhancement  of  the  recep¬ 
tor-target  contrast  agent,  PF  can  detect  cancer  at  much  deeper 
tissue  layer.  Previous  spectral  polarized  imaging  experiments 


performed  in  IUSL  show  that  ICG-stained  small  object  hid¬ 
den  inside  the  host  prostate  tissues  in  the  rectum-membrane- 
prostate  structures  at  depths  of  7.5  mm  can  be  imaged  and 
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identified  using  fluorescence  imaging  methods  (25).  ICG  is 
the  fluorescent  part  of  Cytate.  It  is  reasonable  to  assume  that 
PF  can  detect  cancer  in  deep  prostate  tissue  layer  more  than 
7.5  mm  under  optimal  condition  with  the  enhancement  of  the 
receptor-target  contrast  agent.  The  PF  will  achieve  diagnos¬ 
tic  function  as  a  screening  tool  other  than  TRUS.  It  is  well 
known  that  TRUS  alone  cannot  be  considered  as  a  screening 
tool  for  prostate  cancer  detection  (9).  It  is  used  as  location 
map  when  numbers  of  biopsy  samples  were  taken  randomly 
from  the  prostate  using  biopsy  needles  of  ~1.2  millimeters  in 
diameter  (25).  The  number  of  biopsies  that  should  be  per¬ 
formed  is  debated.  Twelve  or  eighteen  biopsies  suggested  by 
different  protocols  were  reported  (22).  The  oversampling  and 
overtreatment  of  prostate  cancer  is  now  a  great  concerning 
in  prostate  cancer  research  (22,  25)  not  only  because  of  the 
increased  numbers  of  sampling  will  cause  bleeding  and  pain, 
but  also  because  the  penetration  and  bleeding  may  cause  the 
metastasis  of  prostate  cancer  to  other  organ  (22,  25).  In  con¬ 
trast  to  TRUS,  Photonic  Finger  can  locate  the  3D  positions 
of  cancerous  prostate  sites.  Therefore,  it  may  function  as  a 
detecting  and  screening  instrument.  Photonic  Finger,  an  opti¬ 
cal  scanning  imaging  system  may  have  positive  impact  in 
clinic  applications  to  reduce  significant  pain  and  amount  of 
bleeding  for  patients  in  prostate  cancer  detection. 

Conclusion 

In  conclusion,  a  portable  rectal  NIR  scanning  polarization  imag¬ 
ing  unit  with  an  optical  fiber-based  rectal  probe,  named  Pho¬ 
tonic  Finger,  was  designed,  developed,  assembled  and  tested 
on  ex  vivo  tissues  at  IUSL  in  CCNY.  This  instrument  combined 
with  OPTICA  was  achieved  to  characterize  and  3D-locate 
cancerous  prostate  tissue  embedded  in  normal  prostate  tissue 
in  backscattered  geometry  for  the  first  time.  The  retrieved  3D 
locations  are  in  good  agreement  with  the  known  position  of  the 
embedded  object.  This  technique  was  able  to  sense  weak/small 
absorptive,  scattering,  or  fluorescent  inhomogeneities.  This  new 
developed  Photonic  Finger  instrument  may  provide  an  alterna¬ 
tive  imaging  technique  for  prostate  cancer  detection,  which  is 
accurate,  of  higher  spatial  resolution  and  non-or-less  invasive, 
and  can  be  applied  in  near-real-time  for  in-vivo  detection. 
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Abstract:  The  decreasing  relative  content  of  tryptophan,  NADH  and  flavin  of  the  retinoic  acid-treated  breast 
cancerous  cells  were  demonstrated  using  native  fluorescence  spectra.  The  changes  may  be  used  to  evaluate 
the  chemotherapeutic  effects  on  cancer. 

OCIS  codes:  170.0170,  170.4580,  170.6280 

1.  Introduction 

Optical  spectroscopy  has  been  considered  as  an  alternative  technique  for  cancer  detection  since  Alfano  and 
coworkers  were  the  first  to  measure  native  fluorescence  spectra  of  malignant  breast  tissues  [1].  The  aim  of  this 
present  research  is  to  determine  if  the  native  fluorescence  spectroscopy  is  effective  enough  to  detect  changes  of 
fluorophore  compositions  related  to  the  treatment  of  chemotherapy  on  malignant  cell  lines  using  Nonnegative  Matrix 
Factorization  (NMF)  analysis. 

2.  Samples  and  Methods 

The  human  breast  cells  were  purchased  from  American  Type  of  Culture  Collection  (ATCC),  Rockville,  MD.  The  cell 
lines  used  in  the  study  were:  malignant  human  breast  cell  lines  ATCC  FITB22  (Adenocarinoma  Pleural  Effusion)  and 
ATCC  FITB126  (Dutal  Carcinoma).  The  cells  were  cultured  with  90%  modified  Dulbeco’s  minimal  essential  medium 
(Eagle)  (GIBO)  and  10%  of  the  fetal  bovine  serum  (GIBCO)  in  a  humidified  atmosphere  of  4%  CO2  at  37°C.  Half  of 
the  cells  were  treated  with  ~10"6  M  retinoic  acid  and  the  remainder  served  as  controls  while  the  cell  lines  are  cultured. 
The  measured  cell  concentration  was  estimated  as  —105  -  110  x  106 cells/ml. 

The  fluorescence  emission  spectra  of  cells  were  measured  using  Perkin-Elmer  LS  50  spectrometer.  The  excitation 
wavelengths  of  300nm  and  340  nm  were  selected  based  on  our  previous  spectral  study  of  the  main  fluorophores  in 
breast  cells  [2].  Combined  with  NMF,  the  native  fluorescence  measurements  were  used  to  distinguish  the  RA-treated 
breast  cancer  cells  from  the  untreated  cell.  The  main  advantage  of  NMF  is  that  the  NMF-recovered  spectral  data  and 
concentrations  of  constituents  are  positive  values.  This  makes  the  calculated  results  represent  more  physically  and/or 
biologically  meaning  because  the  spectra  and  contents  of  the  fluorephores  should  have  non-negative  values. 

3.  Experimental  results  and  discussion 

Fourteen  pairs  of  the  RA-treated  and  untreated  breast  cancer  cell  samples  were  investigated  using  fluorescence 
spectroscopy  with  excitation  wavelength  of  300  nm.  The  average  fluorescence  spectral  profiles  with  standard 
deviation  error  bars  at  key  wavelengths  for  the  RA-untreated  and  treated  cancerous  breast  cells  are  shown  in  Figs. 
1(a)  and  1(b),  respectively.  Each  spectral  profile  was  normalized  to  unit  value  of  1  (i.e.,  the  sum  of  squares  of  the 
elements  in  each  emission  spectra  data  was  set  as  1)  before  taking  average  and  calculation. 


Wavelength(nm)  Wavelength(nm) 


Fig.  1  Average  fluorescence  spectra  of  (a)  untreated  and  (b)  Retinoic  Acid  (RA)-treated  cancerous 
breast  cells  obtained  with  the  excitation  of  300  nm. 
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The  main  emission  peaks  of  both  untreated  and  treated  cancerous  breast  cells  were  found  at  337nm.  The  major 
difference  of  the  fluorescence  profiles  between  the  untreated  and  treated  cells  is  that  there  exists  a  much  lower  local 
peak  at  ~  460nm  for  the  untreated  cells  while  no  such  a  peak  for  the  treated  cells.  A  tiny  local  peak  at  529  nm  may 
be  observed  in  the  emission  of  the  untreated  cells,  which  can  be  more  clearly  seen  by  enlarging  the  fluorescence 
profile  at  the  spectral  range  of  700  nm  to  900  nm  as  shown  in  insert  of  Fig.  1  (a). 

To  investigate  the  relative  content  changes  of  Principal  Components  (PCs)  and  compare  with  their  corresponding 
fluorophores,  NMF  was  used  to  extract  spectra  of  PCs  (fluorophores)  of  the  1st  component  (tryptophan),  the  2nd 
(NADH)  and  the  3ld  (flavin)  for  the  cells,  which  are  shown  as  dash,  dot  and  dash-dot  lines  in  Fig.  2(a),  respectively. 
The  measured  spectra  of  individual  tryptophan,  NADH  and  flavin  in  solutions  are  plotted  as  thicker,  thick  and  thin 
solid  lines,  respectively,  as  references.  The  two  groups  of  spectra  (the  extracted  and  the  measured)  for  each 
fluorophore  show  reasonable  agreement,  which  demonstrates  that  NMF  model  accounts  for  the  major  spectroscopic 
feature  observed,  and  indicates  that  the  measurements  are  reasonable  [2]. 
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Fig.  2  (a)  Comparison  of  the  extracted  spectra  of  the  three  principal  components  (PCs)  of  the  lsl  PC  -  tryptophan  (dash  line),  2nd  PC  - 
NADH  (dot  line)  and  3  rd  PC  -  flavin  (dash-dot  line)  for  the  untreated  and  treated  breast  cancerous  cells,  and  the  measured  spectra  of 
individual  tryptophan  (thicker  line),  NADH  (thick  line)  and  flavin  (thin  line)  in  solution;  (b)  Fractional  content  of  the  1st  PC  - 
tryptophan  vs.  that  of  the  2nd  PC  -  NADH;  (c)  Fractional  content  of  the  3rd  PC  -  flavin  vs.  that  of  the  1st  PC  -  tryptophan;  (d)  Fractional 
content  of  the  3rd  PC  -  flavin  vs.  that  of  the  2nd  PC  -  NADH  by  analyzing  emission  spectra  excited  at  300  nm  using  NMF  method.  The 
data  for  the  untreated  and  treated  cells  are  displayed  with  the  solid  circles  and  solid  squares  in  2  (b)  -  2  (d),  respectively. 


To  investigate  the  relative  content  changes  of  PCs,  the  contents  of  tryptophan,  NADH  and  flavin  in  two  types  of 
cells  were  extracted  from  the  measured  total  fluorescence  spectra  using  the  NMF  analysis.  Fig.  2(b)  shows  the 
relative  content  of  the  1st  PC  (tryptophan)  vs.  the  2nd  PC  (NADH);  2(c)  displays  the  3rd  PC  (flavin)  vs.  the  1st  PC 
(tryptophan)  and  2(d)  exhibits  the  3rd  PC  (flavin)  vs.  the  2nd  PC  (NADH)  of  the  untreated  (solid  circles)  and  treated 
(solid  squares)  cancerous  breast  cells.  The  most  salient  feature  of  Fig.  2(b)  is  that  all  data  points  for  the  untreated 
cells  are  located  on  the  upper  side  over  the  data  points  for  the  treated  cells,  indicating  that  the  relative  contribution  of 
tryptophan  is  higher  in  the  untreated  cancerous  cells  in  comparison  with  the  treated.  Fig.  2(c)  shows  again  that  the 
relative  fluorescence  contribution  of  tryptophan  is  higher  in  the  untreated  cancerous  cells  in  comparison  with  the 
treated.  Fig.  2(d)  shows  that  most  data  points  for  the  untreated  cells  are  located  in  the  up-right  side  in  comparison 
with  the  data  points  for  the  treated  cells,  indicating  that  the  relative  contents  of  NADH  and  flavin  are  higher  in  the 
untreated  cancerous  cells  in  comparison  with  the  treated. 

The  following  biomedical  studies  may  help  to  understand  these  changes.  Brown  et  al  studied  chemotherapy 
effects  for  breast  cancer  [3].  After  chemotherapy,  the  decreased  plasma  amounts  of  tryptophan  in  plasma  were 
observed  [3]  since  increased  tryptophan  levels  are  somehow  related  to  increased  tumor  cell  proliferation  [3].  NADH 
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is  one  of  the  most  important  coenzyme  for  the  production  of  Adenosine-Tri-Phosphate  (ATP)  [4].  The  results  of 
decrease  of  NADH  in  the  RA-treated  cancerous  cells  may  reflect  a  reduction  in  the  production  of  ATP  in  these 
treated  cells.  Since  the  decrease  of  tryptophan  and  NADH  somehow  reflect  the  less  activity  of  cells  [3,  4],  it  can  be 
hypothesized  that  the  RA-treatment  results  in  cancerous  cell  death  or  reduction  of  the  activity  of  cancer  cells. 

Emission  spectra  of  nine  pairs  of  the  RA-treated  and  untreated  breast  cancer  cell  samples  were  measured  with  the 
excitation  of  340  nm.  The  average  fluorescence  spectral  profiles  of  the  untreated  (solid  line)  and  treated  (dash  line) 
cells  are  displayed  in  Figs.  3(a),  which  shows  that  the  emission  profiles  of  the  untreated  and  treated  cells  were 
mainly  contributed  from  NADH  corresponding  to  a  peak  at  ~  459nm,  and  flavin  corresponding  to  a  peak  of  ~ 
530nm.  The  major  difference  of  the  profiles  between  the  untreated  and  treated  cells  is  that  the  shoulder  peak  at  ~ 
530  nm  for  the  treated  cells  is  little  stronger  than  that  for  the  untreated  cells.  Fig.  3  (b)  displays  the  relative  content 
of  the  flavin  vs.  NADH  of  the  untreated  (solid  circle)  and  treated  (solid  square)  cancerous  breast  cells  analyzed 
using  the  NMF  analysis.  The  obvious  feature  of  Fig.  3(b)  is  that  the  relative  contribution  of  NADH  is  higher  in  the 
untreated  cancerous  cells  in  comparison  with  the  treated,  which  also  reflects  the  decrease  of  the  relative  content  of 
NADH  in  the  treated  cancerous  breast  cells  in  comparison  with  the  untreated  cells. 
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Fig.  3  (a)  Average  fluorescence  spectra  of  the  untreated  (solid)  and  RA-treated  (dash)  cancerous  breast  cells  obtained  with  the 
excitation  wavelength  of  340  nm.(b)  Relative  content  of  the  NADH  vs.  flavin,  obtained  from  the  measured  fluorescence 
spectroscopy  of  the  breast  cells  with  340  nm  excitation  using  the  NMF  analysis. 


4.  Conclusions 

This  study  reports  on  detecting  the  chemotherapeutic  effects  on  malignant  breast  cells  by  measuring  the  nature 
fluorescence  spectra  of  cells  and  extracting  the  fluorescence  spectral  features  of  principal  components  using  an  NMF 
method.  The  presence  and  contents  of  the  detectable  fluorophores  in  breast  malignant  cells  such  as  tryptophan, 
NADH  and  flavin  were  analyzed  as  the  Principal  Components  for  the  NMF.  We  demonstrated  that  fluorescence 
intensities  of  tryptophan,  NADH  and  flavin  decreased  in  the  RA-treated  cells  in  comparison  with  those  in  the 
untreated  cells.  This  work  shows  the  change  of  relative  contents  of  tryptophan,  NADH  and  flavin  studied  using 
native  fluorescence  spectroscopy  with  the  NMF  analysis  may  present  potential  criteria  for  evaluation  of 
chemotherapeutic  effects  on  malignant  breast  cells. 
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ABSTRACT 

Time-resolved  fluorescence  polarization  spectroscopy  of  various  molecular  volume  fluorescein  dye-labeled  polymers 
influenced  by  the  molecular  rotation  motion  was  investigated  by  picosecond  dynamics  measurements.  The  time-resolved 
picosecond  polarization  spectroscopy  was  modeled  using  a  system  of  first-order  linear  differential  equations  containing 
two  main  parameters:  the  decay  rate  of  emission  and  the  rate  of  one  orthogonal  emission  component  transferring  to 
another.  Both  experimental  observation  and  theoretical  calculation  indicate  that  the  fluorescent  molecular  rotation  is  the 
cause  of  the  depolarization.  The  dipole’s  rotational  time  was  extracted  and  shows  increasing  with  the  molecular  volume 
obeying  the  Einstein-Stokes  relation. 

Keywords:  Time-resolved,  picosecond  spectroscopy,  fluorescence  depolarization,  fluorescence  anisotropy,  molecular 
orientation  and  rotation,  fluorescence  lifetime,  rotation  time 


1.  INTRODUCTION 

The  depolarization  of  fluorescent  molecule  in  solution  was  first  realized  by  A.  Jablonski  [1],  This  fluorescence 
depolarization  technique  can  be  used  as  a  probe  in  chemistry  and  biology  to  determine  the  rotational  rate  of  molecules  by 
measuring  the  decay  of  the  different  polarized  components  [2,  3], 

In  this  study,  time-resolved  fluorescence  spectroscopy  of  fluorescein  dye-labeled  polymer  chains  with  molecular  volume 
ranged  from  -410  angstrom3  to  -5000  angstrom3  at  two  polarization  directions,  parallel  and  perpendicular  to  that  of  the 
linearly  polarized  exciting  light,  were  measured.  It  was  observed  in  all  measurements  that  the  decay  rate  of  the 
fluorescence  component  polarized  parallel  to  the  exciting  beam  is  faster  than  that  of  the  fluorescence  component 
perpendicular  to  the  excitation.  A  theoretical  fluorescence  depolarization  model  developed  by  G.  Weber  [4]  was 
extended  from  micro-second  to  pico-second  regimes  to  investigate  and  demonstrate  the  origin  of  the  different  decay  rate 
of  the  two  fluorescence  components.  This  model  describe  the  time-resolved  fluorescence  dynamics  of  fluorescein- 
polymer  using  a  set  of  the  linear  differential  equations  containing  two  main  parameters:  the  decay  rate  of  emission  and 
the  rate  of  one  orthogonal  emission  component  transferring  to  another.  Experimental  tests  of  this  theory  were  verified  by 
measuring  the  time-resolved  parallel  and  perpendicular  components  of  the  fluorescence  using  a  streak  camera  system. 
Both  experimental  observation  and  theoretical  calculation  indicate  that  the  fluorescent  molecular  rotation  is  the  cause  of 
the  depolarization.  The  dipole’s  rotational  time  was  extracted  and  shows  increasing  with  the  molecular  volume  obeying 
the  Einstein-Stokes  relation. 


2.  EXPERIMENTAL  METHOD 


2.1  Samples  and  setup 

The  dyes  used  in  the  time-resolved  fluorescence  polarization  measurements  are  fluorescein  and  fluorescein  conjugated 
dextrans  with  different  molecular  volume.  Fluorescein  is  a  synthetic  organic  compound  available  as  a  dark  orange/red 
powder  soluble  in  water  and  alcohol.  It  is  widely  used  as  a  fluorescent  tracer  for  many  applications  including 
biochemical  research  and  health  care  applications.  The  chemical  structure  of  Fluorescein  is  shown  in  Fig.l.  The 
fluorescence  of  this  molecule  is  very  high,  and  excitation  occurs  at  494  nm  and  emission  at  521  nm.  In  this  study, 
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fluoresceins  were  provided  by  Malinkrodt  Inc  (St.  Louis,  MO).  The  dimensions  of  fluorescein  molecules  were  estimated 
as  oblate  with  semi-axes  of  7  angstrom  and  2 
angstrom  for  bond  length  and  van  der  Waals  radii 
[5],  respectively.  They  were  also  conjugated  to 
different  dextrans  with  molecular  size  range  from 
IK  angstrom  to  700K  angstrom  by  Malinkrodt  Inc. 

The  molecular  weight  and  length  per  unit  dextran  is 
—162  dalton  and  ~3  angstrom,  respectively  [5,  6]. 

The  concentration  of  aqueous  solution  used  for  the 
fluorescence  depolarization  measurements  was  in 
the  range  of~50|iM  for  fluorescein  derivatives  at 
room  temperature. 

The  experimental  arrangement  of  time -resolved 
fluorescence  polarization  measurements  has  been 
given  elsewhere  [7].  Pulses  of  100  fs  at  400  nm 
obtained  from  the  second  harmonic  of  a  Ti:Sapphire 
laser  were  used  to  pump  samples.  The  temporal 
profdes  of  emission  from  the  dye-polymer 
conjugates  of  I||(t)  and  l_(t)  corresponding  to  the 
linearly  polarized  excitation  light  were  recorded  by  a 
synchroscan  streak  camera.  The  fluorescence  was 
collected  by  a  large  diameter  lens  (L2)  with  a  focal  length  of  5  cm  and  directed  onto  a  synchroscan  streak  camera  with 
temporal  resolution  of  10  ps.  A  long  pass  filter  (LP)  was  used  to  block  the  illuminating  light  so  that  only  the  emission 
from  the  samples  was  detected.  Polarizer  Pi  was  used  to  ensure  the  linear  polarization  of  the  input  laser  pulses  and  P2  was 
used  as  an  analyzer.  The  polarization  of  Pj  was  rotated  between  0°  and  90°  with  respect  to  that  of  P2  to  record  the 
intensity  profiles  of  the  parallel  and  perpendicular  polarization  components  of  the  fluorescence.  The  temporal  profiles 
recorded  by  a  silicon  intensified  target  (SIT)  of  the  streak  camera  were  analyzed  to  obtain  temporal  and  polarization 
information. 


Fig.  1.  Molecular  structure  of  fluorescein.  The  dimensions  of 
fluorescein  molecules  were  estimated  as  oblate  with  semi-axes  as 
7  angstrom  and  2  angstrom  from  bond  length  and  van  der  Waals 
radii  by  its  structure. 


2.2  Analysis  model 

The  time-resolved  fluorescence  depolarization  technique  is  widely  used  as  a  probe  for  re-orientation  motion  of  molecules 
in  solution  since  1970’s  [2-4].  A  theory  developed  by  G.  Weber  et  al  described  the  equalization  of  the  polarized  intensity 
among  the  directions  in  the  nanosecond  scale  [4]  under  specified  conditions  of  temperature,  viscosity  and  ordinary  rate 
constants.  We  extend  this  model  to  the  pico-second  regimes  to  explore  the  effects  of  the  molecular  rotation  to  the 
fluorescence  polarization  spectroscopy. 


For  a  single  type  of  fluorephore,  temporal  profiles  of  the  parallel  and  perpendicular  components:  Ij_(t)  and  I||(t)  can  be 
expressed  as  [8]: 


^11  (0  =  ~exp( — — ) 
3  r. 


1  +  2 r0  exp( - ) 

r  , 
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3  r 
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rot  J 


where  /«  is  the  initial  emission  intensity,  Zj 4s  the  fluorephore’s  lifetime,  r0  is  anisotropy  value  at  t=0  indicating  the  initial 
distribution  of  the  excited  polarized  dipoles,  and  zra,  is  re-orientation  time  indicating  the  depolarization  rate  of  the  excited 
polarized  molecules,  which  is  determined  by  the  microenvironments  of  the  medium.  In  a  simple  case  where  the 
fluorephore  molecule  undergoes  Brownian  rotation  as  an  Einstein  sphere,  the  re-orientation  time  zrot  can  be  expressed  in 
terms  of  the  medium  viscosity  (q),  the  absolute  temperature  of  medium  (T)  and  the  molecular  volume  (V)  by  [2,  3]: 

_  TjV 

Tm  ~  kT 

where  k  is  the  Boltzmann  constant.  It  can  be  seen  from  Eequation  (1)  that  measured  temporal  profiles  of  the  parallel  and 
perpendicular  components  of  emission  should  include  the  information  of  the  fluorephore’s  decay  time  (xf)  and  the  re- 
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orientation  (rotational)  time  (Trot).  To  investigate  the  emission  intensity  and  life  time  of  fluorephore  if,  the  total  temporal 
emission  intensity  can  be  calculated  using  an  equation: 


1  total  (0  =  T||  ( t )  +  2 1 L  (t)  =  70  exp( — -) 

Tf 


(3) 


Equation  (3)  indicates  that  the  life  time  zyof  the  excited  fluorephores  can  be  obtained  from  a  temporal  profile  of  the  total 
emission  (Itotai(t)). 

The  time-resolved  polarization  anisotropy  of  emission  is  defined  as  [2,  3]: 


/|,  (0-/l(0 
711(0  +  27^(0 


=  'o  exp(- 


^ rot 
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where  I||(t)  and  1 1  (t)  are  measured  intensities  of  the  parallel  and  perpendicular  polarization  components  of  the  time- 
resolved  fluorescence,  respectively.  Equation  (4)  indicates  that  the  re-orientation  time  (xrot)  can  be  extracted  from  the 
time-resolved  polarization  anisotropy  (r(t))  of  the  emission. 

Using  equations  (3)  and  (4),  we  can  separately  study  the  life  time  and  rotational  time  of  the  emission  of  the  fluorephores 
from  the  time-resolved  total  emission  and  polarization  anisotropy,  respectively. 


3.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

Time-resolved  fluorescence  intensity  Time-resolved  fluorescence  intensity 


onnnn 


Time-resolved  polarization  anisotropy  Time-resolved  polarization  anisotropy 


Fig.  2.  The  time  dependent  fluorescence  intensity  profiles  of  parallel  (solid  line)  and  perpendicular  (dot  line)  components  of 
emission  from  conjugated  fluorescein-polymer  of  MV  at  (a)  ~707  and  (b)  ~3,200  angstrom3.  The  temporal  profiles  of 
the  polarization  anisotropy  of  conjugated  fluorescein-polymer  of  MV  at  (c)  ~707  and  (d)  ~3,200  angstrom3. 
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The  time-resolved  fluorescence  polarization  profiles  were  recorded  for  fluorescein  dye,  and  the  fluorescein  dye-labeled 
dextrans  with  different  Molecular  Volume  (MV)  from  400  to  5,000  angstrom3.  The  dimensions  of  fluorescein  were 
estimated  as  oblate  with  semi-axes  of  7  and  2  angstrom  [5].  As  an  example,  Figs.  2  (a)  and  2  (b)  show  the  measured 
temporal  profiles  of  parallel  and  perpendicular  components  of  the  emission  from  the  conjugated  dye -polymer  with 
different  MV  of  ~707  and  ~3,200  angstrom3,  respectively.  The  solid  line  profiles  show  the  parallel  component,  and  the 
dot  lines  show  the  perpendicular  one.  Fig.  2  (a)  and  2  (b)  show  that  I|j(t)  is  greater  than  1 1  (t)  at  all  of  the  decay  times 
indicating  that  this  dye -polymer  keeps  property  of  the  polarization  preservation.  The  time-resolved  polarization 
anisotropy  r(t)  can  be  obtained  by  substituting  the  data  of  I||(t)  and  1 ,  (t)  shown  in  Figs.  2  (a)  and  2  (b)  into  the  time 
dependent  equation  (4).  The  calculated  results  are  shown  in  Figs.  2  (c)  and  2  (d)  displaying  the  temporal  profiles  of 
polarization  anisotropy  for  Fluorescein-polymer  with  different  MV  of  ~707  and  ~3,200  angstrom3,  respectively. 


Comparison  of  experimental  and  calculated  rotation  time 


Fig.  3.  (a)  The  rotation  time  obtained  by  theoretical  calculation  (dot  line)  and  by  fitting  the  experimental  data 
(solid  line)  of  the  fluorescein  dye-polymer  conjugate  as  a  function  of  molecular  volume. 

The  salient  features  are:  (1)  the  temporal  polarization  anisotropy  shown  in  Fig.  2  (c)  exhibits  a  faster  decay  than  the 
profile  shown  in  Fig.  2  (d),  indicating  that  smaller  molecular  volume  gives  faster  Brownian  rotation  as  described  by 
Equation  (2);  (2)  The  initial  anisotropy  r0  shown  in  Fig.  2  (c)  is  greater  than  that  shown  in  Fig.  2  (d).  In  order  to  study 
rotation  time  changes  with  the  molecule  volume,  the  experiments  were  performed  for  a  set  of  Fluorescein-polymer 
conjugates  with  attached  chain  having  the  volume  varying  from  ~400  to  5,000  angstrom’.  The  rotation  times  can  be 
obtained  by  fitting  the  experimental  data  of  r(t)  with  Equation  (4)  and  the  results  are  displayed  as  solid  lines  in  Fig.  3. 
The  theoretical  values  of  the  rotation  times  with  different  attached  polymer  length  were  calculated  using  Equation  (2) 
and  compared  with  the  experimental  values,  shown  as  the  dot-line  profile  in  Fig.  3.  Since  the  fluorescein  dye  appears  as 
approximated  oblate  symmetric  rotors  [5],  the  volume  were  calculated  to  be  ~  410  angstrom3  as  oblate  ellipsoid  with 
semi-axes  7  and  2  angstrom  while  no  polymer  attached. 
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The  salient  features  of  Fig.  3  are:  (1)  the  rotation  times  raise  with  the  volume  of  the  attached  polymer  chain,  which  are  in 
good  agreement  with  experimental  values  and  exhibit  a  linear  dependence  on  the  volume  of  molecules  as  predicted  by 
Equation  (2)  when  the  volume  of  conjugated  Fluorescein-polymers  is  less  than  3,200  angstrom3.  This  can  be  understood 
because  the  Brownian  rotational  motion  for  larger  volume  conjugated  dye -polymer  is  slower  given  rise  to  a  longer 
rotation  time;  (2)  The  rotation  time  becomes  flat  when  the  volume  of  the  dye-polymer  increases  beyond  ~3,000 
angstrom3.  This  indicates  that  for  large  size  molecules,  where  the  volume  exceeds  a  critical  value,  the  rotation  time 
changes  very  slow  with  volume  of  the  fluorescent  molecules.  The  incorrect  prediction  of  Equation  (2)  may  be  because 
the  Fluorescein-polymer  can  no  longer  be  estimated  as  a  simple  dimension  as  a  sphere  or  an  oblate  since  it  is  difficult  to 
control  the  polymer  chain  extending  direction.  The  irregular  shapes  of  molecules  result  in  the  different  rotation  rate  of  Z),y 
[2-4],  and  therefore  the  depolarization  will  not  show  simple  exponential  decay  as  described  by  time-resolved  polarization 
anisotropy  r(t)  in  form  of  Equation  (3).  In  addition,  the  rapid  internal  motions  originated  by  flexible  structure  of  the 
conjugated  fluorescein-polymers  cause  the  “extra”  depolarization  other  than  rotation  of  the  molecule  [2]. 

4.  CONCLUSION 

We  have  experimentally  and  theoretically  investigated  that  when  an  isotropic  fluorescent  molecular  system  is  excited 
with  a  linearly  polarized  light  beam,  the  intensity  decay  of  their  fluorescence  component  polarized  parallel  to  the  exciting 
beam  is  faster  than  that  of  the  fluorescence  component  perpendicular  to  the  exciting  beam.  This  phenomenon  can  be 
explained  with  the  time-resolve  fluorescence  depolarization  caused  by  the  Brownian  rotation  of  the  dipoles.  For  a  simple 
and  stable  molecule,  the  effect  of  Brownian  rotation  is  the  main  cause  of  the  depolarization.  Time-resolved 
depolarization  measurements  were  used  to  obtain  the  rotation  time  and  initial  polarization  anisotropy.  We  have  also 
shown  that  over  a  the  range  of  the  molecular  volume  of  dye-polymer  when  r(t)  exhibit  the  single  exponential  decay,  the 
accuracy  of  the  theoretical  calculation  is  quite  comparable  with  value  obtained  by  fitting  the  experimental  value. 
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